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I 

Introduction 


Dark matter (DM) ^ has not yet been observed in particle physics 
experiments, and there is not yet any evidence for non-gravitational 
interactions between Dark Matter and Standard Model (SM) par¬ 
ticles. If such inferacfions exisf, parficles of Dark Maffer could be 
produced af fhe LHC. Since Dark Maffer parficles fhemselves do 
nof produce signals in fhe LHC defecfors, one way fo observe fhem 
is when fhey are produced in association with a visible SM particle 
X{=g, q, 7, Z, W, or h). Such reactions, which are observed at collid¬ 
ers as particles or jets recoiling against an invisible state, are called 
"mono-X" or tj+X reactions (see e.g Refs. [BMP04; FST06; PQZ08; 
Bel+10; BFHio]), where is fhe missing fransverse momenfum 
observable in fhe defector. 

Early Tevafron and LHC Run-i searches for Bj+X signafures af 
CDF [Aal+12], ATLAS [ATLi5d; ATL15C; ATL14C; ATLi4b; ATLi4a; 
ATLi5b; ATLi5a; ATLi4d] and CMS [CMSi5b; CMSi4b; CMSi5e; 
CMSi5d; CMSi5f; CMS14C; CMSi5a], employed a basis of confacf 
inferacfion operators in effecfive field fheories (EFTs) [Goo+11; 
Goo+10] fo calculate fhe possible signals. These EFTs assume fhaf 
producfion of Dark Maffer fakes place fhrough a confacf interacfion 
involving a quark-antiquark pair, or two gluons, and two Dark 
Matter particles. In this case, the missing energy distribution of 
fhe signal is determined by fhe nafure and fhe mass of fhe Dark 
Maffer particles and the Lorentz structure of fhe interacfion. Only 
the overall production rate is a free paramefer fo be consfrained 
or measured. Provided fhat fhe confacf inferacfion approximafion 
holds, fhese EFTs provide a sfraighfforward way fo compare fhe 
resulfs from differenf collider searches wifh non-collider searches 
for Dark Maffer. 

The EFT describes fhe case when fhe mediator of fhe inferacfion 
befween SM and DM parficles are very heavy; if fhis is nof fhe case, 
models fhaf explicifly include fhese mediafors are needed [Goo+11; 
SV12; BFHio; Kopii; Fox+11; Fox+12; SV12; Bus+iqa]. Some "sim¬ 
plified models" [AST09; GSii; Alv+12] of Dark Maffer producfion 
were consfrucfed, including parficles and inferacfions beyond fhe 
SM. These models can be used consisfently af LHC energies, and 
provide an exfension fo fhe EFT approach. Many proposals for such 
models have emerged (see, for example Refs. [AJW12; AHW13; 


’ Many theories of physics beyond the 
Standard Model predict the existence 
of stable, neutral, weakly-interacting 
and massive particles that are putative 
Dark Matter candidates. In the follow¬ 
ing, we refer to such matter as Dark 
Matter, even though the observation 
of such matter at a collider could only 
establish that it is neutral, weakly- 
interactive, massive and stable on the 
distance-scales of tens of meters. 
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DiF+13; BDM14; BB13; BB14; AWZ14; Abd+14; Mal+14; Har+15; 

BFG15; HR15; BT13; Car+13; Bel+12; PS14; Car+14]). At the LHC, 
the kinematics of mono-X reactions occurring via a TeV-scale me¬ 
diator can differ subsfanfially from fhe prediction of fhe confacf 
inferacfion. The mediator may also produce qualifafively differenf 
signals, such as decays back info Sfandard Model particles. Thus, 
appropriate simplified models are an imporfanf componenf of fhe 
design, optimization, and inferprefafion of Dark Maffer searches af 
ATLAS and CMS. This has already been recognized in fhe CDF, AT¬ 
LAS and CMS searches quofed above, where bofh EFT and selected 
simplified model resulfs are presented. 

1.1 The ATLAS/CMS Dark Matter Forum 

To undersfand whaf signal models should be considered for fhe 
upcoming LHC Run-2, groups of experimenfers from bofh AT¬ 
LAS and CMS collaborations have held separate meetings with 
small groups of fheorisfs, and discussed furfher af fhe DM@LHC 
workshop [Mal+14; Abd+14; Abd+15]. These discussions idenfified 
overlapping sefs of simplified models as possible benchmarks for 
early LHC Run-2 searches. Following fhe DM@LHC workshop, 

ATLAS and CMS organized a forum, called fhe ATLAS-CMS Dark 
Matter Forum, fo form a consensus on fhe use of fhese simplified 
models and EFTs for early Run-2 searches wifh fhe parficipafion 
of experfs on fheories of Dark Maffer. This is fhe final reporf of fhe 
ATLAS-CMS Dark Maffer Forum. 

One of fhe guiding principles of fhis reporf is fo channel fhe 
efforfs of fhe ATLAS and CMS collaborations towards a minimal 
basis of dark maffer models fhaf should influence fhe design of 
fhe early Run-2 searches. Af fhe same fime, a fhorough survey of 
realisfic collider signals of Dark Maffer is a crucial inpuf fo fhe 
overall design of fhe search program. 

The goal of fhis reporf is such a survey, fhough confined wifhin 
some broad assumptions and focused on benchmarks for kinematically- 
disfincf signals which are mosf urgently needed. As far as fime and 
resources have allowed, the assumptions have been carefully mo¬ 
tivated by theoretical consensus and comparisons of simulations. 

But, to achieve such a consensus in only a few monfhs before fhe 
sfarf of Run-2, if was imporfanf fo resfricf fhe scope and timescale 
fo fhe following: 

1. The forum should propose a prioritized, compact set of bench¬ 
mark simplified models fhaf should be agreed upon by bofh 
collaborafions for Run-2 searches. The values for fhe scan on 
fhe paramefers of fhe models for which experimenfal resulfs are 
provided should be specified, fo facilifafe fheory reinferprefafion 
beyond fhe necessary model-independenf limifs fhaf should be 
provided by all LHC Dark Maffer searches. 

2. The forum should recommend fhe use of fhe sfate of fhe arf cal- 
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culations for these benchmark models. Such a recommendation 
will aid the standardization the event generator implementation 
of the simplified models and the harmonization of other com¬ 
mon technical details as far as practical for early Run-2 LHC 
analyses. It would be desirable to have a common choice of lead¬ 
ing order (LO) and next-to-leading order (NLO) matrix elements 
corresponding to the state of the art calculations, parton shower 
(PS) matching and merging, factorization and renormalization 
scales for each of the simplified models. This will also lead to 
a common set of theory uncertainties, which will facilitate the 
comparison of results between the two collaborations. 

3. The forum should discuss how to apply the EFT formalism and 
present the results of EFT interpretations. 

4. The forum should prepare a report summarizing these items, 
suitable both as a reference for the internal ATLAS and CMS au¬ 
diences and as an explanation of early Run-2 LHC benchmark 
models for theory and non-collider readers. This report repre¬ 
sents the views of its endorsers, as participants of the forum. 

1.2 Grounding Assumptions 

We assume that interactions exist between Standard Model hadrons 
and the particles that constitute cosmological Dark Matter. If this is 
not the case, then proton collisions will not directly produce Dark 
Matter particles, and Dark Matter will not scatter off nuclei in direct 
detection experiments. 

The Dark Matter itself is assumed to be a single particle, a Dirac 
fermion WIMP, stable on collider timescales and non-interacting 
with the detector. The former assumption is reductionistic. The rich 
particle content of the Standard Model is circumstantial evidence 
that the Dark Matter sector, which constitutes five times as much 
of the mass of the universe, may be more complex than a single 
particle or a single interaction. But, as was often the case in the 
discoveries of the SM, here only one mediator and one search chan¬ 
nel might play a dominant role in the opening stages of an LHC 
discovery. The latter assumption focuses our work on early LHC 
searches, where small kinematic differences between models will 
not matter in a discovery scenario, and with the imminent re-start 
of the LHC our report relies heavily on a large body of existing 
theoretical work which assumed Dirac fermionic Dark Matter. 

Different spins of Dark Matter particles will typically give sim¬ 
ilar results. Exceptions exist: For example, the choice of Majo- 
rana fermions forbids some processes that are allowed for Dirac 
fermions [Goo+ii]. Aside from these, adjusting the choice of Dirac 
or Majorana fermions or scalars will produce only minor changes 
in the kinematic distributions of the visible particle and is expected 

to have little effect on cut-and-count^ analysis. Thus the choice of Cut-and-count refers to an analysis 

Dirac fermion Dark Matter should be sufficient as benchmarks for ^ selection 

and checks the inclusive number 

of events which pass. This is to be 
contrasted with a shape analysis, 
which compares the distribution of 
events. 
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the upcoming Run-2 searches. 

One advantage of collider experiments lies in their ability to 
study and possibly characterize the mediator. A discovery of an 
anomalous B j signafure af fhe LHC would nof uniquely imply 
discovery of dark maffer, while af fhe same fime e.g. discovery of 
an anomalous and annually-modulafed signal in a direcf-defecfion 
experimenf would leave unanswered many quesfions abouf fhe 
nafure of fhe inferacfion fhaf could be resolved by fhe simulfaneous 
discovery of a new mediator parficle. Collider, direcf, and indirecf 
defection searches provide complemenfary ways to approach fhis 
problem [Bau+13], and if is in fhis spirif fhaf much of our focus is 
on fhe mediator. 

We sysfemafically explore fhe basic possibilities for mediators 
of various possible spins and couplings. All models considered are 
assumed fo produce a signafure wifh pairs of Dark Maffer parficles. 
Though more varied and interesting possibilities are added fo fhe 
liferafure almosf daily, these basic building blocks account for much 
of fhe physics sfudied af hadron colliders in fhe pasf fhree decades. 

We also assume fhaf Minimal Flavor Violation (MFV) [CG87; 
HR90; Bur+01; D'A+02] applies fo fhe models included in fhis re¬ 
port. This means that the flavor sfrucfure of fhe couplings between 
Dark Maffer and ordinary parficles follows fhe same sfrucfure as 
fhe Sfandard Model. This choice is simple, since no additional fhe- 
ory of flavor is required, beyond whaf is already presenf in fhe SM, 
and if provides a mechanism fo ensure fhaf fhe models do nof vi¬ 
olate flavor consfrainfs. As a consequence, spin -0 resonances musf 
have couplings fo fermions proportional fo fhe SM Higgs couplings. 
Flavor-safe models can still be consfrucfed beyond fhe MFV as¬ 
sumption, for example [ABG14], and deserve furfher sfudy. For a 
discussion of MFV in fhe confexf of fhe simplified models included 
in fhis reporf, see Ref. [Abd+15]. 

In fhe parameter scan for fhe models considered in fhis reporf, 
we make fhe assumption of a minimal decay widfh for fhe parficles 
mediating fhe interacfion befween SM and DM. This means fhaf 
only decays sfricfly necessary for the self-consistency of fhe model 
(e.g. fo DM and fo quarks) are accounted for in fhe definition of 
the mediator width. We forbid any furfher decays fo ofher invisible 
parficles of fhe Dark Secfor fhaf may increase fhe widfh or produce 
sfriking, visible signafures. Sfudies wifhin fhis reporf show fhaf, for 
cuf-and-counf analyses, fhe kinematic disfribufions of many mod¬ 
els, and fherefore fhe sensifivify of fhese searches, do nof depend 
significanfly on fhe mediator widfh, as long as fhe widfh remains 
smaller fhan fhe mass of fhe parficle and fhaf narrow mediafors are 
sufficienfly lighf. 

The parficle confenf of fhe models chosen as benchmarks is 
limifed fo one single kind of DM whose self-inferacfions are nof 
relevanf for LHG phenomenology, and fo one t5rpe of SM/DM in¬ 
feracfion af a fime. These assumptions only add a limited number 
of new parficles and new interactions fo fhe SM. These simpli- 
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fied models, independently explored by different experimental 
analyses, can be used as starting points to build more complete 
theories. Even though this factorized picture does not always lead 
to full fheories and leaves ouf defails fhaf are necessary for fhe 
self-consisfency of single models (e.g. fhe mass generafion for me¬ 
diator parficles), if is a sfarfing poinf fo prepare a sef of disfincf 
buf complemenfary collider searches for Dark Maffer, as if leads fo 
benchmarks fhaf are easily comparable across channels. 


1.3 Choices of benchmarks considered in this report and param¬ 
eter scans 


Confacf inferacfion operafors have been ouflined as basis sef of fhe- 
orefical building blocks represenfing possible f5rpes of inferacfions 
befween SM and DM parficles in [Goo+io]. The approach followed 
by LHC searches (see e.g. Refs. [CMSi5b; ATLi5d] for recenf jef+^j 
Run-1 searches wifh fhe 8 TeV dafasef) so far has been fo simu- 
lafe only a priorifized sef of fhe possible operafors wifh disfincf 
kinemafics for fhe inferprefafion of fhe consfrainfs obfained, and 
provide resulfs fhaf may be reinferprefed in terms of fhe ofher op¬ 
erafors. This reporf inf ends fo follow fhis sfrafegy, firsfly focusing 
on simplified models fhaf allow fhe explorafion of scenarios where 
fhe mediafing scale is nof as large. In fhe limif of large mediator 
mass, fhe simplified models map onto fhe EFT operafors. Secondly, 
fhis reporf considers specific EFT benchmarks whenever neifher 
a simplified model complefion nor ofher simplified models yield¬ 
ing similar kinemafic disfribufions are available and implemenfed 
in one of fhe evenf generators used by bofh collaborations. This 
is fhe case for dimension-5 or dimension-7 operafors wifh direcf 
DM-elecfroweak boson couplings 3 . Considering fhese models as 
separate experimenfal benchmarks will allow fo fargef new sig¬ 
nal regions and help validate fhe confacf inferacfion limif of new 
simplified models developed fo complete fhese specific operafors. 
Resulfs from fhese EFT benchmarks should include fhe condition 
fhaf fhe momenfum fransfer does nof probe fhe scale of fhe infer¬ 
acfion; whenever fhere is no model fhaf allows a direcf mapping 
befween fhese two quanfifies, various options should be tested fo 
ensure a given fracfion of evenfs wifhin fhe range of applicabilify of 
fhe EFT approach. Experimenfal searches should in any case deliver 
resulfs fhaf are independenf from fhe specific benchmark tested, 
such as fiducial cross-sections fhaf are excluded in a given signal 
region. 

When choosing fhe poinfs fo be scarmed in fhe parameter space 
of fhe models, fhis reporf does nof quanfifafively consider con¬ 
sfrainfs fhaf are exfernal fo fhe MET+X analyses. This is fhe case 
also for resulfs from LHC experimenfs searching for mediator de¬ 
cays. The main reason for nof doing so in fhis reporf is fhe diffi- 
culfy of incorporafing fhese consfrainfs in a rigorous quanfifafive 
way wifhin fhe fimescale of fhe Forum. However, even if fhe pa- 


3 An example of a dimension-5 op¬ 
erator for scalar DM is described in 
Appendix A. Dimension-7 operators of 
DM coupling to gauge bosons exist in 
the literature, but they require a larger 
particle spectrum with respect to the 
models studied in this report. 
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rameter scans and the searches are not optimized with those con¬ 
straints in mind, we intend to make all information available to 
the community to exploit the unique sensitivity of colliders fo all 
possible DM signafures. 


1.4 Structure of this report and dissemination of results 

The reporf provides a brief fheorefical summary of fhe models con¬ 
sidered, sfarfing from fhe sef of simplified models and confacf in- 
feracfions puf forward in previous discussions and in fhe liferafure 
cifed above. Ifs main body documenfs fhe sfudies done wifhin fhis 
Forum fo idenfify a kinemafically disfincf sef of model paramefers 
fo be simulafed and used as benchmarks for early Run-2 searches. 
The implemenfafion of fhese sfudies according fo fhe sfafe of fhe arf 
calculafions is defailed, including insfrucfions on how fo esfimafe 
fheorefical uncerfainfies in fhe generafors used for fhese sfudies. 

The presenfafion of resulfs for EFT benchmarks is also covered. 

Chapfer 2 of fhis reporf is dedicafed fo simplified models wifh 
radiation of a hard objecf eifher from fhe inifial sfafe or from fhe 
mediafor. These models produce primarily monojef signafures, 
buf should be considered for all Bj+X searches. Chapfer 3 con- 
fains sfudies on fhe benchmark models for final sfafes specifically 
confaining an elecfroweak boson (W/Z/ j/H). In fhis case, bofh 
simplified models leading fo mono-boson signafures and confacf 
inferacfion operators are considered. Defails of fhe sfafe of fhe arf 
calculafions and on fhe implemenfafion of fhe simplified models 
in Monfe Carlo generafors are provided in Chapfer 4. Chapfer 5 
is devofed fo fhe freafmenf of fhe presenfafion of resulfs for fhe 
benchmark models from confacf inferacfion operators. Chapfer 6 
prescribes how fo esfimafe fheorefical uncerfainfies on fhe simula- 
fion of fhese models. Chapfer 7 concludes fhe reporf. 

Furfher models fhaf could be sfudied beyond early searches 
and fheir implemenfafion are described in Appendix A. For fhese 
models, eifher fhe implemenfafion could nof be fully developed by 
fhe fime of fhis reporf, or some of fhe grounding assumptions were 
nof fully mef. Some of fhese models have been used in previous 
ATLAS and CMS analyses and discussed fhoroughly wifhin fhe 
Forum. They are fherefore worfh considering for furfher sfudies 
and for Run-2 searches, since fhey lead fo unique £t+X signafures 
fhaf are nof shared by any ofher of fhe models included in fhis 
reporf. Appendix B confains fhe necessary elemenfs fhaf should be 
included in fhe resulfs of experimenfal searches fo allow for furfher 
reinferprefafion. 

If is crucial for fhe success of fhe work of fhis Forum fhaf fhese 
sfudies can be employed as cross-check and reference fo fhe fhe¬ 
orefical and experimenfal communify interested in early Run-2 
searches. For fhis reason, model files, parameter cards, and cross- 
sections for fhe models considered in fhese sfudies are publicly 
available. The SVN reposifory of fhe Forum [Fork] confains fhe 
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models and parameter files necessary to reproduce the studies 
within this report. Details and cross-sections for these models, as a 
function of their parameters, will be published on HEPData [Hep]. 
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2 


Simplified models for all Zt +X analyses 


In this Chapter we review models that yield X+£ j signatures, 
where X is a QCD parton or j,W,Z or h. 

The primary simplified models for Dirac fermion DM sfudied 
and recommended by fhis Forum for early LHC Run-2 searches are 
defailed in fhis Chapfer, comprising spin -0 and spin -1 mediators. 

Section 2.1 covers fhe s-charmel exchange of a vecfor mediator ' Colored vector mediators can be 

while we consider bofh s-charmel and f-channel exchange for scalar exchanged m the t-channel, but there 

are no examples in literature so far. 

mediafors in Section 2.2 and 2.3 respectively Spin -2 mediafors are 
briefly mentioned in Secfion 2.4. While fhese models are general 
and cover a broad sef of signafures, fhe discussion and sfudies 
are focused on fhe monojef final sfafe. Defails on final sfafes wifh 
elecfroweak (EW) boson radiation and wifh heavy flavor quarks 
from diagrams arising wifhin fhese models are also discussed in 
fhis Chapfer. 

A summary of fhe sfafe of fhe arf calculations and implemenfa- 
fions for fhese models is provided in Table 6.1. Secfion 4 defails fhe 
implemenfafion of fhese models fhaf have been used for fhe sfud¬ 
ies in fhis Chapfer and fhaf will be employed for fhe simulafion of 
early Run-2 benchmark models for LHC DM searches. 


2.1 Vector and axial vector mediator, s-channel exchange 

A simple exfension of fhe Sfandard Model (SM) is an additional 
Lr(l) gauge symmefry, where a Dark Matter candidafe particle 
has charges only under fhis new group. Assuming fhaf some SM 
particles are also charged under fhis group, a new gauge boson can 
mediafe interactions between fhe SM and DM. 

We consider fhe case of a DM particle x of mass fhaf is a 
Dirac fermion and where fhe producfion proceeds via fhe exchange 
of a spin-1 mediafor of mass io the s-charmel, illusfrafed in 

Fig. 2.1. 

We consider two models with vector and axial-vector couplings 
between the spin -1 mediator Z' and SM and DM fields, with the 
corresponding interaction Lagrangians: 
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q s xi^x) 



Figure 2.1: Representative Feynman 
diagram showing the pair production 
of Dark Matter particles in association 
with a parton from the initial state via 
a vector or axial-vector mediator. The 
cross section and kinematics depend 
upon the mediator and Dark Matter 
masses, and the mediator couplings to 
Dark Matter and quarks respectively: 
(ATmed/ SX’ S'q)- 


^vector =^q ^ q + X (2-l) 

q—ii,d,s,c,b,t 

Axial-vector = ^q E H + gx^iXX^X- (2-2) 

q—ii,d,s,c,b,t 

The coupling is assumed to be universal to all quarks. It is also 
possible to consider other models in which mixed vector and axial- 
vector couplings are considered, for instance the couplings to the 
quarks are axial-vector whereas those to DM are vector. As men¬ 
tioned in the Introduction, when no additional visible or invisible 
decays contribute to the width of fhe mediafor, fhe minimal widfh 
is fixed by fhe choices of couplings and g^. The effecf of larger 
widfhs is discussed in Secfion 2.5.2. For fhe vector and axial-vecfor 
models, fhe minimal widfh is: 




127T ' m2 


^DM0(Mmed(2-3) 


^ ^ ^S^q-^med 


127T 


med / 

2 mi \ 

A0(M^ed - 2m^), 


M- 


med,, 




(24) 


V 4m^ 

1 — 2'' 

^med 

is fhe velocify of fhe fermion / wifh mass nif \t\ fhe mediafor 
resf frame. Note fhe color facfor 3 in fhe quark ferms. Figure 2.2 
shows fhe minimal widfh as a funcfion of mediafor mass for bofh 
vecfor and axial-vecfor mediators assuming fhe coupling choice 
Sq_ ~ Sx ~ fhis choice of fhe couplings, fhe dominanf con- 

fribufion to fhe minimal widfh comes from fhe quarks, due fo fhe 
combined quark number and color facfor enhancemenf. We specif¬ 
ically assume fhaf fhe vecfor mediafor does nof couple fo leptons. 
If such a coupling were presenf, if would have a minor effecf in in¬ 
creasing fhe mediafor widfh, buf if would also bring in consfrainfs 
from measuremenfs of fhe Drell-Yan process fhaf would urmeces- 
sarily resfricf fhe model space. 










rSED [GeV] r™"D [GeV] 
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Figure 2.2: Minimal width as a func¬ 
tion of mediator mass for vector and 
axial-vector mediator assuming cou¬ 
plings of 1. The total width is shown 
as solid lines for Dark Matter masses 
of 10 GeV, 30 GeV, 100 GeV and 
300 GeV in black, red, brown and 
green, respectively. The individual 
contributions from Dark Matter are 
indicated by dotted lines with the 
same colors. The contribution from all 
quarks but top is shown as magenta 
dotted line and the contribution from 
top quarks only is illustrated by the 
dotted blue line. The dotted black line 
shows the extreme case T^in = Mmed- 


•^MED [GgV] 
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Therefore, the minimal set of parameters under consideration for 
these two models is 


{ ^med. } ■ (2.5) 

together with the spin structure of their couplings. 

A thorough discussion of these models and their parameters can 
also be found in [Buc+15]. 

These simplified models are known and available in event gen¬ 
erators at NLO + PS accuracy, as detailed in Section 4.1.1. Results 
in this Section have been obtained using the model implementa¬ 
tion within the fowheg generator (V3359) [HKR13], interfaced to 
PYTHiA 8 [SMS08] for the parton shower. 

In addition, for the vector models considered, initial and final 
state radiation of a Z' can occur which can appear as a narrow jet if 
it decays hadronically and may not be distinguishable from a QCD 
jet, thus accounting for some fraction of the monojet signal. The 
ISR and FSR of Z' becomes more important at large values of the 
couplings [BBL15]. 

2.1.1 Parameter scan 

In order to determine an optimal choice of the parameter grid for 
the simulation of early Run-2 benchmark models, dependencies 
of the kinematic quantities and cross sections on the model pa¬ 
rameters have been studied. Only points that are kinematically 
distinct will be fully simulated, while instructions on how to rescale 
the results according to models with different cross sections are 
presented in Section 2.5. The following paragraphs list the main 
observations from the scans over the parameters that support the 
final proposal for the benchmark signal grid. 

2.1.1.1 Scan over the couplings 

To study the dependence of kinematic distributions on the coupling 
strength, samples were generated where a pair of = 10 GeV 
Dark Matter particles is produced on-shell from the mediator of 
Mmed = 1 TeV. Figure 2.3 compares the shapes of the Bj distri¬ 
bution for the different choices of the coupling strength. This is a 
generator-level prediction with no kinematic selections or detec¬ 
tor simulation. Coupling values in the scan range 0.1-1.45, fixing 
yq = g^, correspond to a rough estimate of the lower sensitivity 
of mono-jet analyses and a maximum coupling value such that 
Tmin < Mjned- observe that the shapes of the £7 or jet pj dis¬ 
tributions do not depend on the couplings (and consequently the 
width) in the ranges considered. A large width of the mediator im¬ 
plies a broad integral over the contributing parton distributions, 
which might not be well approximated by the midpoint of this in¬ 
tegral. This study shows that the effect, in the pj distribution of the 
observed gluon, is not important. 
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[GeV] 

Based on similar findings for differenf choices of Mj^ed ^^id nt^, 
we conclude fhaf fhe shapes of kinemafic disfribufions are nof 
alfered by coupling variations, neifher for fhe on-shell mediafor 
case where M^ed > off-shell case where M^ed < 

2m^. Only fhe production cross sections change. Differences in 
kinemafic disfribufions are expecfed only close fo fhe fransifion 
region befween on-shell and off-shell mediators. 

Special care needs fo be faken when coupling sfrengfhs are com¬ 
bined wifh exfremely heavy mediators. Figure 2.4 suggesfs a change 
in fhe shape of fhe E j disfribufion for a Mj^gd = 5 TeV mediafor 
once Fmin/ Mined is of fhe order of a percenf or lower. 
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Figure 2.3: Scan over couplings. The 
Bt distribution is compared for the 
vector mediator models using the 
parameters as indicated. Ratios of 
the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote 
the acceptance of the Bt > 300 GeV 
and Bt > 500 GeV cut, respec¬ 
tively. All figures in this Section have 
been obtained using the model im¬ 
plementation within the powheg 
generator (V3359) [HKR13], interfaced 
to PYTHiA 8 [SMS08] for the parton 
shower. 
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Such heavy mediators, although inaccessible with early LHC 
data, are interesting since they provide a good approximation for 
benchmark EFT models. The observed difference among the sim¬ 
plified models in the plot arises from the fact that the region of 


Figure 2.4: Comparison of the Bt 
distributions from the D5 EFT sample 
and the vector models with 5 TeV 
heavy mediator of various widths. 
Ratios of the normalized distributions 
with respect to the first one are shown. 
A300 and A500 in the table denote the 
acceptance of the Bt > 300 GeV and 
Bt > 500 GeV cut, respectively. 
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low invariant masses of the Dark Matter pair, is suppressed 
due to narrow Breit-Wigner peak that only probes a narrow win¬ 
dow of parton distribution functions. For wider mediators, the 
low mass region is significantly enhanced by parton distribution 
functions at low Bjorken x, as illustrated in Fig. 2.5(a). This ex¬ 
plains why the sample with the narrowest mediator in Fig. 2.4 is 
heavily suppressed in terms of production cross section and also 
gives different £7 shape. Furthermore, Fig. 2.4 compares the vector 
model with 5 TeV mediator to the D5 EFT sample and reveals that 
the simplified models with larger mediator widths (e.g. for cou¬ 
plings of 1 where Tmin/Mmed ^ 0 . 5 ) are the ones resembling the 
kinematics of contact interactions. This reflects the fact that in an 
EFT there is no enhancement due to on-shell mediators, leading to 
a closer resemblance to an off-shell regime where no peak in the 
distribution is present. In case of narrow width mediators, e.g. 
Tmin/Mmed ~ 0 . 05 , even larger mediator masses need to be chosen 
in order to significantly suppress the peak in the nixx distribution 
and reproduce the kinematic shapes of an EFT model. Figure 2.5(b) 
verifies that the choice of 10 TeV mediator mass is sufficient to 
achieve that. 

Since kinematic distributions are robust to changes in the specific 
values of coupling the choice of =0.25 and =1 is reasonable 
to reduce the parameter space to be scarmed. There are no com¬ 
plications associated with small couplings, but, also, the early part 
of Run 2 will not be sensitive to them. The range of couplings we 
recommend to generate limit the calculated width of the mediator 
to be near or below Mmed- 

For direct mediator searches, such as qq ^ Z' ^ qq, different 
couplings (yq 7^ g^) might also be considered. A scan in g^ vs ^q 
can then be performed for a fixed mediator mass. Such searches 
may restrict yq to a greater degree than g-^. 

2.1.1.2 Scan over 

For a fixed mediator mass M^ed arid couplings, the Dark Matter 
mass falls into three regimes: 

On-shell: When Mj^gd ^ rnost mediators are on-shell. The 
hardness of the ISR is set by Mj^ed/ and the kinematic distribu¬ 
tions do not strongly depend on m^. This is illustrated in Fig. 2.6 
for an example of M^ed =1 TeV 10 GeV < < 300 GeV. The 

cross section decreases as the approaches Mined/ 2 . A coarse 
birming along is sufficient. 

Threshold: When M^ed ~ the production is resonantly 

enhanced, and both the cross section and kinematic distributions 
change more rapidly as a function of the two masses, and finer 
birming is needed in order to capture the changes. 


’ This applies as long as heavy narrow 
mediators are generated without any 
truncation of low-mass tails at the 
generator-level. 


Off-shell: When M^ed ^ the Dark Matter pair is produced by 
an off-shell mediator. The mediator propagator gives an explicit 
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Figure 2.5: Invariant mass of the Dark 
Matter pair in the vector mediator 
samples with = 10 GeV, Mmej = 

5 TeV and different coupling strengths 
(a). A similar comparison is shown for 
the samples with different mediator 
masses considering Tmin/Mmed = 0.05 
and 0.1 (b). An EFT sample is also 
displayed in the latter case. The 
distributions are normalised to unit 


area. 
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suppression of (Mmed/Q)^ that suppresses hard ISR. The = 

1 TeV case, shown in Fig. 2.6, and Figure 2.7 demonstrates that 
the £7 spectrum hardens with increasing m^, accompanied by 
the gradual decrease of the cross section. Due to the significant 
cross section suppression, it is not necessary to fully populate the 
parameter space. Imminent LFIC searches are not expected to be 
sensitive to these signals. 
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Figure 2.6: Scan over Dark Matter 
mass. The distribution is compared 
for the vector mediator models using 
the parameters as indicated. Ratios 
of the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote the 
acceptance of the E7 > 300 GeV and 
Bt > 500 GeV cut, respectively. 


2.1.1.J Scan over the mediator mass 

Changing the mediator mass for fixed Dark Matter mass and cou¬ 
plings leads to significant differences in cross section and shapes of 
the kinematic variables for the on-shell regime, as shown in Fig. 2.8. 
As expected, higher mediator masses lead to harder £7 spectra. On 
the other hand, the £7 shapes are similar for off-shell mediators. 
This is illustrated in Fig. 2.9. Therefore, a coarse birming in M^ed is 
sufficient in the off-shell regime. 

2.1.1.4 Spin structure of the couplings 

This section compares the kinematic properties of vector, axial- 
vector and mixed vector/axial-vector models. The samples with 
pure vector and pure axial-vector couplings are compared for 
Mmed = 100 GeV and different Dark Matter masses in Fig. 2.10. 

No differences in the shape of the £7 distributions are observed 
between the samples with coincident masses. In the case of the on- 
shell mediators, where 2 m^ <C M^ed/ the cross sections of the pure 
vector and pure axial-vector models are similar. With increasing 
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Figure 2.7: Scan over Dark Matter 
mass. The Er distribution is compared 
for the vector mediator models using 
the parameters as indicated. Ratios 
of the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote the 
acceptance of the > 300 GeV and 
Bt > 500 GeV cut, respectively. 
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Figure 2.8: Scan over mediator mass. 
The Bt distribution is compared for 
the vector mediator models using 
the parameters as indicated. Ratios 
of the normalized distributions with 
respect to the first one are shown. 
A300 and A500 in the table denote the 
acceptance of the Bt > 300 GeV and 
Bt > 500 GeV cut, respectively. 
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Dark Matter mass towards the 2m^ = Mj^ed transition and fur¬ 
ther into the off-shell regime, fhe relative difference between fhe 
cross sections of fhe fwo samples is increasing, wifh fhe vector ones 
having larger cross sections. 
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Figure 2.9: Scan over mediator mass. 
The Hr distribution is compared for 
the vector mediator models using 
the parameters as indicated. Ratios 
of the normalized distributions with 
respect to the first one are shown. 
A300 and A500 in the table denote the 
acceptance of the Hr > 300 GeV and 
Hr > 500 GeV cut, respectively. 
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Figure 2.11 shows fhe samples generafed wifh pure and mixed 
couplings for = 100 GeV and Mmed = 1 TeV, i.e. where fhe 
mediafor is on-shell. The mediator widfh befween fhe pure vecfor 
and pure axial-vecfor couplings differ only by 2 % in fhis case, and 
< 10% agreemenf befween fhe cross sections is found. The media¬ 
for widfhs for fhe samples wifh fhe same type coupling to quarks 
agree at better than 1 % since the width is dominated by the quark 
contribution, as expected from Eq. 2 . 3 . No significant differences be¬ 
tween the samples with same type Dark Matter coupling are seen, 
given the statistical precision of the generated samples. This is ex¬ 
pected since the mediator is on-shell, and the details of the invisible 


Figure 2.10: Comparison of the pure 
vector and pure axial-vector couplings. 
The Hr distribution is shown for the 
samples generated with Mmed = 

100 GeV and different Dark Matter 
masses. Ratios of the normalized 
distributions are shown for between 
the samples with coincident masses. 
A300 and A500 in the table denote the 
acceptance of the Hr > 300 GeV and 
Hr > 500 GeV cut, respectively. 
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decay are unimportant in cut-and-count searches. 

For the off-shell case, shown in Fig. 2.12 for = 100 GeV and 
^med = 100 GeV, fhere is approximafely a facfor 2 difference 
befween fhe cross-secfions of fhe samples wifh pure couplings is 
observed. As in fhe previous case, fhe samples wifh fhe same type 
coupling to Dark Matter are similar both in terms of cross secfions 
and Bt shape. Since fhe confribufion fo fhe mediafor widfh from 
Dark Maffer is closed in fhis case, only fhe quark couplings define 
fhe widfh. Only couplings fo lighf quarks are opened in fhe case 
of Mjngd = 100 GeV for which fhe differences befween fhe partial 
widfhs of vecfor and axial-vecfor couplings are marginal. This 
explains fhe similar minimal widfhs for all four samples sfafed in 
Fig. 2.12. 

In general, fhe coupling fo quarks is nof expecfed fo play an 
imporfant role in fhe kinematics as if is only needed fo produce 
the mediator which is confirmed by fhe observations above. Based 
on this argument and on the observations above, we recommend 
to consider only the models with pure vector couplings or pure 
axial-vector couplings for simulafion. 
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Figure 2.11: Comparison of the pure 
vector, V-V, and pure axial-vector. A-A, 
couplings with mixed couplings, A-V 
and V-A where the first (second) letter 
indicates the Standard Model (dark 
sector) vertex. The tj distribution is 
shown for the samples generated with 
ni^ = too GeV and Mmej = 1 TeV. 
Ratios of the normalized distributions 
are shown for A-V over V-V and for 
V-A over A-A. A300 and A500 in the 
table denote the acceptance of the 
Bt > 300 GeV and Bt > 500 GeV cut, 
respectively. 


2.1.1.^ Proposed parameter grid 

The final sfep in proposing a paramefer grid is fo evaluafe fhe sen- 
sifivify of Run-2 LHG dafa wifh respecf fo rafe and/or kinemafics. 
The paramefer scan focuses on fwo imporfanf regions, fhe lighf 
mediafor region and fhe heavy mediafor limif fo reproduce fhe 
EFT limif, and fakes info accounf fhe projecfed sensitivities for fhe 
mono-jef analysis. 

Gonsidering simplified models also allows fo discuss consfrainfs 
from differenf search channels. In fhe case of fhe s-channel ex¬ 
change, fhe resulfs from fhe mono-jef final sfafes, where fhe medi¬ 
afor decays fo a DM pair, one can also fake info account dijef con¬ 
sfrainfs on fhe processes where fhe mediafor decays back to Stan- 
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dard Model particles. The importance of the dijet results depend on 
the magnitude of fhe coupling ^q. We recommend to keep the two 
channels rather independent by choosing yq = 0.25 and = 1 , 
based on the findings given in Ref. [Cha+15]. Furfhermore, if is also 
imporfanf fo menfion fhis choice leads fo kmin/^ 0 . 06 . ISIofe 
fhaf fhe usual choice of gq = g;;( = 1 used in liferafure leads fo 
Tmin/Mj^ed ^ 0 - 5 / quesfioning fhe applicabilify of fhe narrow widfh 
approximafion. 

The expecfed upper limif af 95% confidence level on fhe prod- 
ucf of cross secfion, accepfance and efficiency, cr x A x e, in fhe 
final Run-1 ATLAS mono-jef analysis [ATLi5d] is 51 fb and 7.2 fb 
for Bj > 300 GeV and Bt > 500 GeV, respecfively. Projecfed 
sensifivifies for a 14 TeV mono-jef analysis are available from 
ATLAS [ATLi4d]. These ATLAS sfudies esfimafe a factor of fwo 
increase in sensifivify with the 2015 data. The generator level cross 
section times efficiency fimes accepfance at Bj > 500 GeV for fhe 
model wifh couplings yq = 0.25 and g^ = 1 , a light Dark Matter 
particle of =10 GeV and a M^gd =1 TeV vector mediafor is at 
the order of 100 fb, i.e. fhe early Run-2 mono-jef analysis is going fo 
be sensifive fo heavier mediafors fhan fhis. The value of <7 x e x A 
af E T > 500 GeV for a 5 TeV vecfor mediafor is af fhe order of 
0.1 fb, fherefore fhis model lies beyond fhe reach of fhe LHG in fhe 
early Run-2. However, models wifh high enough mediafors are sfill 
useful fo reproduce fhe EFT resulf. 

Following fhese argumenfs, M^ed grid poinfs are chosen, roughly 
equidisfanf in a logarifhmic scale: 10 GeV, 20 GeV, 50 GeV, 100 GeV, 
200 GeV, 300 GeV, 500 GeV, 1000 GeV and 2000 GeV. In fhe 
fhreshold regime fhe grid poinfs are faken af ap¬ 

proximately Mingd/ 2 , namely: 10 GeV, 50 GeV, 150 GeV, 500 GeV 
and 1000 GeV. Poinfs on the on-shell diagonal are always chosen to 
be 5 GeV away from fhe fhreshold, fo avoid numerical insfabilifies 
in fhe event generation. The detailed studies of fhe impacf of fhe 
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Figure 2.12: Comparison of the pure 
vector, V-V, and pure axial-vector. A-A, 
couplings with mixed couplings, A-V 
and V-A where the first (second) letter 
indicates the Standard Model (Dark 
Sector) vertex. The Er distribution is 
shown for the samples generated with 
= too GeV and = 100 GeV. 
Ratios of the normalized distributions 
are shown for A-V over V-V and for 
V-A over A-A. A300 and A500 in the 
table denote the acceptance of the 
Bt > 300 GeV and Bt > 500 GeV cut, 
respectively. The suppression by for 
~ Aimed can be seen for the curves 
representing axial DM coupling. 
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parameter changes on the cross section and kinematic distributions 
presented earlier in this section support removing some of the grid 
points and relying on interpolation. The optimized grids proposed 
for fhe vecfor and axial-vecfor mediafors are given in Table. 2.1. One 
poinf af very high mediafor mass (10 TeV) is added for each of fhe 
DM masses scanned, fo aid fhe reinferprefafion of resulfs in ferms 
of confacf inferacfion operafors (EFTs), as discussed in Secfion 5.2. 


m^/ GeV 
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10 
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1995 

10000 


Table 2.1: Simplified model bench¬ 
marks for s-channel simplified models 
(spin-l mediators decaying to Dirac 
DM fermions in the V and A case, 
taking the minimum width for gq = 
0.25 and = 1) 


Tables 2.2 and 2.3 give the Tmin/Mjned ratio for the parameter 
grid proposed for vector and axial-vector s-charmel models, respec¬ 
tively. The numbers range from ^ 0.02 in the off-shell regime at 
2 m^ > Mjngd to ^ 0.06 in the on-shell regime for heavy mediators 
where all coupling charmels contribute. 
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0.057 
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Table 2.2: Minimal width of the vector 
mediator exchanged in s-channel di¬ 
vided by its mass, assuming gq = 0.25 
and gx = The numbers tabulated 
under 2m^ = Mmed correspond to the 
width calculated for — 5 GeV. 


nix / GeV 
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0.045 
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0.051 

0.051 
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0.056 
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0.020 







0.029 
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Table 2.3: Minimal width of the 
axial-vector mediator exchanged in 
s-channel divided by its mass, as¬ 
suming g^q = 0.25 and gx = 1. The 
numbers tabulated under 2m^ = Mmed 
correspond to the width calculated for 
Mmed - 5 GeV. 


2.1.2 Additional considerations for V+Bt signatures 

All models detailed in this Section are applicable to signatures 
where a photon, a W boson, a Z boson or a Higgs boson is radiated 
from the initial state partons instead of a gluon. The experimental 
signature is identified as V+Bt and it has been sought by ATLAS 
and CMS in Refs. [CMSi4b; ATL15C; CMSi5e; ATL14C; ATLi4a; 
ATLi4b]. This signature is also produced by the models described 
in Section 3. 

Monojet searches are generally more sensitive with respect to 
final states including EW bosons, due to the much larger rates of 
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signal events featuring quark or gluon radiation with respect to ra¬ 
diation of bosons [ZBW13], in combination wifh fhe low branching 
ratios if lepfons from boson decays are required in fhe final sfafe. 
The rafes for fhe Higgs boson radiation is foo low for fhese models 
fo be considered a viable benchmark [Car+14]. However, fhe pres¬ 
ence of phofons, lepfons from W and Z decays, and W or Z bosons 
decaying hadronically allow backgrounds fo be rejecfed more ef¬ 
fectively, making Z/ j/W+Bj searches sfill worfh comparing wifh 
searches in fhe jef+^j final sfafe (see e.g. Ref. [Ger+08]). 

In fhe case of a spin -1 mediafor, an example Feynman diagram 
for fhese processes can be consfrucfed by faking Fig. 2.1 and replac¬ 
ing fhe gluon wifh 7, W or Z. 

When fhe inifial sfafe radiafion is a W boson, Run-i searches 
have considered fhree benchmark cases, varying fhe relafive cou¬ 
pling of fhe W fo M and d quarks. The simplified model wifh a 
vecfor mediafor mediafor exchanged in fhe s-charmel includes only 
fhe simplesf of fhese cases, in which fhe W coupling fo u and d 
quarks is identical, as required naively by SU( 2 ) gauge invariance. 
Wifh some more complex model building, ofher cases are possible. 
The case in which fhe u and d couplings have opposife sign is par- 
ficularly inferesfing, since fhis enhances fhe W + Bt signal over fhe 
jef-h^T signal [Bel+i5b; BT13; Ham+14]. An example of a model of 
fhis type is discussed in Appendix A. 2. 

Simulations for fhe models in fhis Section have been done af fhe 
LO+PS level using MadGraph5_aMC@NLO 2.2.3 inferfaced fo 
PYTHiA 8, and fherefore no special runtime configuration is needed 
for pyfhia 8. Even fhough merging samples wifh differenf par ton 
mulfiplicifies is possible, fhis has nof been deemed necessary as 
fhe visible signal comes from fhe production of a heavy SM bo¬ 
son whose fransverse momenfum disfribufion is sufficienfly well 
described af LO+PS level. 

In fhese V+Bt models, as in fhe case of fhe jef+^p models, pj of 
fhe boson or fhe B t does nof depend sfrongly on fhe widfh of fhe 
mediafor. An example of fhe parficle-level analysis accepfance us¬ 
ing fhe generator-level cufs from Ref. [ATL15C] for fhe phofon+^j 
analysis, buf raising fhe phofon pj cuf fo 150 GeV, is shown in Fig¬ 
ure 2.4, comparing a widfh fhaf is sef fo F = / 3 fo fhe minimal 

widfh (fhe ratio befween fhe fwo widfhs ranges from 1.05 fo 1.5 
wifh increasing mediafor masses). 


Accepfance ratio for F = Fjnin vs F 

-^med/^ 


nij^/GeV 

Mmed/GeV 

10 

50 200 

400 

50 

0.96 

0.99 

0.95 

100 

0.97 



300 

1.00 

1.02 


600 


0.96 


1000 

1.01 

1.02 1.03 


3000 

1.02 

1.03 

1.01 


Table 2.4: Analysis acceptance ratios 
for the photon+ET" analysis when 
varying the mediator width, in the 
case of a vector mediator exchanged 
in the s-channel. The figures shown 
in this Section have been obtained 
using a LO UFO model in Mad- 
Graph5_aMC@NLO 2.2.3 interfaced 
to PYTHIA 8 for the parton shower. 
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Examples of relevant kinematic distributions for selected bench¬ 
mark points are shown in Fig. 2.13. 




p’ (pre-selection) [GeV] p([. (pre-selection) [GeV] 


(a) Leading photon transverse momentum distribution for the (b) Leading photon transverse momentum distribution for the 

photons-Hr final state, for different mediator mass choices, for photon-r^r final state, for different DM mass choices, with Mmed 

=10 GeV. =1 TeV. 



Emiss jQgyj 


(c) Missing transverse momentum distribution for the leptonic 
Z+'Bt final state, for different mediator mass choices, for 
=15 GeV 



(d) Missing transverse momentum distribution for the hadronic 
W+Bt final state. 

Figure 2.13: Kinematic distributions 
relevant for searches with W, Z and 
photons in the final state, for the sim¬ 
plified model with a vector mediator 
exchanged in the s-channel. 


2.2 Scalar and pseudoscalar mediator, s-channel exchange 


In this section, we consider a parallel situation to the vector and 
axial-vector mediators in the previous sections: a real scalar or a 
pseudoscalar where the associated scalar is decoupled at higher 
energies^. This section is largely based on Refs. [BFG15; Har-H5; 
HR15] which contain a thorough discussion of these models. 

Assuming MFV, spin -0 resonances behave in a similar fashion as 
the SM Higgs boson. If the mediators are pure singlets of the SM, 
their interactions with quarks are not SU{ 2 )i invariant. To restore 
this invariance, one could include the mixing of such mediators 
with the Higgs sector. This leads to extra interactions and a more 


3 Ttiis assumption does not hold in a 
UV-complete model where the two 
components of the complex scalar 
mediator would be approximately 
degenerate. The complex scalar case 
could be studied separately in the case 
of heavy flavor final states given the 
sufficiently different kinematics. 
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Figure 2.14: One-loop diagrams of 
processes exchanging a scalar (S) or 
pseudoscalar (P) mediator, leading to a 
mono-jet signature. 




complex phenomenology with respect to what considered in this 
Section (for a more complete discussion, see Refs. [BFG15; HR15]). 
In fhe inferesf of simplicify, we do nof sfudy models including 
fhose inferacfions in fhis reporf as early Run-2 benchmark models, 
buf we give an example of a model of fhis kind in Appendix A.4. 

Relafive fo fhe vector and axial-vecfor models discussed above, 
fhe scalar models are disfinguished by fhe special consequences 
of fhe MFV assumption: fhe very narrow widfh of fhe mediafor 
and ifs exfreme sensifivify fo which decays are kinematically avail¬ 
able, and fhe loop-induced coupling fo gluons. The interaction 
Lagrangians are 


= 


gx 4 ’XX + ^ E {suy“uiUi + gay^didi + geyfliii'j , 


U = igx^XlbX + E {guyiUijsUi + gdy‘!dij5di+ 


giy 




(2.6) 


( 2 - 7 ) 


where (p and a are respecfively fhe scalar and pseudoscalar media- 
fors, and fhe Yukawa couplings i/| are normalized fo fhe Fliggs vev 
as /v. 

The couplings fo fermions are proportional fo fhe SM Higgs 
couplings, yef one is sfill allowed fo adjusf an overall sfrengfh of fhe 
coupling fo charged leptons and fhe relafive couplings of u- and d- 
fype quarks. As in fhe preceding secfions, for fhe sake of simplicify 
and sfraighfforward comparison, we reduce fhe couplings fo fhe 
SM fermions fo a single universal parameter gq = gu = gd = gl- 
Unlike fhe vector and axial-vecfor models, fhe scalar mediators are 
allowed fo couple fo leptons .4 

The relafive discovery and exclusion power of each search can 
be compared in fhis framework. However, we again emphasize fhe 
imporfance of searching fhe full sef of allowed channels in case vio- 
lafions of fhese simplifying assumpfions lead fo significanf modifi¬ 
cations of fhe decay rafes fhaf unexpecfedly favor differenf channels 
fhan fhe mix obfained under our assumpfions. The coupling g-^ 
paramefrizes fhe entire dependence on fhe sfrucfure between fhe 
mediafor and fhe dark secfor. 


+ This contribution plays no role 
for most of the parameter space 
considered. The choice to allow 
lepton couplings follows Refs. [BFG15; 
Har+15]. 
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Given these simplifications, the minimal set of paramefers under 
consideration is 


{m^, m^/a = Mmed. gx' ^q} ' (2-8) 

Fig. 2.14 shows fhe one-loop diagrams producing a jef+X signafure. 
The full calculation of fhe fop loop is available af LO for DM pair 
production in association wifh one parfon. 

The minimal mediafor widfh (neglecting fhe small confribufions 
from quarks ofher fhan fop in fhe loop) is given by 


T 


(p,a — 


f 


y)g\^<pA 

\ 6 n 



2 2 1 ^ 




2 


Slm<PA 

8 n 



where x = 3 for scalars and x = 1 for pseudoscalars. The loop 
infegrals, wifh / as complex functions, are 


/«(t) 


1 -|- (1 — t) arcfan^ 


Vt- 1 J\ 


T arcfan 


Vt- 1 


where t = 

The minimal widfhs for scalar and pseudo-scalar mediafors 
wifh gq = gx — ^ ^re shown in Fig. 2.20, illusfrafing fhe effecf of 
choosing fhe SM Higgs-like Yukawa couplings for fhe SM fermions. 
For fhe mediafor mass above fwice fhe fop quark mass nit, the 
minimal width receives the dominant contribution from the top 
quark. For lighter mediator masses. Dark Matter dominates as the 
couplings to lighter quarks are Yukawa suppressed. 

As shown in the diagram of Fig. 2.14, the lowest order process of 
these models already involves a one-loop amplitude in QCD, and 
only LO predictions are currently available. The generator used 
for the studies for the jet+^j signature is fowheg [HKR13; HR15; 
Ali+10; Naso4; FNO07], with pythia 8 [SMS08] for the parfon 
shower; within this implementation, the scalar and pseudoscalar 
mediator benchmark models are known at LO+PS accuracy. 


2.2.1 Parameter scan 

Similarly as in the case of the vector and axial-vector couplings of 
spin-1 mediators, scans in the parameter space are performed also 
for the scalar and pseudo-scalar couplings of the spin-0 mediators 
in order to decide on the optimized parameter grid for the pre¬ 
sentation of Run-2 results. Figures 2.15- 2.19 show the scans over 
the couplings. Dark Matter mass and mediator mass and the same 
conclusions apply as in Section 2.1. 


( 2 - 9 ) 


(2.10) 

(2.11) 
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A scan over the mediator mass is shown in Fig. 2.19 where M^ed 
= 300 GeV and 500 GeV are chosen to be below and above 2 mf. The 
off-shell case is assumed by taking an extreme limit {m^ = 1 TeV) 
in order to study solely the effects of fhe couplings fo quarks. No 
differences in fhe kinemafic disfribufions are observed and also fhe 
cross secfions remain similar in fhis case. No significanf changes 
appear for mediator masses around fhe 2 mt fhreshold. 
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If can be seen in Fig. 2.21 fhaf fhe kinematics for fhe scalar and 
pseudoscalar models coincides when considering fhe diagrams in 
Fig. 2.14. For fhis reason, we recommend fo fully simulate only 
one of fhe fwo models. No preference is given befween fhe two 
models as they have the same kinematics, although it is worth not¬ 
ing that the pseudo-scalar model has been used for a Dark Maffer 
inferprefafion of fhe DAMA signal and of fhe galactic cenfer ex¬ 
cess [ADNP15]. Like in the case of fhe vecfor and axial-vecfor mod¬ 
els described in Section 2.1.1.4, fhe differences befween fhe cross 
secfions for fhe scalar and pseudo-scalar samples wifh fhe same 
and Mmed ^re increasing with the Dark Matter mass for fixed 
mediafor mass, wifh fhe pseudo-scalar model yielding larger cross 
secfions. There is an increasing difference befween fhe minimal 
widfhs close fo fhe 2 m^ = Mj^ed threshold. 


Figure 2.15: Scan over couplings. The 
Bt distribution is compared for the 
scalar mediator models using the 
parameters as indicated. Ratios of 
the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote 
the acceptance of the Bt > 300 GeV 
and Bt > 500 GeV cut, respectively. 
Studies in all figures for the jet+^r 
signature is powheg, with pythia 8 for 
the parton shower; 


2.2.1.1 Proposed parameter grid 

The optimized paramefer grid in fhe Mjned~^x plsTie for scalar and 
pseudo-scalar mediators is mofivafed by similar argumenfs as in 
fhe previous secfion. Therefore, a similar pattern is followed here, 
wifh fhe excepfion of faking gq = gx = 1 - The choice of = 0.25 
for fhe vecfor and axial-vecfor models is mofivafed by suppress¬ 
ing consfrainfs from di-jefs, which is nof a concern in fhe scalar 
and pseudo-scalar mediafor case. Here a di-jef signal emerges only 
af fhe 2-loop level fhrough diagrams where fhe mediafor is pro¬ 
duced via gluon-gluon fusion and decays back info fwo gluons 
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Figure 2.16: Scan over Dark Matter 
mass. The Er distribution is compared 
for the scalar mediator models using 
the parameters as indicated. Ratios 
of the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote the 
acceptance of the > 300 GeV and 
Bt > 500 GeV cut, respectively. 
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Figure 2.17: Scan over Dark Matter 
mass. The Bt distribution is compared 
for the scalar mediator models using 
the parameters as indicated. Ratios 
of the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote the 
acceptance of the Bt > 300 GeV and 
Bt > 500 GeV cut, respectively. 


Figure 2.18: Scan over mediator mass. 
The Bt distribution is compared for 
the scalar mediator models using the 
parameters as indicated. Ratios of 
the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote the 
acceptance of the Bt > 300 GeV and 
Bt > 500 GeV cut, respectively. 
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Figure 2.19: Scan over mediator mass. 
The Bt distribution is compared for 
the scalar mediator models using the 
parameters as indicated. Ratios of 
the normalized distributions with 
respect to the first one are shown. 

A300 and A500 in the table denote the 
acceptance of the Bt > 300 GeV and 
Bt > 500 GeV cut, respectively. 
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•^MED [GsV] 



Figure 2.20: Minimal width as a func¬ 
tion of mediator mass for scalar and 
pseudo-scalar mediator assuming cou¬ 
plings of 1. The total width is shown 
as solid lines for Dark Matter masses 
of =10 GeV, 30 GeV, 100 GeV and 
300 GeV in black, red, brown and 
green, respectively. The individual 
contributions from Dark Matter are 
indicated by dotted lines with the 
same colors. The contribution from all 
quarks but top is shown as magenta 
dotted line and the contribution from 
top quarks only is illustrated by the 
dotted blue line. The dotted beige line 
shows the contribution from the cou¬ 
pling to gluons. The dotted black line 
shows the extreme case T^m = Mmed- 
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DARK MATTER BENCHMARK MODELS FOR EARLY LHC RUN-2 SEARCHES: 


REPORT OF THE ATLAS/CMS DARK MATTER FORUM JJ 


[GeV] 

-•- S 10 

p 10 

-■- S 30 

p 30 

-- S 100 

----A---- p 100 


400 600 800 1000 1200 

[GeV] 


> 

CD 

o 


10 r 


c 
CD 
> 

LH 

LU 

T 3 

^ 10-2 


Is 




iS: 


^med = 300 GeV 

(s = 13TeV 
j Ldt = 1 fb‘^ 


:=& = 




10^^ r 


3 _ 


=•= 


—A- 



through a top loop. The strong loop suppression renders such sig¬ 
nals unobservable at the LHC. Further constraints on the scalar 
and pseudo-scalar mediators may emerge from searches in ti final 
sfafes. Sfudies of fhe elecfroweak effecfs fo tt producfion suggesf 
fhaf one can only expecf percenf level confribufions for yq r^j 0(1) 
[HHR14]. Therefore, keeping yq = = 1 is a reasonable choice 

in fhe case of fhe scalar and pseudo-scalar mediators. Confrary fo 
fhe vecfor and axial-vecfor models, nofe fhaf couplings of 1 lead 
to Tmin/Mined ~ 0-1/ ensuring fhe narrow widfh approximafion is 
applicable. Furfhermore, fhe sensifivify fo fhe highesf mediator 
masses has fo be re-evaluafed. The generafor level cross secfion 
fimes fhe accepfance af > 500 GeV for fhe model wifh cou¬ 
plings gq = gx = 1 / light Dark Matter of =io GeV and a Mmed 
=500 GeV scalar mediafor is af fhe order of 10 fb, i.e. jusf af fhe 
edge of fhe early Run-2 sensifivify. Increasing fhe mediafor mass fo 
1 TeV pushes fhe producf a x A down fo approximately 0.1 fb, be¬ 
low fhe LHC sensifivify. Therefore, we choose fo remove fhe 2 TeV 
mediafor mass from fhe grid and presenf fhe final grid wifh 33 
mass poinfs only, as shown in Tab. 2.5. One poinf af very high me¬ 
diafor mass (10 TeV) is added for each of fhe DM masses scanned, 
fo aid fhe reinferprefafion of resulfs in ferms of confacf inferacfion 
operators (EFTs). 
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Figure 2.21: Comparison of the Er 
distributions for the scalar and 
pseudoscalar models for different 
Atmed = 300 GeV and different Dark 
Matter masses. Ratios of the normal¬ 
ized distributions with respect to the 
first one are shown. A300 and A500 in 
the table denote the acceptance of the 
Bt > 300 GeV and Bt > 500 GeV cut, 
respectively. 


Table 2.5: Simplified model bench¬ 
marks for s-channel simplified models 
(spin-O mediators decaying to Dirac 
DM fermions in the scalar and pseu¬ 
doscalar case, taking the minimum 
width for gq = 1 and = 1) 


For fhe paramefer grid for scalar and pseudo-scalar mediafor 
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s-channel exchange, the Tmin/Mmed ratio is given in Tables 2.6 
and 2.7, respectively. In the on-shell regime, the ratio is between 
0.04 and 0.1. Very narrow resonances with Tmin/Mmed < 0.001 
correspond to the mass points where the mediator is off-shell. Nofe 
fhaf fhe loop-induced confribufion from gluons is ignored in fhe 
widfh calculation. 
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Table 2.6: Minimal width of the scalar 
mediator exchanged in s-channel 
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Table 2.7: Minimal width of the 
pseudo-scalar mediator exchanged 
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2 . 2.2 Additional considerations for V + Bt signatures 

The discussion of paramefers for fhe model wifh a color-singlef, 
spin -0 mediator parallels fhaf in Secfion 2. 

Even fhough fhe sensifivify of mono-boson searches fo fhis 
model is low and it may not be in reach of early LHC searches, 
this model can be generated for W, Z and photon searches in order 
fo reproduce fhe kinemafics of confacf interaction operators fhaf are 
further described in Section 3.2.1, to aid later reinterpretation. 

Other models of dark matter fhaf couple dominanfly to elec- 
troweak gauge bosons through either pseudo-scalar or vector medi¬ 
ators can be found in Ref. [LPS13]. 

2 . 2.3 Additional considerations for tt and bh+Br signatures 

Wifh fhe MFV assumpfion, fhe fop and bottom quark can play an 
imporfanf role in fhe phenomenology. The scalar and pseudoscalar 
mediafor models predicf nof only fhe monojef process described 
in Secfion 2.2, buf also production of Dark Matter in associafion 
wifh fop (or bottom) pairs, as illusfrafed in Fig. 2.22. Dedicated 
searches including jefs from heavy flavor quarks in fhe final sfafe 
can be designed for fhis signafure. Anofher class of simplified 
models, which includes a Dark Matter inferprefafion among many 
ofhers, and yields a single fop quark in fhe final sfafe, is defailed in 
Appendix A.i. 

In addition to the ff+DM models illustrated in Fig. 2.22, some 
theoretically motivated scenario (e.g. for high tanf in 2HDM in 
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^ Kb) 



Figure 2.22: Representative Feynman 
diagram showing the pair production 
of Dark Matter particles in association 
with ft (or bh). 


the pMSSM) privilege the coupling of spin -0 mediators to down 
generation quarks. This assumption motivates the study of final 
sfafes involving f;-quarks as a complemenfary search fo the ff+DM 
models, to directly probe the b-quark coupling. An example of such 
a model can be found in Ref. [BFG15] and can be obfained by re¬ 
placing top quarks with b quarks in Fig. 2.22. Note that, because 
of fhe kinemafics feafures of b quark production relative to heavy t 
quark production, a bb+UM final state may only yield one experi¬ 
mentally visible b quark, leading to a mono-b signature in a model 
that conserves b flavor. 

Dedicafed implementations of fhese models for fhe work of 
fhis Forum are available af LO+PS accuracy, even though the state 
of fhe arf is sef fo improve on a fimescale beyond thaf for early 
Run-2 DM searches as defailed in Secfion 4.1.5. The sfudies in fhis 
Section have been produced using a leading order UFO model 
within MadGrafh5_aMC@NLO 2.2.2 [Alw+14; All+14; Deg+12] 
using PYTHiA 8 for the parton shower. 


2.2.3.1 Parameter scan 

The parameter scan for fhe dedicafed tt+tj searches has been stud¬ 
ied in detail to target the production mechanism of DM associafed 
wifh heavy flavor quarks, and shares many defails of fhe scan for 
fhe scalar model wifh a gluon radiafion. The benchmark points 
scarming the model parameters have been selected to ensure that 
the kinematic features of the parameter space are sufficiently rep¬ 
resented. Detailed studies were performed fo idenfify poinfs in fhe 
m^^a, Sx' paramefer space fhaf differ significanfly 

from each other in terms of expected detector acceptance. Because 
missing transverse momentum is the key observable for searches, 
fhe mediator pj spectra is taken to represent the main kinemat¬ 
ics of a model. Anofher considerafion in determining fhe sef of 
benchmarks is fo focus on fhe parameter space where we expect 
the searches to be sensitive during the 2015 LHG run. Based on a 
projected integrated luminosity of 30 fb~^ expecfed for 2015, we 
disregard model poinfs wifh a cross section times branching ratio 
smaller than 0.1 fb, corresponding to a minimum of one expecfed 
event assuming a 0.1% efficiency times acceptance. 

The kinematics is most dependent on the masses m^^ and 
Figure 2.23 and 2.24 show typical dependencies for scalar and 
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pseudoscalar couplings respectively. Typically, the mediator pj 
spectrum broadens with larger The kinematics are also differ¬ 
ent between on-shell (Mmed > and off-shell (Mmed < 
mediafors as discussed in Secfion 2.2. Furfhermore, fhe kinemafic 
differences in fhe B j specfrum between scalar and pseudoscalar are 
larger for lighf mediafor masses wifh respecf fo heavier mediafors. 
If is fherefore imporfanf fo choose benchmark poinfs covering on- 
shell and off-shell mediafors wifh sufficienf granularify including 
fhe fransifion region befween on-shell and off-shell mediafors. 



Figure 2.23: Example of the depen¬ 
dence of the kinematics on the scalar 
mediator mass in the U+tj signature. 
The Dark Matter mass is fixed to be 
=lGeV. 



Ef®® [GeV] 


Figure 2.24: Example of the depen¬ 
dence of the kinematics on the pseu¬ 
doscalar mediator mass in the tt+tj. 
The Dark Matter mass is fixed to be 
=lGeV. All figures concerning the 
ti+Bj signature have been produced 
using a leading order model within 
MadGrafh5_aMC@NLO 2.2.2, 
using PYTfflA 8 for the parton shower. 


Typically only weak dependencies on couplings are observed 
(see Fig 2.25) where fhe variafion wifh widfh of fhe infegral over 
parfon disfribufions is unimporfanf. As shown in Secfion 2.1.1, 
for couplings ^ 0(1) the width is large enough that the pj of fhe 
mediafor is defermined mainly by fhe PDF. 

Af large mediafor masses (~ 1.5 TeV) or very small couplings 
(~ 10“^), widfh effecfs are significanf, buf fhese regimes have pro- 




































DARK MATTER BENCHMARK MODELS FOR EARLY LHC RUN-2 SEARCHES: 


REPORT OF THE ATLAS/CMS DARK MATTER FORUM 4I 

duction cross sections that are too small to be relevant for 30 fb~^ 
and are not studied here. However, with the full Run 2 dafasef, 
such models may be wifhin reach. 


Figure 2.25: Study of the depen¬ 
dence of kinematics on the width of 
a scalar mediator tt+Bj. The width is 
increased up to four times the minimal 
width for each mediator and Dark 
Matter mass combination. 
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Anofher case where fhe widfh can impacf fhe kinemafics is when 
nifp^a is slighfly larger fhan 2 m^. Here, fhe widfh defermines fhe 
relafive confribufion befween on-shell and off-shell mediators. An 
example is given in Fig. 2.26. As fhe minimal widfh choice pursued 
in fhis documenf is fhe mosf conservafive one, fhis effecf can be 
neglected in order fo reduce fhe number of benchmark poinfs fo be 
generated. 



Figure 2.26: Dependence of the kine¬ 
matics on the width of a scalar media¬ 
tor tt+Bj- The width is increased up to 
the mediator mass. Choices of media¬ 
tor and Dark Matter masses such that 
is slightly larger than 2 m^ is the 
only case that shows a sizeable varia¬ 
tion of the kinematics as a function of 
the width. 


The poinfs for fhe parameter scan chosen for fhis model are 
lisfed in Table 2.5, chosen fo be harmonized wifh fhose for ofher 
analyses employing fhe same scalar model as benchmark. Based on 
fhe sensifivify considerafions above, DM masses are only simulated 
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up to 500 GeV (but the 5 TeV mediator point is retained) leading to 
a total of 24 benchmark points. However for fhese searches we rec¬ 
ommend fo generafe and simulafe scalar and pseudoscalar models 
separafely, as fhe kinematics differs due fo fhe differenf coupling of 
fhe mediator fo fhe final sfafe fop quarks in fhe fwo cases, as shown 
in Figs. 2.23 and 2.24. 

Similar sfudies were performed in fhe bb case. If was found fhaf 
fhey show fhe same weak dependence of fhe kinematics of fhe 
evenf on fhe mediafor widfh. The same benchmark parameters of 
fhe tt case could fhen be chosen. 


2.3 Colored scalar mediator, t-channel exchange 

The preceding secfions address models wifh a Dirac fermion cou¬ 
pled fo fhe SM fhrough exchange of a neufral spin -0 or spin -1 par- 
ficle in an s-charmel process. A f-channel process may couple fhe 
SM and DM direcfly, leading fo a differenf phenomenology. For 
complefeness, we examine a model where x is a Sfandard Model 
(SM) single!, a Dirac fermion; fhe mediafing parficle, labeled (p, is 
a charged scalar color friplef and fhe SM parficle is a quark. Such 
models have been sfudied in Refs. [AWZrq; PVZrq; BB13; DiF+13; 
Cha+rq; Bel+r2]. However, fhese models have nof been sfudied as 
extensively as ofhers in fhis Forum. 

Following fhe example of Ref. [PVZrq], the interaction La- 
grangian is written as 


Ant — g {<P{i),LQ(i),L + <P{i),u,RU{i),R + <P{i),d,Rd{i),R)X (2.I2) 


!= 1,2 


where Q(!),L' “(0,R quarks of fhe ;-fh gen- 

erafion and ^{i),u,R 4 ’{i),d,R corresponding me¬ 

diators, which (unlike fhe s-charmel mediafors) musf be heavier 
than X- These mediators have SM gauge representations under 
{SU{ 3 ),SU{ 2 ))y of ( 3 , 2 )_i/ 6 , (3,1)273 and ( 3 , l)_i/3 respectively 
Variations of fhe model previously sfudied in fhe liferature include 
coupling fo fhe leff-handed quarks only [Cha+14; Bus+iqc], fo fhe 
i’{i),u,R [FliF+rq] or (p(i)^d,R [PVZrq; Abd+14], or some combina¬ 
tion [BB13; AWZrq]. 

The minimal width of each mediafor is expressed, using fhe 
example of decay fo an up quark, as 



(2.13) 
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which reduces to 


16^ V ^4-) 


(2.14) 


in the limit » mu^.y 

The generation index i for ^(,3 is linked to the incoming fermion(s), 
and it runs on all three quark generations due to the MFV assump¬ 
tion. Ref. [PVZ14] considers fwo exfreme cases for fhis model in 
ferms of cross-secfions: fhe case in which all mediafor flavors are 
presenf, leading fo fhe maximal cross-secfion, and fhe case in which 
only righf-handed down-f5rpe mediafors are presenf. Neifher of 
fhe models in fhis reference include couplings fo fhe fhird quark 
generafion, leading fo a violafion of fhe MFV assumpfion. In fhe 
case of purely down-f5rpe righf-handed squarks fhis is sfill safe 
from flavor consfrainfs. Furfhermore, reinfroducing fhe fhird gen¬ 
erafion squarks would lead fo models fhaf produce qualifafively 
similar signals in fhe mono-jef and SUSY squark searches, fhe main 
difference being fhe production cross-secfion. Af fhe same time 
fhe presence of fhird generafion squarks will lead fo furfher con¬ 
sfrainfs from ofher searches such as fhose for mono-bjefs, for stops 
and for sboffoms, as discussed in Sec. 2.3.3. The sfudies in fhis 
Section are performed using a model wifh a mediafor coupling fo 
all fhree generafion, following Ref. [Bel+12]. Furfher differences 
befween fhe fwo models (hypercharge, chiralify) only lead fo a 
change in fhe cross-secfion. The LO UFO model is inferfaced fo 
MadGraph5_aMC@NLO V2.2.3, buf if was nof possible fo go 
beyond parfon-level sfudies and interface fhose models fo a parfon 
shower in fime for fhe conclusion of fhis Forum. The sfafe of fhe arf 
for calculating fhese models is LO+PS, and fhe implemenfafion of 
mulfi-parfon merging has been sfudied in defail [Mal+15; Aqu+12; 
AVM09; PVZ14], and furfher sfudies should be underfaken prior fo 
generating signal samples for early Run-2 LHC searches. 

The leading-order processes involved in ^j+jef producfion are 
shown in Fig. 2.27. This model can also give a signal in fhe £7 + di- 
jef charmel when, for example, fhe x is exchanged in fhe f-channel 
and fhe resulfing (p pair each decay fo a jef + x- Fig- 2.28 shows 
fhe leading order diagrams. Excepf for fhe gg induced process, di- 
jef producfion fhrough fhe fhird-generafion mediafor 4 ’{ 3 ),u is riof 
possible, and producfion fhrough (p{x )4 is suppressed. However, 
if fhe coupling g includes a Yukawa coupling proportional fo fhe 
quark mass, and g is sufficienfly large, LHC searches will sfill be 
sensitive fo fhis model, as explained in Secfion 2.3.3. 

The diagram involving fhe f-channel exchange of x is sfrongly 
dependenf upon fhe Dirac fermion assumpfion. For a Majorana 
fermion, qq, qq, and qq producfion would be possible wifh fhe latter 
having a pronounced enhancemenf af fhe LHC. 

This model is similar fo fhe simplified model considered in SUSY 
searches, implemented as fhe MSSM wifh only lighf squarks and a 
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neutralino, except for two distinct points: the is a Dirac fermion 
and the coupling g is not limited to be weak scale {g ^ 1 ). In the 
MSSM, most of these processes are sub-dominant, even if reso¬ 
nantly enhanced, because the production is proportional to weak 
couplings. In the more general theories considered here, g is free 
to take on large values of order i or more, and thus diagrams ne¬ 
glected in MSSM simulation can occur at a much higher rate here. 
While constraints from SUSY jets+^j analyses on MSSM mod¬ 
els can be recast to apply to the specific model in this report, DM 
searches should also directly test their sensitivity to the MSSM 
benchmark models. 



Figure 2.27: Leading order mono-jet 
(-channel processes, adapted from 
[PVZ14]. 


The state of the art calculation for these models is LO and 
they can be interfaced with a parton shower program. The stud¬ 
ies in this Section use a LO model implementation within Mad- 
Grafh5_aMC@NLO V2.2.3, but no parton shower could be em¬ 
ployed in the time-frame of the conclusions of this Forum. Further 
implementation details can be found in Section 4.1.3. 
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Figure 2.28: Leading order two-jet 
(-channel processes, adapted from 
[PVZ14]. 
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2.3.1 Parameter scan 

As for the s-channel models, we adopt the simplifying assumption 
that the mediator masses and couplings are equal for each flavor 
and handedness. The free parameters are then 

{m^, M^, ^}. (2.15) 

Ref. [PVZ14] studies the parameter space and obtains bounds 
on this model from LHC Run-i mono-jet and dijets+^j data. The 
Forum did not exhaustively compare the kinematic distributions of 
the t-charmel models as done in the s-charmel case. In particular, 
the absence of a parton shower simulation can affect some of the 
conclusions on the points and sensitivity chosen. While this means 
the conclusions on the parameter scan below should be taken with 
more caution, the model is plausible and distinctive, and it should 
be included in the design of early Run-2 LHC searches. 

As in the s-channel models, scans should be performed over 
and M,p. The viable ranges of both parameters nearly coin¬ 
cide with the scan proposed for the s-channel. For the early Run-2 
searches, we recommend to generate and fully simulate a sub¬ 
set of the s-charmel mono-jet grid that accounts for the on-shell 
and off-shell regions. In contrast to the s-channel case, the bounds 
one obtains from $j+\ searches depend strongly on the width of 
the mediator, as is visible in Figs. 5 and 6 of Ref. [PVZ14] and in 
Fig. 2.29 (a), except in the heavy mediator limit (M^ w 2 TeV). This 
figure has been obtained applying a simplified analysis selection 
(cuts on the leading jet pj >150 GeV and rj < 2.8, £t> 150 GeV.) 
using MadAnalysis [Con+14; Dum+15]. Figure 2.29 (b) also shows 
that, if the DM mass is low and the mediator is produced on-shell 
and its width is narrow, the cross-section is dominated by qg qxx 
diagram. The mediator energy is then split evenly between the light 
DM particles and the quark, leading to a broad enhancement at 
^med/ 

Points with distinct kinematic distributions for a preliminary 
scan in { ntj^, M,p, g} are selected taking into account the expected 
sensitivity of Run-2 searches, and requiring at least 100 events to 
pass the kinematic cuts outlined for Fig. 2.29 in 25 fb~^ of collected 
data, and respect F/M^ed < 1 . They are outlined in Table 2.8. 

The conclusions in this table may change when a parton shower is 
employed together with multiparton matching. 

2.3.2 Additional considerations for V signatures 

The models and parameters with emission of an EW boson gener¬ 
ally follow those in Section 2.3. even though different diagrams are 
involved. A representative Fe5mman diagram can be constructed 
by replacing a final-state gluon in Fig. 2.27 with a 7, W, Z boson, 
but radiation of electroweak bosons directly from the mediator also 
leads to a mono-boson signature. 


Events (scaled to one) 
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(a) Hr distribution for a 200 GeV (-channel mediator, when varying 
the couplings. 


Figure 2.29: Kinematic distributions 
normalized to unit area from the 
(-channel model from Ref. [Bel+12], 
using MadAnalysis [CFS13; Con+14] 
and simplified analysis cuts on the 
leading jet pj >150 GeV and rj < 

2.8, Hr>i 50 GeV. For these models, 
a LO UFO model is interfaced to 
MadGraph5_aMC@NLO V2.2.3, and 
studies are at parton-level only. 



(b) Leading jet pj distribution for a 2 TeV (-channel mediator with 
small (^=0.5) to large (g=7) couplings with a DM mass of 1 GeV 
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tn ^/ GeV 

Mmed/GeV 

couplings 

1 

10 50 100 300 

0.1, 1, 3, 7 

1 

500 1000 

0.25, 1, 3, 7 

1 

2000 

h 3/7 

50 

55 

0.1, 1, 3, 471 

50 

200 300 

0.1,1, 3, 7 

500 

550 

1/3 

500 

1000 

0.25,1, 3 

500 

2000 

3 

1000 

1100 

3/471 

1000 

2000 

471 


Table 2.8: Simplified model benchmark 
points for f-channel simplified model 
(spin -0 mediators coupling to Dirac 
DM fermions, taking the minimum 
width.) 


The models considered in Section 2.3 present a relevant dif¬ 
ference concerning final sfafes wifh an elecfroweak boson. In fhe 
model in [Bel+12], bofh righf- and leff-handed mediafors can ra- 
diafe a Z boson, while only fhe leff-handed mediafor in [Bel+12] 
allows for W and Z radiafion. 

The sfudies in fhis Section use fhe LO+PS UFO model from [Bel+12] 
in MadGraph 5 _aMC@NLO V2.2.3, using pythia 8 for fhe parfon 
shower. Figure 2.30 shows fhe £7 disfribufion for fhe hadronic 
Z+^7 final sfafe, wifh varying DM and mediafor mass, before any 
selecfion. The accepfance for a series of basic analysis selections (£7 
>350 GeV, leading jef pj > 40 GeV, minimum azimufhal angle 
befween jef and £7 > 0.4) applied af fhe generafor level is shown in 
Figure 2.31. 



[GeV] 


Figure 2.30: Missing transverse mo¬ 
mentum distribution for the hadronic 
Z+^Bt final state, for the simplified 
model with a colored scalar mediator 
exchanged in the f-channel. 


The discussion of fhe paramefer scan for fhe f-charmel model 
in fhe case of signafures including EW bosons parallels fhaf of 
fhe monojef case for mediafor and DM masses, buf no kinemafic 
dependence on fhe widfh is observed, so a coupling scan is nof 
needed. 
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moM [GeV] 


Figure 2.31: Acceptance for the 
hadronic Z+Bt final state, for the 
simplified model with a colored scalar 
mediator exchanged in the t-channel. 


2.3.3 Additional considerations for signatures with b—quarks + 

Models of bottom-flavored Dark Matter that are closely related to 
the t-charmel mediated model from this Section have been pro¬ 
posed in Refs. [LKW13; Agr+i4b]. We describe the b-FDM model 
of Ref. [Agr+i4b], created to explain the Galactic Center (GC) 
gamma-ray excess observed in data collected by the Fermi-LAT 
collaboration [Day+14; CCW15]. This model favors couplings to 
third-generation quarks via Yukawa couplings, therefore respecting 
the MFV assumption. 

The model contains a Dirac fermion transforming as a flavor 
triplet, exclusively coupling to right-handed down-type quarks. The 
third component of the triplet Xb comprises the cosmological DM. 

Within the MFV framework, the other fermions in the flavor triplet 
can be made sufficiently heavy and weakly-coupled that they can 
be neglected in the analysis. A flavor singlet, color triplet scalar 
field O mediates the interactions between the DM and the Standard 
Model quarks. The model is similar to the MSSM with a light bot¬ 
tom squark and neutralino, and is thus a flavor-specific example 
of a t-charmel model. Similar top-flavored models can exist, as e.g. 
in Refs. [KT13; BLWiqa]. In the case where the top coupling is the 
main DM coupling, the signal is very similar to a signal from a stop 
quark, since unlike the other t-charmel cases there is no top in the 
initial state parton distribution functions (PDFs). This is the reason 
why it wasn't considered as an additional model. More recent lit¬ 
erature shows that other flavor states could also contribute to LHG 
signals, as shown in Ref. [KKY15], but such models will have to be 
investigate on a longer timescale with respect to that of this Forum. 

The Lagrangian considered is given by 

-L D g^*XbbR + h.c. (2.16) 

This model is known at LO+PS accuracy, and the studies in this 
Section use a LO model implementation within MadGraph 5 _aMG@NLO 
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V2.2.3 interfaced to pythia 8 for the parton shower. Further imple¬ 
mentation details can be found in Section 4.1.5. 


2.3.3.1 Parameter scan 

In this model, the interference of diagrams with QCD production 
of the mediator (which scale as g^) with diagrams that are propor¬ 
tional to the coupling g in the b+tj and bb+tj final states. In the 
case of large couplings, this is not conducive to a simple scaling 
behavior that would allow us to reduce the number of points to be 
simulated. This can be seen in Fig. 2.33. 

A full study of the parameter scan for this model was not avail¬ 
able for this report; thus for early Run-2 searches we recommend 
scanning a range of possible widths as discussed in a more limited 
way than for the t-channel mono-jet, spanning from the minimal 
width to a value approaching the particle limit, e.g. g = 0 . 5 , 1 , 2 , 3 . 

A coupling benchmark such as ^ = 1 should be considered for each 
mass point since this would be a distinctive feature of this bench¬ 
mark from SUSY models with sbottom squarks (see Section 2.3 for 
further discussion). 

A scan of Dark Matter and mediator masses should be done 
in the on-shell region Mo > m^ + mi,, since the cross-sections 
in the off-shell region are too small to be probed with early LHC 
data, sparming from 10 to 500 GeV in m^^ and from 10 to 1300 GeV 
in Mo- Examples of the kinematic distributions produced by this 

model are shown in Fig. 2.32 5. ^poUowing the grounding assump¬ 

tions in this report, the normalization 
tn khp rplir dpnsitv is rnrisidprpd nnlv 




2.4 Spin -2 mediator 

In models with extra dimensions, the Kaluza-Klein excitations of 
the graviton could also serve as a mediator between the Standard 
Model and dark sector physics. This kind of model was not studied 
in the forum and is not included in the recommendations, but mod¬ 
els such as Ref. [LPSiqa; LPSiqb] may warrant further study on a 
longer timescale. 


Figure 2.32: t.j (left) and jet multiplic¬ 
ity (right) for various DM and media¬ 
tor masses and couplings normalized 
to the relic density observed in the 
early universe. Studies in this section 
use a LO UFO model implementation 
within MadGraph5_aMC@NLO 
V2.2.3 interfaced to pythia 8 for the 
parton shower. 
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Figure 2.33: £7- (left) and jet multi¬ 
plicity (right) for = 35 GeV and 

2.5 Presentation of results for reinterpretation of s-channel me- alo = soo Gev for varying couplings 
diator models 

The aim of the parameter grid optimization done for the s-charmel 
models in the previous sections is to reduce the parameter space 
that must be simulated. We then need a procedure for populating 
the full parameter space by using the simulated grid points. We 
recommend doing this as follows: 

• When the dependences on parameters are known, the cross 
sections and efficiencies at general points can be calculated from 
the grid data. 

• In other cases, this information can be obtained by interpolation 
between the grid points. We have chosen the grid points so that 
the dependence is sufficiently smooth that this will be possible. 

The results of the scan over the couplings presented in the previ¬ 
ous sections indicate that there are no changes in kinematic distri¬ 
butions for different choices of the coupling strengths. This means 
that the acceptance remains the same in the whole gq-gx plane and 
it is sufficient to perform the detector simulation only for one sin¬ 
gle choice of gq,gx- The resulting truth-level selection acceptance 
and the detector reconstruction efficiency can then be applied to all 
remaining grid points in the gq-gx plane where only the generator- 
level cross section needs to be known. This significantly reduces the 
computing time as the detector response is by far the most CPU¬ 
intensive part of the Monte Carlo sample production. However, the 
number of generated samples can be reduced even further if a pa¬ 
rameterization of the cross section dependence from one grid point 
to another exists. In this section, we describe the details of a cross 
section scaling procedure that can be used to reinterpret results for 
a fixed coupling for s-channel mediator models. The studies in this 
section employ the powheg [HR15] generator. 

The propagator for the s-charmel exchange is written in a Breit- 
1 

Wiener form as ^ -- ——where q is the momentum 
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transfer calculated from the two partons entering the hard process 
after the initial state radiation, which is equivalent to the momen¬ 
tum of the Dark Matter pair The size of the momentum transfer 
with respect to the mediator mass allows us to identify three cases: 

• off-shell mediator, when ^ ^med leading to suppressed cross 
sections, 

• on-shell mediator, when ~ ^med leading to enhanced cross 
sections, 

• effective field theory (EFT) limit when 

In the case of the off-shell mediator and the EFT limit, the first 
and second term in the propagator dominate, respectively, which 
reduces the dependence on the mediator width. Therefore, in these 
cases one can approximate the cross section as 


^ Using a running width and replacing 
the denominator of the propagator 
with + ' Q" mL 

be considered in the case of wide 
mediators [Bar+89]. 


« g^g^x- ( 2 . 17 ) 

The on-shell regime is the most interesting one as it gives the best 
chances for a discovery at the LHC given the cross section enhance¬ 
ment. The propagator term with the width carmot be neglected in 
this case and, in the narrow width approximation which requires 
T ^ Mined (this is not necessarily the case in the benchmarks con¬ 
sidered in the scans), one can integrate 


/ 



n 

Mined r 


which further implies the cross section scaling 


(2.18) 


cr <x 


2 2 

Sqgl 

T 


(2.19) 


The narrow width approximation is important here as it ensures 
an integration over parton distribution functions (PDFs) can be 
neglected. In other words, it is assumed the integrand in Eq. 2.18 
is non-zero only for a small region of s, such that the PDFs can be 
taken to be constant in this range. By simplifying the dependence 
of the minimal width on the couplings as T ~ -h g^, one can 
approximate this scaling rule in the extreme cases as follows 


a oc 


cr <x 


Sqgx 

Sq+gx 

^qSx 

Sq+gx 


Sq^Sx 


Sq^Sx 


>4 

(2.20) 


(2.21) 


However, it is important to keep in mind that this formula omits 
color and multiplicity factors as well as possible Yukawa suppres¬ 
sion, and there is no simple scaling rule for how the cross section 
changes with the Dark Matter mass and the mediator mass, or for 
mediators with a large width, because PDFs matter in such cases as 
well. Therefore, the scaling procedure outlined above is expected to 
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work only for fixed masses and fixed mediator widfh, assuming fhe 
narrow widfh approximafion applies. 

Figure 2.34 shows fhe minimal widfh over fhe mediator mass 
in fhe gc^gx plane for vecfor and scalar mediafors for M^ed = 

100 GeV and looo GeV, faking = 10 GeV. The individual 
colors indicafe fhe lines of consfanf widfh, along which fhe cross 
secfion scaling may work for narrow mediafors. The limifing case 
Tmin = Mined defines fhe upper values of fhe couplings below 
which fhe narrow widfh approximafion can be considered and 
provides more sfringenf consfrainf fhan fhe perfurbafive limif = 
g-^ = 4/1. For vecfor and axial-vecfor mediafors, fhe minimal widfh 
is predominanfly defined by due to fhe number of quark flavors 
and fhe color factor. On fhe confrary, bofh fhe Sfandard Model and 
Dark Matter parfial widfh have comparable confribufions in case of 
scalar and pseudo-scalar mediafors if fhe fop quark charmel is open 
(Mined > 2 OTt). However, mosfly g^ defines fhe minimal widfh for 
Mmed < due to fhe Yukawa-suppressed lighf quark couplings. 



The performance of fhe cross secfion scaling is demonsfrafed 
in Fig. 2.35 where two mass poinfs Mmed = 100 GeV and 1 TeV 
wifh = 10 GeV are chosen and rescaled from fhe sfarfing poinf 
gq = gx = ^ according to Eq. 2.19 to populafe fhe whole gq-gx 
plane. This means fhe widfh is nof kepf consfanf in fhis fesf and 
fhis is done in purpose in order to poinf ouf deviations from fhe 
scaling when fhe widfh is altered. For each mass poinf, fhe rescaled 


Figure 2.34: Minimal width over the 
mediator mass for vector (top) and 
scalar (bottom) mediators as a function 
of the individual couplings gq and 
assuming Mmed = 100 GeV (left) and 
Mmed = 1 TeV (right), = 10 GeV is 
considered in all cases. Only the cases 
with Tmin < Mmed ^re shown. 
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cross section is compared to the generator cross section and the 
ratio of the two is plotted. For the given choice of fhe mass poinfs, 
fhe scaling seems fo work approximafely wifhin fhe precision of 
~ 20 % in fhe region where Tmin < Mmed- Consfanf colors indicafe 
fhe lines along which fhe cross secfion scaling works precisely and 
fhere is a remarkable resemblance of fhe pafferns shown in fhe 
plofs of fhe mediafor widfh. To prove fhe scaling along fhe lines 
of consfanf widfh works, one such line is chosen in Fig. 2.36 for a 
scalar mediafor, defined by M^ed = 300 GeV, = 100 GeV, 

= gx = 1/ and fhe rescaled and generafed cross secfions are 
found fo agree wifhin 3%. 


Q 

O) 



Figure 2 . 35 : Ratio of the rescaled and 
^ g generated cross sections in the gq-gx 

plane. The point at g^ = g^ = 1 , 

1 .4 taken as a reference for the rescaling, 

is denoted by a star symbol. Scalar 
■0 model with M^ed = 100 GeV (left) 

"I 2. srid 1 TeV (right) is plotted for = 

10 GeV. The limiting case Tmin = 

' ■ ' Aimed is indicated by a black line and 

-j no results are shown beyond. 

0.9 

0.8 

0.7 
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^SM 
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2.5.1 Proposed parameter grid for cross-section scaling 

We propose fo deliver collider resulfs in fhe gq-gx plane using 
the following prescripfion, fo ease reinterprefation fhrough cross- 
secfion scaling: 




















wag 
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gSM 


Figure 2 . 36 : Scaling along the lines of 
constant width. The line of constant 
width for M^ed = 300 GeV and = 
100 GeV, intercepting gq = = 4 

is shown on left. The generated and 
rescaled cross sections are compared in 
the middle, the corresponding ratio is 
shown on right. 
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• Since the shapes of kinematic quantities do not change for differ¬ 
ent couplings, use the acceptance and efficiency for the avail¬ 
able = 50 GeV, M^ed = 300 GeV grid point from the 

plsriE for the scalar and pseudo-scalar mediator. In 
case of the vector and axial-vector mediator, use the grid point 
= 150 GeV, M^ed = 1 TeV. 

• Generate additional samples in order to get generator cross 
sections only. For scalar and pseudo-scalar mediator, choose 

= 50 GeV, Mmed = 300 GeV with the following values for 
yq = g^: 0.1,1, 2, 3. For vector and axial vector mediator, choose 
= 150 GeV, M^ed = 1 TeV with the following values for 
yq = g^: 0.1, 0.25, 0.5, 0.75,1,1.25,1.5. The upper values are 
defined by the minimal width reaching the mediator mass. 


• Rescale the generator cross sections for on-shell resonance pro¬ 
duction along the lines of constant width in order to populate 
the whole gq-gx plane in the region Fmin < Mmed- The scaling 
follows from Eq. 2.19 which for the constant width implies: 


a' = cr X 


SqSx 

s\sl ■ 


(2.22) 


2.5.2 Rescaling to different mediator width 

In general it is also important to consider a larger mediator width 
than Fmin in order to accommodate additional interactions of the 
mediator with the visible and hidden sector particles [BFG15; 
FIar+15]. If the narrow width approximation applies, the cross 
section scaling method described above can be used to reinterpret 
the results presented for the minimal width, since multiplying the 
width by factor n is equivalent to changing the coupling strength by 
factor i.e. 

The cross section for the sample with couplings yq and g^ and 
modified mediator width F = nFmin can therefore be rescaled from 
a sample generated with the minimal width corresponding to the 
couplings scaled by as described in the following formula. 


Sx' ”Tmin (S'q' Sx ) ) 


1 


(^{Vngq, Vngx, Tmin {Vngq,Vngx)) 


(2.24) 

The advantage of doing this is in the fact that no event selection 
and detector response needs to be simulated since the changes in 
couplings do not have an effect on the shapes of kinematic distribu¬ 
tions. 

It should be noted again that this procedure is only useful when 
the narrow width approximation applies. Care must be taken to 
ensure that is the case. For example, in the vector and axial-vector 
cases, one quickly breaks this approximation even for small n. 
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2.5.3 Additional considerations for tt and bb+tj signatures 

The cross-section scaling considerations shown in Sec. 2.5 still apply 
for the reactions in the scalar and psuedoscalar models with explicit 
b and t quarks. Here we detail the specific sfudies done for fhe ti 
model. 

Given fhaf fhe kinemafics are similar for all couplings g 1 , 
we recommend fo generafe only samples wifh = 1- It 

follows from fhis fhaf fhese benchmark poinfs should be a good 
approximafion for non-unify couplings and for g^ f gq, provided 
fhaf fhe sample is rescaled fo fhe appropriafe cross secfion times 
branching ratio. 

While fhe simple scaling function 

cr' X BR' = [ct X BR] x (— ^ x (— ^ x (2.25) 

is sufficienf for a limifed range of coupling values (see Fig. 2.37 for 
example), fhis scaling is only approximafe (up fo 20%) and relies on 
fhe narrow widfh approximafion, ignoring PDFs effecfs. 


M0 = 400, Mx = 1, S 




Figure 2 . 37 : An example compar¬ 
ing a simple cross section scaling 
versus the computation from the 
MadGraph 5 _aMC@NLO gener¬ 
ator, for a scalar tt+Bj model with 
= 400 GeV, = 1 GeV and all 

couplings set to unity. In this exam¬ 
ple, the scaling relationship holds for 
Tf/m^ below 0 . 2 , beyond which finite 
width effects become important and 
the simple scaling breaks down. 


I'.M« 
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Specific models for signatures with EW bosons 


In this Section, we consider specific models with a photon, a W bo¬ 
son, a Z boson or a Higgs boson in the final sfafe (V+Et signafure), 
accompanied by Dark Matter parficles fhat eifher couple direcfly fo 
fhe boson or are mediafed by a new particle. The common feafure 
of fhose models is fhaf fhey provide differenf kinemafic disfribu- 
fions wifh respecf fo fhe models described in Secfion 2. 



Figure 3 . 1 : Sketch of benchmark 
models including a contact interac¬ 
tion for V+MET searches, adapted 
from [Nett 14 ]. 


The models considered in fhis Secfion can be divided info two 

cafegories: 

V-specific simplified models These models posfulafe direcf couplings 
of new mediators fo bosons, e.g. fhey couple fhe Higgs boson fo 
a new vector or fo a new scalar [Car+14; BLWi4b]. 

Models involving a SM singlet operator including a boson pair that couples to Dark Matter through a contact interaction 
Shown on fhe righf-hand side of Figure 3.1, fhese models allow 
for a confacf interaction verfex fhaf direcfly couples fhe boson fo 
Dark Matter [Cof+13; Car+13; CHH15; BLWi4b]. These models 
are included in fhis reporf devofed fo simplified models since 
UV complefions for mosf of fhese operafors proceed fhrough 
loops and are nof available fo dafe. These models provide a 
benchmark fo mofivafe signal regions fhaf are unique fo searches 
wifh EW final sfafes and would ofherwise nof be sfudied. How¬ 
ever, we recommend fo use fhese models as placeholders and 
emphasize model-independenf resulfs especially in signal re¬ 
gions tailored to these models. Wherever results are interpreted 
in terms of fhese operafors, a fruncation procedure fo ensure fhe 
validify of fhe EFT should be employed, as defatted in fhe nexf 
Secfion (Sec. 5). 
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The following Sections describe the models within these cate¬ 
gories, the parameters for each of the benchmark models chosen, 
the studies towards the choices of the parameters to be scanned. 

3.1 Specific simplified models including EW bosons, tailored to 
Higgs+MET searches 

Three benchmark simplified models [Car+14; BLWi4b] are recom¬ 
mended for Higgs+^T searches: 

• A model where a vector mediator (Zg) is exchanged in the 
s-channel, radiates a Higgs boson, and decays into two DM par¬ 
ticles (Fig. 3.2 (a)). As in Section 2.1, we conservatively omit 
couplings of the Zg to leptons. 

• A model where a scalar mediator S is emitted from the Fliggs 
boson and decays to a pair of DM particles (Fig. 3.3). 

• A model where a vector Z' is produced resonantly and decays 
into a Higgs boson plus an intermediate heavy pseudoscalar 
particle A^, in turn decaying into two DM particles (Fig. 3.2 (b)). 



Figure 3.2: Examples of Feynman 


diagrams leading to Higgs+Er events: 
(a) a model with a vector mediator (Z') 
coupling with DM and with the Higgs 
boson h, and (b) a 2HDM model with 
a new invisibly decaying pseudoscalar 
from the decay of an on-shell 


resonance Z' giving rise to a Higgs+Er 
signature . 


These models are kinematically distinct from one another, as 
shown in the comparison of the spectra in Fig. 3.4 for high and 
low masses of the pseudoscalar mediator. Figure 3.4 (a) shows the 
Bt distribution for models with high mediator masses (jng = 1 TeV, 
niz' = 1 TeV, ?n^o = 1 TeV) and DM mass of either 50 (Zg and 
models) or 65 GeV (scalar mediator model). Figure 3.4 (b) shows 
the B T distribution for models with low pseudoscalar mediator 
masses = 1^0 GeV, mz' = 1 TeV, = 100 GeV) and DM 
mass of 1 TeV for all models. 

Predictions for this class of models have been so far considered 
at LO+PS, even though they could be extended to NLO+PS in the 
near future. The studies in this Section have been performed using 
a model within MadGraph5_aMG@NLO V2.2.3, interfaced to 
PYTHiA 8 for the parton shower. The implementation details for 
these models are discussed in Section 4.2.1.2. 
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Figure 3.3: Examples of Feynman 
diagrams leading to Higgs+^r events 
for a model with a scalar mediator (S) 
coupling with DM and with the Higgs 
boson h. 



3.1.1 tj+Higgs from a haryonic Z' 

The model shown in Fig. 3.2 (a) postulates a new gauge boson 
Z' corresponding to a new LI( 1 )b baryon number symmetry. The 
stable baryonic states included in this model are the DM candidate 
particles. The mass of the Z' boson is acquired through a baryonic 
Higgs hg, which mixes with the SM Higgs boson. 

The interactions between the Z', the quarks and the DM are 
described by the following Lagrangian: 


^ + g^x 7 ^‘xzf (3.1) 

The quark couplings are fixed fo be equal fo one fhird of fhe 
gauge coupling gs, while fhe DM coupling fo fhe Z' are propor- 
fional fo fhe baryon number and fo fhe gauge coupling (g^ = Bgs). 
No lepfonic couplings of fhe Z' are allowed, fhus evading dilep- 
fon consfrainfs. Affer incorporafing fhe mixing of fhe baryonic 
and SM Higgs bosons, fhis model is is described by fhe following 
Lagrangian ferm af energies below mz' 

^eff = , (3.2) 

The firsf ferm of fhis equafion is fhe sfandard DMV model in fhe 
large M^/ limif. This ferm can lead fo a monojef signafure, which 
can be also used fo consfrain fhis model. The second ferm describes 
fhe inferacfion between the Z' and the SM Higgs boson, via the 
coupling gfjz'z' — where sin 9 is the mixing angle between 

the SM Higgs and the baryonic Higgs hg, and vg is the Baryonic 
Higgs vacuum expectation value. 

In its most general form, fhis model can contribufe fo mono-Z 
signals due fo fhe Z' mixing wifh fhe Z or phofon. Nofe fhaf EWSB 


’ The operator in Eqn. 3.2 is an effec¬ 
tive one, to highlight the two main 
terms. The full dimension-4 simplified 
model is used in the model for event 
generation. 







Events/20 GeV Events/20 GeV 
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Figure 3.4: Comparison of the missing 
transverse momentum distributions 
at generator level in different simpli¬ 
fied models leading to a Higgs+^r 
signature. The model parameter set¬ 
tings are detailed in the text. The 
figures in this Section have been ob¬ 
tained using LO UFO models within 
MadGrafh5_aMC@NLO V2.2.3, 
interfaced to pythia 8 for the parton 
shower. 


(a) High mediator mass 



(b) Low mediator mass 
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and LI( 1 )b breaking do not lead to this mixing at tree-level. In¬ 
stead, kinetic mixing occurs between the 11 ( 1 ) y and LI( 1 )b gauge 
bosons due to the gauge invariant term Fy^FBfjv This mixing is a 
free parameter which we assume to be small in order to focus on 
fhe mono-Higgs signafure. Mixing may also occur due fo radia- 
five corrections, however fhis is model dependenf so we choose fo 
ignore fhis here. 

The predictions of fhe model depend upon fhe fwo additional 
paramefers beyond an s-charmel simplified model, namely fhe 
mixing angle between baryonic Higgs /ig and fhe SM-like Higgs 
boson sin 0 and the coupling of fhe mediator fo SM-like Higgs 
boson, ghz'Z'- Thus, a full model is specified by: 

{Mmed/ g^, gq, sin0, g^z'Z'} ■ (3-3) 


■j.i.i.i Parameter scan 

The widfh of fhe Z' mediator is calculated using all possible decays 
fo SM particles (quarks) and to pairs of DM particles if kinemafi- 
cally allowed as in the DMV model. 

The dependence of fhe missing fransverse momenfum (£7) on 
fhe model paramefers is sfudied by varying fhe paramefers one af 
a fime. The variafion of paramefers ofher fhan Mmed 3ud does 
nof resulf in significanf variations of fhe £ j specfrum, as shown in 
Figures 3.5. Figure 3.6 shows fhaf for an on-shell mediafor, varying 
m-z wifh fhe ofher paramefers fixed does nof affecf fhe £7 disfri- 
bufion, while fhe disfribufion broadens significanfly in fhe case of 
an off-shell mediafor. For fhis reason, fhe same grid in M^ed/ 
as for fhe vecfor mediafor of fhe jef+£7 search (Table 2.1) is chosen 
as a sfarfing poinf. The coupling ghz'Z'i along wifh gq and g^, are 
subjecf fo perfurbativify bounds: 

<47T 

and 


ghZ'z' < V^mz' sine 

The value ghz'Z' / ^z' = 1 is chosen as a benchmark value for fhe 
generation of Monte Carlo samples since if maximizes fhe cross 
secfion (as shown in fhe following paragraph) wifhouf violafing 
fhe bounds. The mediafor-DM coupling g^ is fixed fo 1, and fhe 
mediafor-quark gq coupling is fixed fo 1/3. The kinematic disfri- 
bufions do nof change as a funcfion of fhese paramefers, so resulfs 
for ofher values of ghZ'Z'/^z'> gx obfained fhrough 

rescaling by fhe appropriate cross sections. 

Figs 3.7 and 3.8 show fhe kinematic disfribufions for fhe fwo 
leading jefs in fhe H ^ bb decay channel, for fwo values of fhe 
mediafor mass and varying fhe DM mass. 

Analyses should perform furfher sfudies, beyond fhose sfudies 
performed for the forum, fo estimate fhe reach of fhe analysis wifh 
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respect to all points in the grid and therefore decide on a smaller 
set of grid poinfs fo be generafed. 




Figure 3.5: Missing transverse momen¬ 
tum distributions at generator level 
in the vector mediator scenario for 
different values of: the mediator-dark 
matter coupling (left), and the cou¬ 
pling between the mediator and the 
SM-like Higgs boson, scaled by the 
mediator mass, ghz'z' (right). 


3.1.2 Bt +Higgs from a scalar mediator 

A real scalar single! S coupling fo DM can be infroduced as a porfal 
befween SM and fhe dark secfor fhrough fhe Higgs field. The mosf 
general scalar pofenfial is defailed in Ref. [ORMW07], including 
ferms fhaf break Z2. The Z2 symmefry, which causes fhe new 
scalar fo also be a DM candidafe, is nof covered in fhis reporf, buf 
follows Ref. [Car-H4] inf reducing an addifional coupling fo DM 
fhaf breaks Z2 and leads fo a new invisible decay of S. For fhis 
reason, no symmefry is broken and no new inferacfions arise, so 
there is no dependence on the vacuum expectation value of S: a 
shift in the field leads fo a redefinition of fhe model couplings. 

The new scalar S mixes wifh fhe SM Higgs boson, and couples fo 
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Figure 3.6: Missing transverse momen¬ 
tum distributions at generator level 
in the vector mediator scenario: for 
different values of the dark matter 
mass and a mediator mass of Mmed 
= 100 GeV (left) and Mmed = 1 TeV 
(right). 
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Figure 3.7: Comparison of the kine¬ 
matic distributions for the two leading 
fo—jets (from the Higgs decay) in the 
vector Z' simplified model, when 
fixing the Z' mass to 100 GeV and 
varying the DM mass. 




(c) Angular distance between the two leading b—jets 
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Figure 3.8: Comparison of the kine¬ 
matic distributions for the two leading 
jets from the Higgs decay in the vector 
Z' simplified model, when fixing the 
Z' mass to 1000 GeV and varying the 
DM mass. 




(c) Angular separation of the two leading fo-jets 
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DM through a Yukawa term y^. The relevant terms in the scalar 
potential are: 


V D a\H\^S + b\H\^S^ + 

— la{h + i;)^S +lb{h + v)^S^ + ^{h + (3.4) 

where a, b are new physics couplings and A;, is the Higgs quartic 
coupling. 

The additional Lagrangian terms for this model are: 

fH 

£ D —y^XXicosB S — sinOh) — ■^qq{cos 6 h + sm 6 S) (3.5) 

where 9 is the mixing angle between the Higgs boson and the new 
scalar. 

Mono-Higgs signals in this second model arise through pro¬ 
cesses shown in Fig. 3.3 (a,b), or through the radiation of a Higgs 
boson from fhe t quark in fhe producfion loop, in Fig. 3.3 (c). 

The firsf two processes depend on fhe h^S and hS^ cubic terms 
in Eq. (3.4). At leading order in sinF, these terms are: 

Kubic- ^^i 2 ml + ml)h^S + bvhS^ + ... (3.6) 

with a and A;, expressed in terms of sin 9 and respecfively. At 
leading order of sinF, fhe h^S ferm is fixed once fhe mass eigen¬ 
values nif,, nig and mixing angle are specified. The h ferm is nof 
fixed and remains a free paramefer of fhe model, depending on fhe 
new physics coupling b. 

This model also has mono-X signafures fhrough h/S mixing. 

This model is relafed fo fhe scalar model discussed in Sec. 2.2 in fhe 
case of ntg ntf, or m^i :$> nig and Mmed equal fo fhe lighfer of fhe 
fwo masses, albeif wifh differenf mono-Higgs signafures due fo fhe 
hS^ verfex. 

3.1.2.1 Parameter scan 

The model is described by five paramefers: 

1. fhe Yukawa coupling of heavy scalar fo dark maffer, (also 
referred fo as y^) 

2. the mixing angle between heavy scalar and SM-like Higgs bo¬ 
son, sinF; 

3. the new physics coupling, F; 

4. mass of heavy scalar, mg, also fermed Mmedi 

5. mass of dark maffer. m^; 

The mixing angle is consfrained from currenf Higgs data to sat¬ 
isfy cos F = 1 within 10% and therefore sinF < 0 . 4 . This provides a 
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Figure 3.9: Missing transverse momen¬ 
tum distributions at generator level 
in the scalar mediator scenario, for 
different values of: the new physics 
coupling gh (left), and the mediator- 
dark matter coupling (right). 
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Starting point for the parameter scan in this model: we recommend 
to set sin 0 = 0 . 3 . 

Figure 3.10 shows that there is no dependence of fhe kinemafics 
from fhe value of fhis angle, and differenf values can be obfained 
via rescaling fhe resulfs for fhis mixing angle according fo fhe rele- 
vanf cross-secfion. If can also be observed from Figures 3.11 and 3.9 
fhaf fhe kinemafics of fhis model follows fhaf of fhe equivalenf 
jef+^T model: only small changes are observed in fhe on-shell re¬ 
gion, while fhe relevanf disfribufions diverge when fhe mediafor 
is off-shell. For fhis reason, fhe same grid in Mj^ed/ 
scalar mediafor of fhe jef+^j search (Table 2.5) is chosen as a sfarf- 
ing poinf. The Yukawa coupling fo DM i/dm is sef fo 1, fhe new 
physics coupling between scalar and SM Higgs b = 3. Resulfs for 
ofher values can be obfained via a rescaling of fhe resulfs for fhese 
paramefers. 
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Figure 3.10: Missing transverse mo¬ 
mentum distributions at generator 
level in the scalar mediator scenario: 
for different values of the mixing angle 
sin0. 


Figs. 3.12 and 3.13 show fhe kinemafic disfribufions for fhe two 
leading jefs in fhe H ^ bb decay channel, for two values of fhe 
mediafor mass and varying fhe DM mass. 

3.1.3 Higgs+Br signal from zHDM model with a Z' and a new pseu¬ 
doscalar 

In fhis simplified model [BLWiqb], a new Z' resonance decays fo 
a Higgs boson h plus a heavy pseudoscalar sfafe in fhe 2HDM 
framework, which in furn decays fo a DM pair. This model is repre- 
senfed in fhe diagram in Fig. 3.2 (b). 

The mofivafion for coupling fhe dark maffer fo fhe pseudoscalar 
is fhaf dark maffer coupling fo a Higgs or Z' boson is generically 
consfrained by ofher signal channels and direcf defecfion. A reason 
fo consider fhis model is fhaf if has differenf kinemafics due fo 
fhe on-shell Z' producfion, where for heavy Z' masses fhe Bj and 
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Figure 3.11: Missing transverse mo¬ 
mentum distributions at generator 
level in the scalar mediator sce¬ 
nario: for different values of the 
dark matter mass and a mediator 
mass of Mmed = 100 GeV (left) and 
M^ed = 1 TeV (right). 
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Figure 3.12: Comparison of the kine¬ 
matic distributions for the two leading 
jets from the Higgs decay in the scalar 
simplified model, when fixing the new 
scalar mass to 100 GeV and varying 
the DM mass. 




(c) Angular distance between the two leading b—jets 
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Figure 3.13: Comparison of the kine¬ 
matic distributions for the two leading 
jets from the Higgs decay in the scalar 
simplified model, when fixing the new 
scalar mass to 1000 GeV and varying 
the DM mass. 




(c) Angular distance between the two leading b—jets 
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Pj spectra are much harder. This model can satisfy electroweak 
precision tests and constraints from dijef resonance searches, and 
sfill give a pofenfially observable Higgs+^j signal. 

This model comprises two doublefs, where Oj, couples fo up- 
fype quarks and couples fo down-f5q)e quarks and lepfons: 

-£ D yuQ^uU + PdQ^dd + + h.c. (3.7) 

Affer elecfroweak symmefry breaking, fhe Higgs doublefs affain 
vacuum expecfafion values Vu and Vd, and in unifary gauge fhe 
doublefs are paramefrized as 


^d = ^( -sin/5H+ 

y /2 \vd — sina. h + cos a. H — i sin /3 

cos ^ H+ 

Vu + COS a.h + sin a H + i cos ji 

where h, H are neufral CP-even scalars, is a charged scalar, and 
is a neufral CP-odd scalar. In fhis framework, tan ji = Vu/Vd, 
and cc is fhe mixing angle fhaf diagonalizes fhe h — H mass squared 
mafrix. This model also confains an additional scalar single! cp fhaf 
leads fo spontaneous symmetry breaking. We take a = /3 — nil, 
in the limit where h. has SM-like couplings to fermions and gauge 
bosons as per Ref. [CGT13], and fan/I > 0.3 as implied from fhe 
perturbafivify of fhe fop Yukawa coupling. The Higgs vacuum 
expecfafion values lead to Z — Z' mass mixing, wifh a small mixing 
paramefer given by 





1 


e = 


SSz 


M%, - My 2 cos 9 i 


-{Zdvl+Zuvl) 


m|, — m| 


2 gz cos Bu 


g 


■Zu sin^ ji, 


( 3 - 9 ) 


where Z; are fhe Z' charges of fhe two Higgs doublets, and g and gz 
related to the mass-squared values in absence of mixing (M^)^ = 
^^(^d + ^D/( 4 cos 2 0,„) and = gl{zjvj + zlvl+zlvl). 

The production cross section for this model scales as (gz)^, as the 
decay width for this process to leading order in e (Eq. 3.9) is 

Tz'^hAO = UzSin^cos^)2^j^. (3.10) 


where the center of mass momentum for the decay products \p\ = 
2^^ yj (M|, - (m;, -h?n^o) 2 )(M|, - (m,, The Z' can also 

decay to Zh, leading to the same signature if the Z decays invisibly. 
The partial width for this decay is: 


Tz'^/iZ 


sin ^2)2 


IpI 

247 r 



(3.11) 


. We recommend to generate these two decays separately and com¬ 
bine them at a later stage. 
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3.1.3.1 Parameter scan 

The model is described by five parameters: 

• the pseudoscalar mass M^o, 

• the DM mass 

• the Z' mass, Mz', 

• isn^{=Vu/vd), 

• the Z' coupling strength g^. 

To study the signal production and kinematic dependencies on 
these parameters, we produced signal samples varying each of the 
five parameters through MadGraph5_aMC@NLO for the matrix 
element, pythia 8 for the parton shower, and DELPHES[Fav+i4] for 
a parameterized detector-level simulation. 

As seen in Eig. 3.14, variations of tan /3 does not lead to any 
kinematic difference and the production cross section simply scales 
as a function of tan /I. Hence we recommend to fix tan ^ to unity in 
the signal generation. 


MET 



Figure 3.14: Kinematic distributions 
of the signal process varying tan fti, in 
the case of a Higgs boson decaying 
into two b quarks, after parameterized 
detector simulation: no kinematic 
dependence is observed 


I Delta Phi - Two Leading Jet 

J 2 0.2 F 

C 

tt) 

iSO- 18 : 

0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 


• Mz=1000GeV, M ,=400GeV,tan|3=0.3 

• Mz-1000GeV, M ,=400GeV,tan|3=1 
Mz=1000GeV, M ,=400GeV,tan|3=3 

• Mz-1000GeV, M ,=400GeV,tan|3=5 

■ Mz=1000GeV, M ,=400GeV,tan|3=10 




1.5 2 2.5 3 

Delta Phi Leading jets (radians) 


(b) hip distance between the two b— jets 

Similarly, variations of gz do not lead to any kinematic changes. 
The value of gz for a given M^/ and tan ^ can be set according to 
the maximum value allowed by electroweak global fits and dijet 
constraints, as described in [BLWi4b]. Since this parameter does not 
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influence the kinematics, we leave it up to individual analyses on 
whether they generate benchmark points only according to these 
external constraints. 

Since the DM pair are produced as a result of the decay of A^, 
there are minimal kinematic changes when varying ni^ as long 
as < Mj^o /2 so that production is on-shell, as shown in 
Fig. 3.15 and 3.16 (before detector simulation). 


MET 



Figure 3 . 15 : Kinematic distributions 
of the signal process varying 
minimal kinematic dependency on 
as expected when A® is produced 
on-shell. Plots shown for = 

1000 GeV, = 300 GeV. 
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(b) hip distance between the two b— jets 

We recommend to produce signal events for a fixed gz = 0.8, 
tan j 6 = 1 and = 100 GeV. For these values, we scan the 2 -D pa¬ 
rameter space of M 2 ', M^o with M 2 ' = 600,800,1000,1200,1400 GeV, 
and M^o = 300,400,500,600,700,800 GeV with M^o < M 2 ' — 
for a total of 24 points. The choice of scan is justified by the sensi¬ 
tivity study in [BLWi 4 b]: the expected LHG sensitivity for Run -2 
is up to M 2 ' ~ 1.5 TeV. For the parameter scan, the DM mass is 
fixed to 100 GeV. For two M 2 ', M^o value sets, we vary the DM 
mass to obtain sample cross section for rescaling results. All LO 
cross sections for the various parameter scan points are reported in 
Appendix A. The parameter scan excludes the off-shell region, as 
the cross-sections are suppressed and the LHG would not have any 
sensitivity to these benchmark points in early data. 

The kinematic distributions with varying M 2 ' for fixed M^o 
are shown in Fig. 3 . 17 , while the dependency on M^o is shown in 
Fig. 3 . 18 . 
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Figure 3.16: Missing transverse mo¬ 
mentum distributions at generator 
level in the Z' +2HDM scenario for 
different values of the dark matter 
mass m^, with = 1 TeV and m^o = 
300 GeV (left) and = 1 TeV (right). 
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Figure 3.17: Kinematic distributions 
of the signal process varying M^/, for 
= 100 GeV, = 300 GeV. 
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Figure 3 . 18 : Kinematic distributions 
of the signal process varying M^o, for 
= 100 GeV, Mz' = 1000 GeV. 
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This model also allows for an additional source of Higgs plus 
£7 signal with a similar kinematics (Fig. 3 . 19 , shown with detector 
simulation samples) to the signal process from the decay of 7J 
hZ, where the Z decays invisibly The partial decay width for the Z' 
is: 


Tz'^/iZ = fecosasin^)-^^ 


\P\ I \P? 


Ml, 


+ 3 ^' 

M|, 


(3.12) 


The values for the Z' masses scanned for those samples should 
follow those of the previous samples, namely values of Mz' = 
600,800,1000,1200,1400 GeV. This signal process has no Ma de¬ 
pendence. 


3.2 EFT models with direct DM-boson couplings 

The EFT operators considered in this section do not have an imple¬ 
mentation of a simplified model completion for Dirac fermion Dark 
Matter available to date. They provide kinematic distributions that 
are unique to mono-boson signatures, and that in most cases are 
not reproduced by an equivalent simplified model.^ 

A complete list of effective operators with direct DM/boson cou¬ 
plings for Dirac DM, up to dimension 7 , can be found in [Cot+ 13 ; 
Car+ 13 ; CHH 15 ]. Higher dimensional operators, up to dimen¬ 
sion 8 , leading to Higgs+^j signatures, are mentioned in [Car+ 13 ; 
BLWiqb]. The first part of this Section outlines the main character¬ 
istics for a limited number of these models that could be considered 
in early Run -2 searches. However, the EFT approximation made for 
these operators can be problematic, see Ref. [BLWiqb] for discus¬ 
sion. For this reason, model-independent results as in Appendix B 
should be privileged over considering these operators as realistic 
benchmarks. 

However, the Forum discussion highlighted that the EFT ap¬ 
proach allows more model-independence when reinterpreting 
results, and that it is worth still considering interpretation of the 
results available in terms of these operators. Furthermore, once 
simplified models are available for those operators, EFT results can 
be used as a limiting case for consistency checks. We devote the 
end of this Section to a discussion on the presentation of results 
from this model, including an assessment of their reliability using a 
conservative procedure that is only dependent on EFT parameters. 

The studies in this Section have been performed using a UFO 
model within MadGrafh 5 _aMC@NLO V 2 . 2 . 3 , interfaced to 
PYTHiA 8 for the parton shower. The implementation of these mod¬ 
els is discussed further in Section 4 . 2 . 2 . 

3 . 2.1 Dimension 5 operators 

The lowest dimension benchmark operators we consider are effec¬ 
tive dimension 5 , such as the one depicted in Figure 3 . 20 . 


’ Wherever this is the case, for practical 
reasons one can only generation a 
simplified model result in the limiting 
EFT case, as the results can be rescaled 
and reinterpreted. 
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Figure 3 . 19 : Kinematic distributions 
of Z' ^ h exclusive production, 
Z' —> Zh exclusive production and 
Z' inclusive production for M^i = 
1000 GeV and M^o = 300 GeV 
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Figure 3 . 20 : Diagram for EFT operators 
giving rise to a Higgs+^j signature. 
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Following the notation of [Car+13], models from fhis cafegory 
have a Lagrangian fhaf, affer elecfroweak symmefry breaking, in¬ 
cludes ferms such as: 



(3-13) 


where mz and are fhe masses of fhe Z and W boson, and 
Zt‘ are fhe fields of fhe gauge bosons, x denofes fhe Dark Maffer 
fields and A5 is fhe effective field fheory scale. Nofe fhaf fhese 
operafors are of frue dimension 7, buf reduce fo effective dimension 
5 once fhe Higgs vacuum expecfafion values, confained in fhe W 
and Z mass ferms, are inserfed. As such, one expecfs fhaf fhese 
operafors would nafurally arise in UV complefe models where 
Dark Maffer inferacfs via a Higgs portal where heavy mediators 
couple to the Higgs or other fields in an exfended Higgs secfor. In 
such models fhe full fheory may be expecfed fo confain additional 
operators with Higgs-Dark Matter couplings [Djo+13]. The above 
operator also induces signatures with £7 in conjunction with Z 
and W bosons at tree level, as shown in Fig. 3.1, while at loop level 
it induces couplings to photon pairs and Zj through W loops. 

In these models, a clear relation exists between final sfafes wifh 
phofons, EW bosons and Higgs boson. 

As shown in Fig. 3.21, fhe kinematics of fhis model can be ap- 
proximafed by fhaf of a simplified model including a high-mass 
scalar mediafor exchanged in fhe s-channel described in Sec- 
fion 2.2.2. For fhis reason, fhe lisf of benchmark models wifh direcf 
boson-DM couplings for photon, Z and W only includes dimension 
7 operafors: if fhe scalar model wifh initial sfafe radiation of an EW 
boson is already generated, fhen ifs resulfs can be rescaled. 

The Higgs+^7 analysis, however, will nof consider fhe scalar 
simplified model as benchmark, due fo fhe very tow sensifivify in 
early LHC analyses, and will instead use fhis dimension 5 operafor. 
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Figure 3.21: Comparison of the miss¬ 
ing transverse momentum for the 
simplified model where a scalar medi¬ 
ator is exchanged in the s-charmel and 
the model including a dimension -5 
scalar contact operator, in the lep- 
tonic Z+Bt final state. All figures 
in this Section have been performed 
using a UFO model within Mad- 
Graph5_aMC@NLO V2.2.3, in¬ 
terfaced to PYTHiA 8 for the parton 
shower. 


3 . 2 . 1.1 Parameter scan 

The two parameters of this model are the scale of new physics A 
and fhe DM parficle mass. SM-DM coupling and new physics scale 
are relafed by = (246 GeV)/A. 

The inifial value of fhe new physics scale A chosen for fhe sam¬ 
ple generation is 3 TeV. This is a convention and does not affect the 
signal kinematics: the cross-section of fhe samples can be rescaled 
when deriving fhe consfrainfs on fhis scale. However, more care 
should be given when rescaling Higgs-F^j operafors of higher di¬ 
mensions, as differenf diagrams have a differenf A dependence. 

The DM mass values for fhe benchmark poinfs fo be simulafed 
are chosen fo span a sufficienf range leading fo differenf kinemafics, 
thaf is wifhin fhe LHC sensifivify for early searches and fhat is con- 
sisfenf across fhe various signafures and EFT operafors. We fhere- 
fore sfarf fhe mass scan af =1 GeV, where collider experimenfs 
are complemenfary fo direcf and indirecf defecfion and choose fhe 
lasf poinf corresponding fo a DM mass of 1 TeV. We recommend a 
scan in seven mass poinfs, namely: 

= 1,10,50,100,200,400,800,1300 GeV. 

A sef of kinemafic disfribufions from fhe Higgs-F^p signafure 
where fhe Higgs decays info fwo b—quarks is shown in Fig. 3 . 22 , 
for poinfs similar fo fhose of fhe grid scan proposed. 


3.2.2 Dimension 7 operators 

The dimension -7 benchmark models confain fhe SU{2)i x Lr(l)y 
gauge-invarianf couplings befween DM fields and fhe kinefic ferms 
of fhe EW bosons. The CP-conserving scalar couplings of fhis type 
can be written as 


G_ 






( 3 - 14 ) 
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Figure 3,22: Comparison of the kine¬ 
matic distributions for the two leading 
b— jets (from the Higgs decay) in the 
model with direct interactions between 
the Higgs boson and the DM particle, 
when varying the DM mass. 




(c) Angular distance between the two leading fo—jets 
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are the iJ(l)Y SU{2)i field strength tensor, respectively, and g2 


denotes the weak coupling constant. In the case of the pseudoscalar 
couplings, one has instead 


^ X75X -h ^ X75X W', (3.15) 


where Bpy = 1/2e^yXpB^P and = l/2epy\pW^'^P are the dual 
field strength tensors. In addition to the CP-conserving interactions 
(3.14) and (3.15), there are also four CP-violating couplings that are 
obtained from the above operators by the replacement XX ^ X75X- 

The effective interactions introduced in (3.14) and (3.15) appear 
in models of Rayleigh DM [WY12]. Ultraviolet completions where 
the operators are generated through loops of states charged under 
Lr(l)y and/or SU{2)i have been proposed in [WY13] and their 
LHC signatures have been studied in [Liu+13]. If these new charged 
particles are light, the high-px gauge bosons that participate in 
the Bt processes considered here are able to resolve the substruc¬ 
ture of the loops. This generically suppresses the cross sections 
compared to the EFT predictions [HKU13], and thus will weaken 
the bounds on the interaction strengths of DM and the EW gauge 
bosons to some extent. Furthermore, the light charged mediators 
may be produced on-shell in pp collisions, rendering direct LHC 
searches potentially more restrictive than Bj searches. Making the 
above statements precise would require further studies beyond the 
timescale of this forum. 

Since for A5 = Ap the effective interactions (3.14) and (3.15) 
predict essentially the same value of the mono-photon, mono-Z and 
mono-W cross section [Car+13; CHH15], we consider below only 
the former couplings. We emphasize however that measurements 
of the jet-jet azimuthal angle difference in £7 +2] events may be 
used to disentangle whether DM couples more strongly to the 
combination BpyBP'' {W^yW’P'") or the product BpyW'^ {W^yW’P'^) of 
field strength tensors [Cot+13; CHH15]. 

After EW symmetry breaking the interactions (3.14) induce direct 
couplings between pairs of DM particles and gauge bosons. The 
corresponding Feynman rule reads: 



(3.16) 


where p, (p,) denotes the momentum (Lorentz index) of the vector 
field Vi and for simplicity the spinors associated with the DM fields 
have been dropped. The couplings take the form: 


SyZ ~ ^luCw (ci C2) , 

^ZZ 'X f 


(3-17) 


S'ww — C2 , 
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with Su! {Cw) the sine (cosine) of the weak mixing angle. Note that 
our coefficients Ci and C2 are identical to the coefficients Cg and 
Cw used in [CHH15], while they are related via ki = Cw^Ci and 
^2 = s-w^C2 to the coefficients k\ and k2 introduced in [Car+13]. 

The coefficients ci and C2 appearing in (3.17) determine the rel¬ 
ative importance of each of the £7 charmels and their correlations. 
For example, one observes that: 

• Only C2 enters the coupling between DM and W bosons, mean¬ 
ing that only models with C2 7^ 0 predict a mono-W signal; 

• If Cl = C2 the mono-photon (mono-Z) signal does not receive 
contributions from diagrams involving Z (photon) exchange; 

• Since numerically c^/s^ ~ 3.3 the mono-photon charmel is 
particularly sensitive to Ci. 


3.2.2.1 Parameter scan 

As stated above and shown in Ref. [Nel+14], the kinematic distribu¬ 
tions for dimension-7 scalar and pseudoscalar operators only shows 
small differences. This has been verified from a generator-level 
study: the signal acceptance after a simplified analysis selection (£7 
>350 GeV, leading jet pj > 40 GeV, minimum azimuthal difference 
between either of the two jets and the £7 direction > 0.4) is roughly 
70% for both models, independent from the coefficients ci and C2. 
We therefore only suggest to generate one of the two models. 

The differences in kinematics for the various signatures are neg¬ 
ligible when changing the coefficients Ci and C2, since these coeffi¬ 
cient factorize in the matrix element. Only the case Ci = C2 = 1 is 
generated as benchmark; other cases are left for reinterpretation as 
they will only need a rescaling of the cross-sections. 



Figure 3.23: distribution for the 

dimension-7 model with a hadroni- 
cally decaying Z in the final state, for 
the scalar and pseudoscalar opera¬ 
tors representing direct interactions 
between DM and bosons. The values 
of the coefficients in the legend are 
multiplied by 100. 
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3.2.3 Higher dimensional operators 

Many higher dimensional operators can induce signals of photons 
or YJ 12.1 H bosons in the final state. A complete list can be found 
in Refs. [Car+14; BLWi4b; PS14] and references therein. 

Although with lower priority with respect to the operators 
above, a representative dimension-8 operators can be chosen as 
benchmark, with the form: 

In this case, the new physics scale is A is related to the coupling 
1 

of the DM as i/;^- = .An advantage of this operator is that it 
includes all signatures with EW bosons, allowing to assess the rel¬ 
ative sensitivity of the various charmels with the same model. The 
kinematics for this operator is different with respect to other oper¬ 
ators, leading to a harder £ 7 spectrum, as illustrated by comparing 
the leading b—jet distribution for the dimension 5 operator to the 
dimension 8 operator. 


3.2.4 Validity ofEW contact operators and possible completions 

It is important to remember that the operators described in this 
section may present problems in terms of the validity of the contact 
interaction approach for the energy scales reached at the LHC. 

As outlined in [BLWi4b], designing very high £7 search signal 
regions that are exclusively motivated by the hard £ 7 spectra of 
the dimension 7 and 8 operators will mean that the momentum 
transfer in the selected events is larger. This in turn means that 
processes at that energy scale (mediators, particles exchanged in 
loops) are accessible, and a simple contact interaction will not be 
able to correctly describe the kinematics of these signals. 

Contact interaction operators like the ones in this section remain 
useful tools for comparison of the sensitivity of different search 
charmels, and for reinterpretation of other models under the cor¬ 
rect assumptions. To date, while UV-complete models are known, 
their phenomenology has not been studied in full detail as their 
completion involves loops 3 . 

However, this may be the focus of future theoretical exploration, 
as discussed in Ref. [CHH15]. An example of a complete model for 
scalar DM corresponding to the dimension-5 operator is provided 
in the Appendix A. Providing results for the pure EFT limit of these 
models will prove useful to cross-check the implementation of 
future. 

Given these considerations, we recommend to present results for 
these models as follows: 

• Deliver fiducial limits on the cross section of any new physics 
events, without any model assumption, according to the guide¬ 
lines in Appendix B. 


3 An example case for the need of loop 
completions is a simplified model with 
an additional scalar exchanged at tree 
level. The scalar couples to WW and 
ZZ in a gauge-invariant way. Integrat¬ 
ing out the mediator does not lead to 
the Lorentz structure of a dimension-7 
operator, so it is not possible to gener¬ 
ate dimension-7 operators that satisfy 
gauge and Lorentz invariance at the 
same time. A model with a spin -1 
mediator cannot be considered as an 
candidate for completion either, since 
dimension-7 operators only have scalar 
or pseudoscalar couplings. 


Events / 30 GeV Events / 30 GeV 
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Figure 3.24: Comparison of the trans¬ 
verse momentum for the leading b— jet 
from the Higgs decay for a dimension 
5 and dimension 7 operator with direct 
boson-DM couplings. 



(b) Dimension 7 operator 
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Assess the percentage of events that pass a condition of validify 
for fhe EFT approximafion fhat does nof depend on a specific 
completion, and presenf resulfs removing of fhe invalid evenfs 
using fhe procedure in Section 5 alongside fhe raw EFT resulfs. 



4 

Implementation of Models 


4.1 Implementation of s-channel and t-channel models for tj 
+X analyses 

In the studies to date, a number of different Monte Carlo tools have 
been used to simulate DM signals. In this Chapter, we make rec¬ 
ommendations on the accuracy at which simulations should be 
performed for differenf final sfafes. We also provide explicif exam¬ 
ples of codes and implemenfafions (including specific seffings) fhaf 
have been used to obtain the results in this report. We stress that 
these recommendations are based on the current status of publicly 
available codes and users should always check whefher new resulfs 
af a better accuracy have appeared in fhe meantime. In fhaf case, 
we recommend fo updafe fhe corresponding analyses direcfly us¬ 
ing the new releases and/or codes, and in case this would not be 
possible, to at least take into account the new information in the 
analysis (e.g., via a MC comparison with the latest predictions, or 
by effectively using global/local fC-f actors). For all models included 
in this report, pythia 8 has been used to provide the parton shower 
simulation. Nevertheless, we note that showering matrix element 
events with Herwig [Bah+08; Cor+02; Cor+01; Mar+92] should be 
considered as an equally valid alternative. 

4.1.1 Implementation of s-channel models for mono-jet signature 

These models include those discussed in Secs. 2.1 and 2.2. In mono¬ 
jet analyses, i.e. when final sfafes are selecfed wifh a few jefs and 
observables and in particular fhe Bj specfrum depend upon 
fhe accuracy of fhe simulation of QCD radiation. For fhe vecfor 
and axial vecfor models, fhe currenf sfafe of fhe arf is NLO+PS. If 
is parficularly simple fo obfain simulations for fhese processes af 
NLO+PS and even for merged samples af NLO accuracy, sfarfing 
from SM implemenfafions. We fherefore recommend simulations fo 
be performed af NLO+PS, and in case mulfi-jef observables are em¬ 
ployed, by merging samples wifh differenf mulfiplicifies. Resulfs af 
such accuracy can be obfained eifher in dedicafed implemenfafions, 
such as fhaf of fowheg [HKR13], or via general purpose NLO fools 
like MadGrafh 5 _aMC@NLO employing available UFO models 
af NLO. A tesfing version of fhe full sef of fhese UFO models has 
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been made available only in June 2015 [New], For this reason, it 
was not used as part of the studies of fhis Forum on inifial Run-2 
benchmark models. Nevertheless, we encourage further study of 
fhese UFO models by fhe ATLAS and CMS collaborafions. 

A sfudy using POWHEG [FIKR13; FW13] has shown fhaf fhe 
NLO correcfions resulf in a subsfanfial reduction in fhe dependence 
on fhe choice of fhe renormalizafion and facforizafion scales and 
hence a reduced fheorefical uncerfainfy on fhe signal prediction. 

For fhe cenfral choice of renormalizafion and facforizafion scales, 
fhe NLO correcfions also provide a minor enhancemenf in fhe cross 
secfion due fo fhe jef veto fhaf has been so far employed in Run-i 
analyses. 

For the scalar and pseudoscalar models, the lowest order process 
already involves a one-loop amplitude in QCD. Because of fhe com¬ 
plexify of performing NLO calculafions for fhis class of processes 
and in parficular fhe absence of general mefhods for compufing 
fwo-loop virfual contribufions, only LO predictions are currenfly 
available. These can be inferfaced fo shower programs exacfly as 
usual free-level Born compulations, i.e. by considering one parfon 
mulfiplicify af fhe fime or by merging different parton multiplic¬ 
ities via CKKW or MLM schemes to generate inclusive samples 
with jet rates at LO accuracy. For spin-o mediators in the mono-jet 
final sfafe, fhe fop-quark loop is fhe mosf imporfanf considerafion. 
The mafrix elemenf implemenfafion wifh exacf fop-loop depen¬ 
dence of fhe s-channel spin-o mediafed DM producfion is available 
in MCFM [FW13; FIar+15] ^ af fixed order and in powheg [FIR15] 
and MadGraph5_aMC@NLO [New] for even! generation af 
LO+PS level. The powheg and mcfm implemenfafions include fhe 
finite fop quark mass dependence for DM pair producfion and 
one exfra parfon af LO. The same processes are available in Mad- 
Graph5_aMC@NLO V2.3 and could be made available in fhe fu- 
fure in codes like Sherpa+OpenLoops/GoSam, including up fo fwo 
extra partons in the final sfafe. Samples can be merged employing 
CKKW, Kt-MLM procedures. 

Mosf of fhe resulfs fhat have been presented in this document 
for fhese processes have been obfained wifh powheg inferfaced fo 
PYTHiA 8, mafching fhe sfafe of fhe arf calculation as of Spring 2015. 
For fufure reference, we document the specific settings needed fo 
run fhe powheg generafion for fhe Dark Matter models so fhey can 
serve as nominal benchmarks for fhe early Run-2 ATLAS and CMS 
DM analyses, powheg paramefer cards for all models can be found 
on fhe Forum SVN repository [Fori; Foro; Forn; Form]. 


4.1.1.1 POWHEG configuration for s-channel DM models 

The lafesf powheg release is available for download using fhe in- 
sfrucfions af http://powhegbox.mib.infn.it/. The Forum recom¬ 
mends using at least version 3059 . 


’ Only the scalar mediator is available 
in the public release. 
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powHEG can generate either unweighted (uniformly-weighted) 
or weighted events. The relevant ke5rwords in the input card are 
bornsuppfact and bornktmin. 

1. unweighted events: 
bornsuppfact: negative or absent 
bornktmin PT 

This runs the program in the most straightforward way, buf if 
is likely nof fhe more convenienf choice, as will be explained 
below. POWHEG will generafe unweighfed evenfs using a 
sharp lower cuf (wifh value PT) on fhe leading-jef pj. Since 
fhis is a generafion cuf, fhe user musf check fhaf fhe choice 
of bornktmin does not change the cross section for signal 
events passing analysis selections. It is good practice to use 
as a value in the input card a transverse momentum 10-20% 
smaller than the final analysis selection on £7, and check that 
the final result is independent, by exploring an even smaller 
value of bornktmin. The drawback of using this mode is that 
it is difficult to populate well, and in a single run, both the 
low-pT region as well as the high-px tail. 

2. weighted events: 
bornsuppfact PTS 
bornktmin PT 

POWHEG will now produce weighted events, thereby allowing 
to generate a single sample that provides sufficient statistics 
in all signal regions. Events are still generated with a sharp 
lower cut set by bornktmin, but the bornsuppfact parameter is 
used to set the event suppression factor according to 

^(^t) = , 2 , . • (4-1) 

kj + bornsuppfact^ 

In this way, the events at, for instance, low £7, are suppressed 
but receive higher weight, which ensures at the same time 
higher statistics at high £7. We recommend to set bornsuppfact 
to 1000. 

The bornktmin parameter can be used in conjunction with 
bornsuppfact to suppress the low £7 region even further. It 
is recommended to set bornktmin to one-half the value of the 
lowest £7 selection. For instance, for the event selection used 
in the CMS/ATLAS monojet analyses, assuming the lowest 
£7 region being defined above 300 GeV, the proposed value 
for bornktmin is 150. However, this parameter should be set 
keeping in mind the event selection of all the analyses that 
will use these signal samples, and hence a threshold lower 
than 150 may be required. 

The POWHEG monojet implementations can generate events using 
two expressions for the mediator propagators. The default setup 
(i.e if the keyword runningwidth is absent, commented out or set 
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to o) is such that a normal Breit-Wigner function is used for fhe 
propagafor: in fhis case, fhe expression 

Q^-M^ + iMT 


is used for fhe propagator's denominafor, where Q is fhe virfual- 
ify of fhe mediator, and M and F are ifs mass and widfh, respec- 
fively. This is fhe more sfraighfforward, simple and fransparenf 
opfion, and if was used for fhe Forum sfudies. If should be fhe 
mefhod of choice, unless one approaches regions of paramefer 
space where gamma/M sfarfs to approach order i values. In 
fhose cases, a more accurafe modelling (or af leasf a check of fhe 
validify of fhe fixed widfh approach) can be achieved by using a 
rurming widfh: by seffing fhe runningwidth token to i, powheg 
uses as fhe denominafor of fhe mediaforaAZs propagafor fhe 
expression 


q" 


+ i 


r 

M' 


which is known fo give a more realistic description. See Ref. [Bar+89] 
for a discussion. 


• Sef fhe parameters defining fhe bounds on fhe invarianf mass of 
fhe Dark Maffer pair, mass_low and mass_high, fo -1. In fhis way, 
POWHEG will assign values infernally. 

• The minimal values for ncalll, itmxl, ncall2, itmx2 are 250000, 
5,1000000, 5 for fhe vector model, respecfively. 

• The minimal values for ncalll, itmxl, ncall2, itmx2 are 100000, 
5,100000, 5 for fhe scalar fop-loop model, respecfively. 

• When NLO correcfions are included (as for insfance in fhe vecfor 
model), negafive-weighfed evenfs could happen and should be 
kepf in fhe evenf sample, hence withnegweights should be sef fo 
1. If needed, fheir fraction can be decreased by seffing foldsci 
and f oldy fo bigger value (2 for insfance). fold phi can be kepf fo 
1. 


• One should use fhe aufomafic calculation of systematic uncer- 
fainfies associated wifh fhe choice of hard scale and PDFs as 
described in Section 6. 

• idDM is fhe infeger fhaf identifies fhe DM parficle in fhe Monfe 
Carlo evenf record. This should be chosen so fhaf ofher fools can 
process fhe powheg oufpuf properly. 

POWHEG in ifself is nof an evenf generafor and musf be inferfaced 
wifh a fool fhaf provides parfon showering, hadronizafion, etc. For 
some fime, a pythia 8 [Sjo+15] interface has existed for powheg. 

The PYTHIA 8 runtime configurafion is fhe following: 

POWHEG:veto = 1 
POWHEG:pTdef = 1 
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POWHEG:emitted = 0 
POWHEG:pTemt = 0 
P0WHEG:pThard = 0 
POWHEG:vetoCount = 100 
SpaceShower:pTmaxMatch = 2 
TimeShower:pTmaxMatch = 2 

As always, it is recommended to use the latest pythia 8 release, 
available at http://home.thep.lu.se/~torbjorn/Pythia.html. At 
the time of this report, the latest version is 8. 209. 

4 . 1.2 Merging samples with different parton multiplicities 

For the models discussed in the previous section, it is important 
to calculate the hard process as accurately as possible in QCD. For 
many other signal models, the £ j signature depends more upon the 
production and decay of fhe mediator. In some cases, observables 
builf in ferms of fhe jefs presenf in fhe final sfafe are considered, 
somefhing thaf assumes inclusive samples accurate in higher jet 
multiplicities are available. In these cases, one can employ LO+PS 
simulations where different parton multiplicities are merged and 
then matched to parton shower, using schemes such as CKKW or 
MLM merging. 

Here, we consider the example of an EFT model produced in 
associafion wifh up fo 2 addifional QCD par tons. A Monfe Carlo 
sample based on fhis mefhod could be used in alfernafive fo a 
NLO+PS sample for describing shapes and jef disfribufions (buf 
nof for fhe overall normalisafion which would sfill be af LO). The 
mefhodology described here could also be used for fhe f-charmel 
model discussed in Sec. 2.3. 

For fhe calculation of free-level merged samples for DM signals, 
fools fhaf can read UFO files and implement multi-parton merg¬ 
ing should be employed, such that MadGrafh 5 _aMC@NLO 
(+PYTHIA 8 or HERWIG++) and Sherfa [Hoe+15]. In this report 
we have mostly employed MadGrafh 5 _aMG@NLO. Mad- 
Grafh 5 _aMC@NLO provides a flexible and easy-fo-use frame¬ 
work for implemenfing new models via fhe EeynRules package. 
MadGrafh 5 _aMC@NLO can perform bofh LO and NLO calcu- 
lafions in QCD, mafched/merged fo parfon showers [AVM09]. Eor 
NLO ones, dedicated UEO model implemenfafions af NLO should 
be used. Several UEO models af NLO are publicly available fhaf 
while nof developed specifically for DM, are suifable fo make mode 
independenf simulafions af NLO accuracy, including mulfipar- 
fon merging via fhe ExEx technique [EE12]. A dedicafed DM UEO 
implemenfafion has been developed and if has been released as a 
fesfing version [New]. 

Merging evenfs generafed via mafrix elemenfs wifh differenf 
number of parfons in fhe final sfafe can be achieved by a judicious 
procedure fhaf avoids double counting of fhe parfons from mafrix 
elemenfs and parfon showering. Several merging fechniques are 
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available. Based on some comparative studies [Alw+o8], there is 
some advantage to using the CKKW-L merging scheme [LP12] 
implemented in pythia 8. Alternatively, one can use the fcj-MLM 
scheme also available in pythia 8. 

4.1.2.1 Generation of the LHE file 

The example presented here is a D5 EFT model, and includes 
tree-level diagrams with partons. We stress that Mad- 

Graph 5 _aMC@NLO, like powheg, is not in itself and event gen¬ 
erator, but must be interfaced wifh an event generator through an 
LHE file. The producfion of fhe LHE file proceeds fhrough setting 
fhe process paramefers and fhe run paramefers. 

The process paramefers are: 

import model MODELNAME 

generate p p > chi chi- [QCD] @0 

add process p p > chi chi- j [QCD] @1 

add process p p > chi chi- j j [QCD] @2 

The runtime parameters are more numerous, and define the 
collider properties, PDF sets, etc. The specific parameters needed 
for matching are, for the example of CKKW-L matching: 


ickkw = 0 

ktdurham = matching scale 

dparameter = 0.4 

dokt = T 

ptj=20 

drjj=0 

mmj j=0 

ptjlmin=0 


For different kinds of matching, a different choice of ickkw and 
related parameters would be made. 


4.1.2.2 Implementation of the CKKW-L merging 

To illustrate the settings related to merging different multipliticities, 
the EFT D 5 samples were generated with MadGraph 5 _aMC@NLO 
version 2.2.2 and showered in pythia 8.201, using the Madgraph 
parameters in the previous section (Sec. 4.1. 2.1). 

The PYTHIA 8 parameters for the CKKW-L kj-merging scheme 
are: 


Merging:ktType = 1 

Merging:TMS = matching scale 

1000022:all = chi chi- 2 0 0 30.0 0.0 0.0 0.0 0.0 

1000022:isVisible = false 

Merging:doKTMerging = on 

Merging:Process = pp>{chi,1000022}{chi-, -1000022} 

Merging:nJetMax = 2 
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The matching scales should be the same for the generation and 
parton showering. In the model implementation, the particle data 
group ID 1000022 is used for weakly inferacfing dark maffer can- 
didafes. Since fhis is a Majorana parficle by defaulf (wifh no cor¬ 
responding anfi-parficle), and fhe model produces a DM Dirac 
fermion, fhe parficle properties are changed accordingly. Also, 
fhe DM mass is sef fo 30 GeV. The Merging: Process command 
specifies fhe lowesf parfon emission process generafed in Mad- 
Graph 5 _aMC@NLO and Merging: nJetMax = 2 gives fhe max¬ 
imum number of addifional parfon emissions wifh respecf fo fhe 
lowesf parfon emission process. 

In general, if is desired fo fake fhe hard parfon emissions from 
the matrix element generation in MadGraph 5 _aMG@NLO and 
allow PYTHiA 8 to take care of soff emissions only. The fransifion 
befween fhese two regimes is defined by fhe mafching scale and ifs 
optimal value can be defermined by sfudying fhe cross-secfion as a 
function of fhe number of jefs (differenfial jef rafes). The differential 
rafes —y give fhe number of evenfs which pass from i jefs fo 

j jefs as fhe kj value increases beyond kcut- An optimal mafching 
scale should lead to smooth differential jet rates. 

Two examples of differential jet rates, using matching scale 
30 GeV and 80 GeV, from fhe EFT D5 sample generafed as de¬ 
scribed in fhe previous section are given in Fig. 4.1 and 4.2, respec¬ 
tively. Although a kink is visible around the matching scale value 
in both cases, the 80 GeV scale leads to smoother distributions. 

In order to find fhe opfimal mafching scale, additional samples 
with matching scale 50, 70, and 90 GeV are generated as well and a 
detailed comparison of fhe differenfial jef rafes close fo fhe fransi¬ 
fion region is shown in Fig. 4.3. The largesf differences among fhe 
samples are visible for fhe 1 —>• 2 jefs fransifion where fhe 30 GeV 
and 50 GeV scale lead fo a drop of fhe rafes around fhe mafching 
scale values. On fhe confrary fhere is a hinf of an increased rafe 
around fhe mafching scale value in fhe sample generafed wifh fhe 
90 GeV scale. Therefore, we recommend fo use 80 GeV as fhe base¬ 
line mafching scale. 

The prescription for fhe evenf generafion given in Secfion4.i.2.2 
sfarfs wifh fhe emission of o parfons and ends wifh maxim 2 par- 
fons in addition. Producing the samples separately allows to inves¬ 
tigate the relative composition of fhe individual samples in various 
parfs of fhe phase space. Figure 4.4 shows fhe £7 disfribufion of 
the EFT D5 sample with the matching scale at 80 GeV. The plot 
reveals that the o-parton sample gives the dominant contribution 
in the region below the matching scale value that rapidly decreases 
at higher £7. Assuming the lowest analysis £7 cut in early Run-2 
mono-jet analyses at 300 GeV, the generation of fhe o-parfon emis¬ 
sion sample can be safely omitted as if only gives < 1% confribu- 
fion af £7 > 300 GeV. For fhe 1- and 2-parfon emission samples, 
one can use a generafor cuf on fhe leading parfon pj, pt j Imin, in 
order fo avoid generafing low £7 evenfs fhaf are irrelevanf for fhe 
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vertical line is drawn at the matching 
scale. 
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analysis. 
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In order to describe the signal kinematics correctly and save time 
during MC production, the parton emissions will only be generated 
up to a certain multiplicity. The higher multiplicity samples usually 
have small enough cross sections and the corresponding parts of 
the phase space can be sufficiently approximated by parton show¬ 
ering in PYTHiA 8. A dedicated study comparing samples generated 
with up to 1-, 2-, or 3-parton multiplicities was performed, using 
again fhe seffings for fhe CKKW-L k^-merging wifh fhe 80 GeV 
mafching scale and fhe Merging: nJetMax paramefer adjusfed ac¬ 
cordingly. Figure 4.5 shows fhe Bj disfribufion of fhe samples af 
Bt > 250 GeV. 

With an event selection requiring B t and the leading jet pj being 
larger than 250 GeV, the sample generated with up to 1 parton has 
10.3% larger yield compared to the sample with up to 3 partons, 
while the yield of fhe sample wifh up to 2 partons is only 2.3% 
larger. If an addifional cuf is applied allowing for up fo 3 jefs wifh 
Pj > 30 GeV, the agreement improves to 3.2% larger for up fo 
1 parfon and 0.7% larger for up fo 2 parfons, compared wifh up 
fo 3 parfons. A similar comparison is shown in Fig. 4.6 for fhe jef 
mulfiplicify in fhe evenfs wifh fhe leadning jef pj > 250 GeV, where 
an agreemenf af fhe level of ^ 3 % between fhe samples wifh up 
fo 2 and 3 parfon emissions is observed for number of jefs up fo 7. 
This jusfifies if is sufficienf fo produce samples wifh up fo 2 parfon 
emissions only af fhe generafor level and ignore generating higher 
parfon emissions. 

4.1.3 Implementation oft-channel models for the jet+Bj final state 

The simulafions for f-channel models are available via LO UFO 
implemenfafions, where evenfs are generafed af LO+PS accuracy. 
The UFO file and paramefer cards for fhe f-charmel models wifh 
couplings fo lighf quarks only [PVZ14] can be found on the Forum 
SVN repository [Forj]. The model files from Ref. [Bel+12] can also 


Figure 4.4: Missing transverse mo¬ 
mentum distributions for EFT D5 
sample with CKKW-L matching scale 
at 80 GeV. Individual contributions 
from the o-, 1- and 2-parton emission 
samples are shown. 
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Figure 4.5: Missing transverse mo¬ 
mentum distributions for EFT D5 
sample with CKKW-L matching scale 
at 80 GeV produced with maximum 
1 (black), 2 (red) and 3 (blue) partons 
emitted at the generator level. The 
ratios are shown with respect to the 
latter sample. 


(a) No jet multiplicity cut 
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Figure 4.6: Multiplicity of jets with 
Pt > 30 GeV and |j;| <2.8 for EFT D5 
sample with CKKW-L matching scale 
at 80 GeV produced with maximum 
1 (black), 2 (red) and 3 (blue) partons 
emitted at the generator level. The 
ratios are shown with respect to the 
latter sample. The leading jet pj is 
required to be larger than 250 GeV. 


be found on the repository [Fori]. The latter is the implementation 
that has been used for the studies in this report: in the monojet case 
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there are only cross section differences befween fhis model and fhe 
model in [For]]. 

Mulfi-parfon simulation and merging are necessary and require 
particular care for fhis model: fhis has nof been a topic of defailed 
sfudies wifhin fhe Forum, and we suggesf fo follow fhe procedure 
ouflined in Ref. [PVZ14]. 

4.1.4 Implementation ofs-channel and t-channel models with EW bosons 
in the final state 

Currenfly, simulafions for mosf of fhese models are available via 
LO UFO implemenfafions, allowing even! generafion af fhe LO+PS 
accuracy. We note, however, fhaf inclusion of NLO corrections 
would be possible. In MadGraph5_aMC@NLO, for example, 
fhis amounfs fo simply upgrading fhe currenfly employed UFO 
models fo NLO, where fhe calculafions exisf for fhis class of pro¬ 
cesses. Flowever, fhis was nof available wifhin fhe timescale of fhe 
Forum fowards simulation of early Run-2 benchmarks. As a con¬ 
sequence, in fhis work we have used LO UFO implemenfafions 
wifhin MadGraph5_aMC@NLO 2.2.3 interfaced fo pythia 8 for 
fhe parfon shower. The corresponding paramefer cards used for fhe 
Run-2 benchmark models can be found on fhe Forum SVN repos- 
ifory [Fora]. This is fhe implemenfafion fhaf will be used for early 
Run-2 LHC Dark Matter searches. 

None of fhese models requires merging samples wifh differenf 
parfon multiplicities since fhe visible signal comes from fhe produc- 
fion of a heavy SM boson whose fransverse momenfum disfribufion 
is sufficienfly well described af LO+PS level. As a resulf, no special 
runtime configuration is needed for pythia 8. 

4.1.5 Implementation ofs-channel and t-channel models with heavy 
flavor quark signatures 

Dedicafed implemenfafions for DM signals in fhis final sfafe are 
available af LO+PS accuracy. However, fhe sfafe of fhe arf of fhe 
simulafions for tt and bb wifh a generic scalar and vecfor mediator 
is NLO+PS accuracy. For example, simulafions for tt + scalar can be 
obfained via powheg and sherpa sfarfing from fhe SM implemen¬ 
fafions. In MadGraph5_aMC@NLO, all final relevanf final sfafes, 
spin-0 (scalar and pseudo scalar) and spin-i, (vecfor and axial) are 
available af NLO+PS via fhe dedicafed NLO UFO for DM has been 
released in June 2015 [New]). 

In fhe work of fhis Forum, simulafions for fhe tt and bb sig- 
nafures of fhe scalar mediator model have been generafed sfarf¬ 
ing from a leading order UFO wifh MadGraph5_aMC@NLO 
2.2.2, using PYTHIA 8 for fhe parfon shower. The UFO file and 
paramefer cards fhaf will be used as benchmarks for early Run- 
2 searches in fhese final sfafes can be found on fhe Forum SVN 
repository [Ford]. Mulfi-parfon merging has been used for fhe bb 
case buf if has nof been sfudied in defail wifhin fhis Forum. The 
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b-flavored DM model of Section 2.3.3 simulated at LO+PS us¬ 
ing MadGrafh5_aMC@NLO V2.2.3 and pythia 8 for the parton 
shower. The corresponding UFO and parameter files can be found 
on the Forum SVN repository [Forg]. 


4.1.5.1 Quark flavor scheme and masses 

In the case of bb final state an additional care should be taken when 
choosing the flavor scheme generation and whether quarks should 
be treated as massive or massless. 

The production of DM+bb, Dark Matter in association with b jets 
via a decay of a (pseudo) scalar boson, is dominated in simplified 
mediator models by the gluon-gluon initiated production, similar 
to the production of Z+bb at the LHC. The Z+bb process has been 
studied in detail in the Z(ll)+b-jets final state, which can be used to 
validate both the modeling of DM+bb and, its main background, 
Z{vv)+bb. In this context, the pj of the Z boson is related to the 
observed MET, whereas the b-jet kinematics determines the ratio of 
mono-b/di-t> signatures in the detector. 

For basic kinematic criteria applied to Z+bb production, this 
process leads in ^ 90 % of the events to a signature with only 1 
b-jet in the acceptance ('Z+ib-jet production') and only in ~ 10 % 
of the events to a signature with 2 t>-jets in the detector {'Z+2b- 
jets production). The production cross section of the Z+bb process 
can be calculated in the 'five-flavor scheme', where b quarks are 
assumed massless, and the 'four-flavor scheme', where massive b 
quarks are used [Cam+04; MMW05; Cam+06]. Data slightly favour 
the cross-section predictions in the five-flavor scheme [CMSiqa] 
for the 1 b-jet signature. In this document we have preferred the 
5-flavor scheme due to its simplicity and cross sections and models 
in the 5-flavor scheme are available in the repository. The PDF used 
to calculate these cross section is NNPDF3.0 (lhaid 263000). 

On the other hand, both data [CMSiqa; CMS13; CMS15C] and 
theoretical studies [Fre+11; Wie+15] suggest that the best modelling 
of an inclusive Z+bb sample especially for what concerns fe-quark 
observables, is achieved at NLO+PS using a 4-flavor scheme and a 
massive treatment of the t>-quarks. In Figure 4.7 we show that, at 
LO, as expected, no appreciable difference is visible in the kinemat¬ 
ics between either flavor scheme used for UM+bb. In our generation 
we have used NNPDF3.0 set (lhaid 263400). 




Figure 4.7: Comparison of the jet mul¬ 
tiplicity (left) and angular correction 
AROi/h) (right) for the DM+bb scalar 
model generated in the 4-flavor and 
5-scheme. The samples are generated 
for nir = 1 GeV and = 10 GeV. 
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4.2 Implementation of specific models for V + tj analyses 

4.2.1 Model implementation for mono-Higgs models 

Currently, simulations for most of fhese models are available via 
LO UFO implemenfafions, allowing even! generafion af fhe LO+PS 
accuracy. We nofe, however, fhaf fhe inclusion of NLO correcfions 
would be possible buf nof available in time for fhe conclusion of 
fhese sfudies. In MadGraph5_aMC@NLO, for example, fhis 
amounfs fo simply upgrading fhe currently employed UFO mod¬ 
els to NLO. Simulation of loop-induced associafed production of 
DM and Higgs is also possible wifh fhe exacf fop-quark mass de¬ 
pendence. In MadGraph5_aMC@NLO, for example, fhis can be 
obfained from fhe NLO UFO SM and 2HDM implemenfafions. 

In fhis work all fhree Higgs+^p models have been generafed af 
leading order wifh MadGraph5_aMG@NLO 2.2.2, using pythia 8 
for fhe parfon shower. No merging procedure has been employed. 

The LO UFO implemenfafions of fhe scalar and vector models fhaf 
will be used as early Run-2 benchmarks can be found on fhe Forum 
SVN reposifory [Forh], while fhe 2HDM model can be found af fhis 
link [Forb]. 

As a final fechnical remark, we suggesf always fo lef fhe shower 
program handle fhe h decay (and fherefore fo generafe a sfable h af 
the matrix element level). In so doing a much faster generation is 
achieved and the h branching ratios are more accurately accounted 
for by fhe shower program. 

4 . 2 . 1.1 MadGraph5_aMG@NLO details for scalar mediator Higgs+MET 
model 

The case of fhe associafed producfion of a Higgs and scalar me¬ 
diator via a fop-quark loop can be eifher considered exacfly or 
via an effective Lagrangian where fhe fop-quark is infegrafed ouf. 

While fhis laffer model has been shown nof fo be reliable [HKU13; 
HLVV14; BG90], for simplicify we have chosen fo perform fhe sfudy 
in fhis free-level effective formulafion. A full sfudy of fhe process 
including finife fop-quark mass and parfon shower effecfs is possi¬ 
ble yef leff for fufure work. 

4 . 2 . 1.2 MadGraph5_aMG@NLO details for iHDM Higgs+MET model 

While a 2HDM UFO implemenfafion af NLO accuracy fo be used 
wifh MadGraph5_aMC@NLO has been made available af fhe 
end of fhe work of fhe Forum [New], in fhis work we have only 
considered LO simulafions. 

The two couplings fhaf can be changed in fhe implemented 
model follow fhe nomenclafure below: 

• Tb - fan f 

• gz - yz/ gauge coupling of Z' fo quarks 


DARK MATTER BENCHMARK MODELS FOR EARLY LHC RUN-2 SEARCHES: 


REPORT OF THE ATLAS/CMS DARK MATTER FORUM 

The other couplings are not changed, including gx (the cou¬ 
pling) which has little impact on the signal, sin a is fixed internally 
such that cos(/l — a) = 0 . The width of fhe Z' and A can be com- 
pufed aufomafically wifhin MadGraph5_aMC@NLO. The cou¬ 
plings here don'f affecf fhe signal kinematics, so fhey can be fixed 
fo defaulf values and fhen fhe signal rafes can be scaled appropri- 
afely. 

The nomenclafure for fhe masses in fhe implemenfed model is: 

• MZp - PDG ID 32 - Z' 

• MAO - PDG ID 28 - A 

• MX - PDG ID 1000022 - dark matter particle 

The ofher masses are unchanged and do nof affecf fhe resulf. 

Bofh Z' hZ{vv) and Z' hA{xx) confribufe fo fhe final sfafe, 
scaling differenf wifh model paramefers. We recommend fo gener- 
afe fhem separafely, and fhen add fhe fwo signal processes fogefher 
weighfed by cross sections. 

4.2.2 Implementation of EFT models for EW boson signatures 

The sfafe of fhe arf for fhese models is LO+PS. NLO+PS can be 
achieved as well, buf fhe corresponding implemenfafion is nof yef 
available. In our simulations we have implemenfed fhe models in 
fhe corresponding UFO files and generafed evenfs af LO via Mad- 
Grafh5_aMC@NLO 2.2.2, using pythia 8 for fhe parfon shower. 

UFO files and paramefer cards fhaf will be used as early Run-2 
benchmarks can be found on fhe Forum SVN reposifory: [Forh] for 
operafors wifh Higgs+MET final sfafes and [Fore] for W/Z/7 final 
sfafes. These models do nof require merging. 
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5 

Presentation of EFT results 


Most of this report has focused on simplified models. In this Chap¬ 
ter, we wish to emphasize the applicability of Effective Field Theo¬ 
ries (EFTs) in the interpretation of DM searches at the LHC. Given 
our current lack of knowledge about the nature of a DM particle 
and its interactions, it appears mandatory to provide the neces¬ 
sary information for a model independent interpretation of the 
collider bounds. This approach should be complemented with an 
interpretation within a choice of simplified models. We note that, 
even though EFT benchmarks are only valid in given conditions, 
the results provided by the current list of simplified models cannot 
always characterize the breadth of SM-DM interactions. In at least 
one case, composite WIMPs [NUS85; Kap92; BFTio], the contact 
interaction framework is the correct one to constrain new confine¬ 
ment scales. 

Ideally, experimental constraints should be shown as bounds of 
allowed signal events in the kinematic regions considered for the 
search, as detailed in Appendix B. A problematic situation is the 
attempt to derive a limit on nucleon-dark matter scattering cross 
sections from EFT results based on collider data ^. Experiments 
that directly probe the nucleon-dark matter scattering cross section 
are testing the regime of small momentum transfers, where the 
EFT approximation t5rpically holds. Collider experiments, though, 
are sensitive to large momentum transfers: We first illustrate the 
complications that can arise with EFTs at colliders by considering 
an effective interaction 

■^int — 1^2 ~ i‘l7!i‘l){X7^ X) 


that couples quarks and DM x fields.^ The strength of this inter- 

1 V 

action is parametrized by . A monojet signature can be 

generated from this operator by applying perturbation theory in 
the QCD coupling. An experimental search will place a limit on 
M*. For a fixed M*, a small value of g will correspond to a small 
value of A. The EFT approximation breaks down if Q > A, where 
Q is a typical hard scale of the process. The limit on small g can 
only be reliable if the kinematic region Q > A is removed from 
the event generation. However, if a fraction of events is removed 
from the prediction, the corresponding value of g must increase to 


’ Comparisons between constraints 
from different experiments meant 
to highlight their complementarity 
should be expressed as a function of 
the model parameters rather than on 
derived observables; however this is a 
point that should be developed further 
after the conclusion of the work of this 
Forum. 


’ The exact operator chosen is not 
important: as detailed in the following, 
statements concerning the applicability 
of an EFT can also be made without a 
specific relation to simplified models. 
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match the experimental limit on M*. On the other hand, if, for fhe 
same value of M*, a large A is assumed so fhaf fhe full sef of evenfs 
fulfill fhe EFT validify condifion, a larger value of g is required. For 
large enough g, compufafions based on perfurbafion fheory become 
unreliable. 

In fhe firsf parf of fhis Chapfer, we summarize two mefhods fhaf 
have been advocafed fo fruncafe evenfs fhaf do nof fulfill fhe condi¬ 
fion necessary for fhe use of an EFT. These mefhods are described 
in defail in Refs. [Bus+iqa; Bus+iqb; Bus+iqc; ATLi5d; RWZ15; 
BLWiqb]. We fhen propose a recommendafion for fhe presenfafion 
of EFT resulfs for early Run-2 LHC searches. 

5.1 Procedures for the truncation of EFT benchmark models 

5.1.1 EFT truncation using the momentum transfer and information 
on UV completion 

In fhe approach described in Ref. [Bus+iqb], fhe EFT prediction is 
modified fo incorporafe fhe effecf of a propagafor for a relatively 
lighf mediafor. For a free-level inferacfion befween DM and fhe SM 
via some mediafor wifh mass Mmed/ the EFT approximafion corre¬ 
sponds fo expanding fhe propagafor for fhe mediafor in powers of 
Q^j./ truncating at lowest order, and combining the remaining 

parameters into a single parameter M* (cormected to the scale of 
fhe inferacfion A in fhe liferafure). For an example scenario wifh a 
Z'-fype mediafor (leading fo some combinafion of operators D5 fo 
D8 in fhe nofafion of [Goo+10] for fhe EFT limif), fhis corresponds 
fo seffing 

SxSq ^ gxSq Qtr , .o ( Qtr \\_ 

<edV <ed Wmed/y 

( 5 - 1 ) 

where Qtr is the momentum carried by the mediator, and g^, yq are 
the DM-mediator and quark-mediator couplings respectively.^ A 
minimal condition that must be satisfied for fhis approximafion fo 
be valid is fhaf = gxSq^^- This requiremenf avoids 

fhe regions: Q^j. ~ ^med' which case fhe EFT misses a resonanf 
enhancemenf, and if is conservafive fo ignore fhis enhancemenf; 
and Qjj. ^ ^med' which case fhe signal cross section should 
fall according fo a power of insfead of The laffer is fhe 

problemafic kinemafic region. 

The condifion Q^j. < = gxSq^Tt^ was applied fo resfricf 

fhe kinematics of fhe signal and remove evenfs for which fhe high- 
mediafor-mass approximafion made in fhe EFT would nof be reli¬ 
able. This leads fo a smaller effective cross-section, affer imposing 
fhe evenf selection of fhe analysis. This fruncafed signal was fhen 
used fo derive a new, more conservafive limif on M* as a funcfion 
of (^ X ' gxSq )- 

For fhe example D5-like operafor, where fhe cross secfion a 
scales as there is a simple rule for converting a rescaled 


3 Here, we ignore potential complica¬ 
tions from the mediator width when 
the couplings are large. 
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cross section into a rescaled constraint on M*. if the original limit 
is based on a simple cut-and-count procedure. Defining as 
the cross section truncated such that all events pass the condition 
> Qtr, we have 

1/4 


m; 


rescaled _ | _ ^EFT j 




( 5 - 2 ) 


which can be solved for via eifher iferafion or a scan. 

Similar relafions exisf for a given UV completion of each operator. 

This procedure has been proposed in Ref. [Bus+i4b] and ifs 
applicafion to ATLAS resulfs can be found in Ref. [ATLi5d] for a 
range of operafors. We reiferafe: knowledge of fhe UV completion 
for a given EFT operator was necessary for fhis procedure; fhis 
infroduces a model-dependence fhaf was nof presenf in fhe non- 
fruncafed EFT resulfs. 

Currenfly, simplified models (including the full effecf of fhe 
mediafor propagafor) are available for comparison wifh fhe dafa, 
and since knowledge of fhe simplified models is needed for fhe 
fruncafion procedure, fhere is no reason to apply fhis prescription. 
Insfead, fhe simplified model limif for large M* can be presented 
for inferprefafion in ferms of EFT operafors. 


5.1.2 EFT truncation using the center of mass energy 

The procedure presenfed in fhe previous secfion was predicafed on 
some knowledge of fhe simplified model. This led to fhe idenfifi- 
cation of fhe mass of fhe DM pair as fhe relevant kinematic quan¬ 
tity to use in a truncation procedure. In general, if no assump¬ 
tion is made about the underlying d5mamics, it is more conser¬ 
vative to place a limit on the total center of mass energy E^m of 
fhe DM producfion process. Furfhermore, fhe direcf connecfion 
befween fhe mass scale of fhe EFT validify, Mcut/ and fhe mass 
scale fhaf normalizes fhe EFT operafor, M*, is unknown. For such 
cases, Refs.[RWZi5; BLWiqb] proposed a procedure to exfracf 
model independenf and consisfenf bounds wifhin fhe EFT fhaf can 
be applied fo any effective Lagrangian describing fhe interactions 
befween fhe DM and fhe SM. This procedure provides conservative 
limifs fhaf can be direcfly reinferprefed in any complefion of fhe 
EFT. The condition ensuring fhaf fhe EFT approximation is appro¬ 
priate is: 

Bern < Mcut ■ (5'3) 

The relationship befween Mcut and M* can be parameferized 
by an effective coupling strength such fhaf Mcut = g* M*. A 
scan over values of provides an indication of fhe sensifivify of 
fhe prediction fo fhe fruncafion procedure. In fhe Z'-type model 
considered above, is equal fo ^fgxgq- The resulfing plofs are 
shown in [RWZ15] for a parficular effective operafor. 

The advanfage of fhis procedure is fhaf fhe obfained bounds 
can be direcfly and easily recasf in any complefion of fhe EFT, by 
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computing the parameters M*, Mcut in the full model as functions 
of fhe paramefers of fhe complefe fheory. On fhe ofher hand, fhe 
resulfing limifs will be weaker fhan fhose obfained using Qtj- and a 
specific UV complefion. 


5.1.3 Truncation at the generator level 

The condifions on fhe momenfum fransfer can also be applied di- 
recfly af fhe generator level, by discarding evenfs fhaf are invalid 
and calculafing fhe limifs from fhis fruncafed shape. This pro¬ 
vides fhe necessary rescaling of fhe cross section while keeping fhe 
information on the change in the kinematic distributions due to 
the removal of fhe invalid evenfs. This procedure is more general 
wifh respecf fo rescaling the limit in the two sections above, and it 
should be followed if a search is nof simply a counting experiment 
and exploits the shapes of kinematic disfribufions. 


5.1.4 Sample results of EFT truncation procedures 

An example of fhe application of fhe two procedures fo fhe limif 
on M* from Ref. [ATLi4d] as a function of fhe producf of fhe cou¬ 
plings is shown in Figure 5.3. Only fhe region befween fhe dashed 
and fhe solid line is excluded. If can be seen fhaf fhe procedure 
from [RWZ15] ouflined in Section 5.1.2, shown in blue, is more 
conservafive fhan fhe procedure from Refs. [Bus+i4b; ATLi5d], 
described in Section 5.1.1. 



Figure 5.1: 95% CL lower limits on 
the scale of the interaction of the D5 
operator at 14 TeV, after the two trun¬ 
cation procedures. The procedure 
from [RWZ15] outlined in Section 5.1.2 
is shown in blue, while the proce¬ 
dure from Refs. [Bus+i4b; ATLi5d], 
described in Section 5.1.1 is shown 
in red. Only the region between the 
dashed and the solid lines is excluded. 
Even though the intersection between 
the two lines is not shown in this plot, 
it should be noted that no limit can be 
set anymore for sufficiently low cou¬ 
plings, whatever truncation method is 
used. 


5.1.5 Comments on unitarity considerations 

A furfher considerafion applicable fo EFT operators af hadron 
colliders is fhe pofential violafion of unitarify. An analysis of fhe 
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operator 


q'y}iqxi^x 

Ml 


provides the limit: 


M* > [i{s)Vs 



(54) 


where ^/s is (maximally) the collider energy and / 5 (s) is the DM 
velocity [SV12]. Constraints for other operators have also been de¬ 
rived [EY14]. This constraint on M* still is open to interpretation, 
since the relation to Mcut is not resolved, except for a specific sim¬ 
plified model. Derived limifs on M* should be compared fo fhis 
unifarify bound fo check for consisfency. 


5.2 Recommendation for presentation of EFT results 

In fhis report, we make two recommendations for the presenta¬ 
tion of collider resulfs in ferms of Effecfive Eield Theories for fhe 
upcoming Run-2 searches. A full discussion of fhe presenfation of 
collider resulfs in relafion fo ofher experimenfs is leff fo work be¬ 
yond fhis Eorum, where ATLAS, CMS, fhe fheory communify and 
the Direct and Indirect Detection communities are to be involved. 

We divide the EET operators in two categories: those that can 
be mapped to one or more UV-complete simplified models, such 
as fhose commonly used in LHC searches so far and defailed 
in [Goo+10], and fhose for which no UV complefion is available 
fo LHC experimenfs, such as fhose ouflined in Secfion 3.2. 


5.2.1 EFT benchmarks with corresponding simplified models 

If a simplified model can be mapped fo a given EET, fhen fhe 
model's high-mediafor-mass limif will converge fo fhe EET. 

A sfudy of 14 TeV benchmarks for narrow resonances wifh 
= 0.25 and gx = r (see Secfion 2.1.1) shows fhaf a mediafor wifh a 
mass of af leasf 10 TeV fully reproduces fhe kinemafics of a confacf 
inferacfion and has no remaining dependence on fhe presence of 
a resonance. A comparison of fhe main kinematic variables for 
fhe s-channel vector mediafor model wifh a widfh of 0.1 M^ed is 
shown in Eig. 5.2.4 

As already observed in Secfion 2.1.1, varying fhe DM mass 
changes fhe kinemafics, bofh in fhe simplified model and in fhe 
EET case. This can be seen in Eig. 5.3. 


+ The use of a fixed width rather than 
the minimal width is exclusive of these 
plots. 


Based on fhese sfudies, fhe Eorum recommends experimenfal 
collaborafions fo add one grid scan poinf af very high media¬ 
for mass (10 TeV) fo fhe scan, for each of fhe DM masses for fhe 
s-channel simplified models described in Secfion 2. This will allow 
fo reproduce fhe resulfs of an equivalenf confacf inferacfion as a 
simple exfension of fhe exisfing parameter scan. 

If should be checked fhaf fhe high-mass mediafor case for fhe 
simplified model is correcfly implemented 
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(a) Hr 



(b) Center of mass energy E„„ 



(d) Leading DM transverse momentum 



(c) Mediator transverse momentum 



(e) DM transverse momentum (sub-leading) 

Figure 5.2: Comparison of the kine¬ 
matic distributions at 14 TeV between 
a narrow s-channel mediator and the 
corresponding D5 contact operator, at 
generator level for a jet+Er signature. 
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5.2.2 EFT benchmarks with no corresponding simplified models 

Whenever a UV completion is not available, an EFT still captures 
a range of possible theories beyond the simplified models fhaf we 
already consider. However, in fhe case of fhe dimension-7 operafors 
defailed in Secfion 3.2 we can only roughly confrol how well fhe 
EFT approximafion holds, as described in Secfion 3.2.4. Despife fhe 
facf fhaf a propagafor was infroduced fo mofivafe fhe fruncation 
procedure for s-channel models, fhe prescripfion from Sec. 5.2.1 
depends upon fhe simplified model fo derive fhe energy scaling 
thaf is used for fhe comparison wifh fhe momenfum fransfer. The 
simple facf remains fhaf fhe effecfive coupling of fhe operafor - 
g/ A” - should nof allow momenfum flow Q > A or y > 4/1. Given 
our ignorance of fhe acfual kinematics, the truncation procedure 
recommended for fhis purpose is fhe one described in Secfion 5.1.2, 
as if is independenf from any UV complefion defails. 

Because fhere is no UV complefion, fhe paramefer Mcut can be 
freafed more freely fhan an explicif function of g and A. If makes 
sense fo choose Mcut such fhaf we idenfify fhe fransifion region 
where fhe EFT sfops being a good description of UV complefe 
fheories. This can be done using fhe ratio R, which is defined as 
fhe fracfion of evenfs for which s > For large values of Mcut/ 

no evenfs are fhrown away in fhe fruncation procedure, and R = 1. 

As Mcut becomes smaller, evenfually all evenfs are fhrown away in 
the truncation procedure, i.e. R = o, and the EFT gives no exclusion 
limits for fhe chosen accepfance. 

We propose a rough scan over Mcut/ such fhaf we find fhe values 
of Mcut for which R ranges from 0.1 fo 1. The analysis can fhen 
perform a scan over several values of Mcut/ and show fhe fruncafed 
limif for each one of fhem. 
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(b) Invariant mass of the two WIMPs D5 operator 




(c) Bt, simplified model 


(d) Invariant mass of the two WIMPs simplified model 


Figure 5.3: Comparison of the kine¬ 
matic distributions for a narrow 
s-channel mediator, at generator level 
for a jet+^r signature, for varying DM 
masses. 
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Evaluation of signal theoretical uncertainties 


A comprehensive and careful assessment of signal theoretical un¬ 
certainties plays in general a more important role for the back¬ 
ground estimations (especially when their evaluation is non-entirely 
data-driven) than it does for signal simulations. Nevertheless, also 
for signal samples theoretical uncertainties are relevant, and may 
become even dominant in certain regions of phase space. 

The uncertainties on the factorization and renormalization scales 
are assessed by the experimental collaborations by varying the orig¬ 
inal scales of the process by factors of 0.5 and 2. The evaluation of 
the uncertainty on the choice of PDF follows the PDF4LHC recom¬ 
mendation [Pdf] of considering the envelope of different PDF error 
sets, in order to account for the uncertainty on the various PDFs as 
well as the uncertainty on the choice of the central value PDF. The 
Forum has not discussed the uncertainties related to the merging 
of different samples, nor the uncertainty due to the choice of the 
modeling of the parton shower. This Chapter provides technical de¬ 
tails on how scale and PDF uncertainties can be assessed for events 
generated with powheg and MadGraph 5 _aMC@NLO. 

6.1 POWHEG 

When using powheg [FNO07; Ali+10; Naso4], it is possible to study 
scale and PDF errors for the dark matter signals. A fast reweighting 
machinery is available in powheg-box that allows one to add, after 
each event, new weights according to different scale or PDF choices, 
without the need to regenerate all the events from scratch. 

To enable this possibility, the variable storeinfo_rwgt should be 
set to 1 in the powheg input file when the events are generated for 
the first time^. After each event, a line starting with 

#rwgt 

is appended, containing the necessary information to generate extra 
weights. In order to obtain new weights, corresponding to different 
PDFs or scale choice, after an event file has been generated, a line 

compute-rwgt 1 

should be added in the input file along with the change in param¬ 
eters that is desired. For instance, renscfact and facscfact allow 


’ Notice that even if the variable is not 
present, by default it is set to 1. 
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one to study scale variations on the renormalization and factoriza¬ 
tion scales around a central value. By rurming the program again, a 
new event file will be generafed, named <OriginalName>- rwgt . Ihe, 
wifh one more line af fhe end of each evenf of fhe form 

#new weight,renfact,facfact,pdfl,pdf2 

followed by five numbers and a characfer sfring. The firsf of fhese 
numbers is fhe weighf of fhaf evenf wifh fhe new paramefers cho¬ 
sen. By running in sequence fhe program in fhe reweighfing mode, 
several weighfs can be added on fhe same file. Two remarks are in 
order. 

• The file wifh new weighfs is always named 
<OriginalName>-rwgt.the 

hence care has to be taken to save it as 
<OriginalName>.the 

before each iteration of the reweighfing procedure. 

• Due to the complexity of the environment where the program 
is likely to be run, it is strongly suggested as a self-consistency 
check that the first reweighfing is done keeping the initial pa¬ 
rameters. If the new weights are not exactly the same as the 
original ones, then some inconsistency must have happened, or 
some file was probably corrupted. 

It is possible to also have weights written in the version 3 Les 
Houches format. To do so, in the original run, at least the token 
lhrwgt_id 'ID' 

must be present. The reweighting procedure is the same as de¬ 
scribed above, but now each new run can be tagged by using a 
different value for the lhrwgt_id ke5rword. After each event, the 
following lines will appear: 

<rwgt> 

<wgt id='ID'> 

<wgt id='IDl'> 

</rwgt> 

A more detailed explanation of what went into the computation 
of every single weight can be included in the <header> section of 
the event file by adding/changing the line 
lhrwgt_descr 'some info' 

in the input card, before each "reweighting" run is performed. 

Other useful ke5rwords to group together different weights are 
lhrwgt_group_name and lhrwgt_group_combine. 

More detailed information can be obtained by inspecting the 
document in /Docs/V2- paper, pdf under the common powheg-box- 
V2 directory. 

6.2 The SysCalc package in MadGraph5_aMC@NLO 

SysCalc is a post-processing package for parton-level events as ob¬ 
tained from leading-order calculations in MadGraph5_aMC@NLO. 
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It can associate to each event a series of weights corresponding to 
the evaluation of a cerfain class of fheorefical uncerfainties. The 
event files in inpuf and oufpuf are compliant with the Les Houches 
V3 format. For NLO calculations, PDF and scale uncertainties are 
instead evaluated automatically by setting corresponding instruc¬ 
tions in the run_card.dat and no post-processing is needed (or 
possible). 

The requirements of fhe package as inpufs are : 

• A sysfemafics file (which can be generafed by MadGraph 5 v. 

1.6.0 or lafer) [Alw+14; Alw+ii]. 

• The Pyfhia-PGS package (v. 2.2.0 or lafer) [SMS06]. This is 
needed only in fhe case of mafching scales variations. 

• The availabilify of LHAPDF5 [WBG05]. 

• A configurafion file (i.e. a fexf file) specifying fhe paramefers fo 
be varied. 

SysCalc supporfs all leading order compulations generafed 
in MadGraph 5 _aMC@NLO including fixed-order compufa- 
fion and mafched-merged compulation performed in fhe MLM 
scheme [Man+07]. MadGraph 5 _aMC@NLO stores additional 
information inside fhe even! in order fo have access fo all fhe infor¬ 
mation required fo compufe fhe convolution of fhe PDFs wifh fhe 
mafrix elemenf for fhe various supporfed sysfemafics. 

Below follows an example configurafion file which could serve as 
an example: 

# Central scale factors 
scalefact: 

0.5 1 2 

# Scale correlation 

# Special value -1: all combination (N**2) 

# Special value -2: only correlated variation 

# Otherwise list of index N*fac_index + ren_index 

# index starts at 0 
scalecorrelation: 

-1 

# tts emission scale factors 
alpsfact: 

0.5 1 2 

# matching scales 
matchscale: 

30 60 120 

# PDF sets and number of members (optional) 

PDF: 

CTlO.LHgrid 53 
MSTW2008nlo68cl.LHg rid 


Wifhouf mafching/merging, SysCalc is able fo compufe fhe 
variation of renormalisafion and facforisafion scale (paramefer 
scalefact) and fhe change of PDFs. The variation of fhe scales 
can be done in a correlafed and / or uncorrelafed way, basically 
following fhe value of fhe scalecorrelation parameter which can 
take the following values: 

• -1 : to account for all combinations. 
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• -2 : to account only for the correlated variations. 


• A set of positive values corresponding fo fhe following enfries 
(assuming 0.5, i, 2 for fhe scalefact enfry): 


o: = hT^/2, Hr = H°r^/2 

1: HF = Hf^/ 2 ,HR = Hr 

2 : HF = F7'^/^>m = Fr*2 

„ _ ..“"g „ _ „°rig /9 

3. JAp — Up , jiR — 12 . 

4: JAp = 

5 : HF=Hr, hr = nr *2 

6: hf = hT^ *2,hr = H°r^/2 
7: HF = Hf^* 2 ,HR = Hr 
8: Hf = hT^ *2,Hr = hT^ * 2 


Wifhouf correlation, fhe weighf associafed fo fhe renormalisafion 
scale is fhe following: 


^ *^new 


^(Fr) 


(6.1) 


where A is fhe scale variation considered, Worig and Wnew are re- 
specfively fhe original/new weighfs associafed fo fhe evenf. N is 
fhe power in fhe sfrong coupling for fhe associafed evenf (inferfer- 
ence is nof faken accounf on an evenf by evenf basis). The weighf 
associafed fo fhe scaling of fhe factorisation scale is: 




/l,orig('^l/^ * Ff) */2,orig(■^2/^ * Ff) 
/l,orig(^l/ Ff) * /2,orig(^2/ Ff) 


* Worig, 


(6.2) 


where /i,orig are fhe probabilifies from fhe original PDF sef asso¬ 
ciafed fo fhe incoming parfons, which hold a proton momenfum 
fraction xi and X2 for fhe firsf and second beam respectively 
The variations for fhe PDF are given by fhe corresponding 
weighfs associafed fo fhe new PDF sefs: 


W, 


iPDF 


/l,new(^l,FF) */2 ,new {X2,Ff) 
/^l,orig(^l/F f) */2,orig(^2 /Ff) 


* Worig, 


(6.3) 


where /i^new is the new PDF probability associated to parton i. 

In presence of mafching, MadGraph 5 _aMC@NLO associafes 
one hisfory of radiafion (inifial and / or final sfafe radiation) ob- 
fained hy a kj clustering algorifhm, and calculates ccg af each verfex 
of fhe hisfory fo a scale given by fhe aforemenfioned clusfering 
algorifhm. Furfhermore, MadGraph 5 _aMC@NLO reweighfs 
fhe PDF in a fashion similar fo whaf a parfon shower would do. 
SysCalc can perform fhe associafed re-weighfing (paramefer 
alpsfact) by dividing and multiplying by fhe associafed facfor. 

For each sfep in fhe hisfory of fhe radiafion (associafed fo a scale 
jii = kji), fhis corresponds fo fhe following expression for a Final 
Sfafe Radiafion (FSR): 


W, 


FSR 


«s(^*Fi) 

XsiFi) 


* Worig, 


(6.4) 
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and to the following expression for Initial Sfafe Radiation (ISR), 
associafed fo a scale and fraction of energy x,: 


W. 


ISR 


Its (A * jij) 

as(F/) 


fbi^um+i) 


* Worig, 


(6.5) 


where is fhe scale of fhe nexf sfep in fhe (inifial sfafe) history 
of radiafion. 

SysCalc can include fhe weighf associafed fo differenf merging 
scales in fhe MLM mafching/merging mechanism (for outpuf of fhe 
pythia6 package or pythia-pgs package). 

In fhaf case, fhe parfon shower does nof vefo any evenf accord¬ 
ing fo fhe MLM algorifhm, alfhough in fhe outpuf file fhe scale of 
fhe firsf emission is refained. Having fhis information, SysCalc 
can fesf each value of fhe specified mafching scales under fhe 
matchscale paramefer block. SysCalc will fhen fesf for each of fhe 
values specified in fhe paramefer matchscale if fhe evenf passes fhe 
MLM criferia or nof. If if does nof, fhen a zero weighf is associafed 
fo fhe evenf, while if if does, fhen a weighf 1 is kepf. As a reminder, 
fhose weighfs are fhe equivalenf of having a (approximate) Sudakov 
form-facfor and removing af fhe same time fhe double counting 
befween fhe evenfs belonging fo differenf mulfiplicifies. 

Finally, we give an example of fhe SysCalc outpuf which fol¬ 
lows fhe LHEF V3 formaf. The following block appears in fhe 
header of fhe oufpuf file: 


<heacler> 

<initrwgt> 

<weightgroup type="Central scale variation" combine="envelope"> 
<weight id="l"> mur=0.5 muf=0.5 </weight> 

<weight id="2"> mur=l muf=0.5 </weight> 

<weight id="3"> mur=2 muf=0.5 </weight> 

<weight id="4"> mur=0.5 muf=l </weight> 

<weight id="5"> mur=l muf=l </weight> 

<weight id="6"> mur=2 muf=l </weight> 

<weight id="7"> mur=0.5 muf=2 </weight> 

<weight id="8"> mur=l muf=2 </weight> 

<weight id="9"> mur=2 muf=2 </weight> 

</weightgroup> 

<weightgroup type="Emission scale variation" combine="envelope"> 
<weight id="10"> alpsfact=0.5</weight> 

<weight id="ll"> alpsfact=l</weight> 

<weight id="12"> alpsfact=2</weight> 

</weightgroup> 

<weightgroup type="CT10nlo.LHgrid" combine="hessian"> 

<weight id="13">Member 0</weight> 

<weight id="14">Member l</weight> 

<weight id="15">Member 2</weight> 

<weight id="16">Member 3</weight> 


<weight id="65">Member 52</weight> 
</weightgroup> 

</initrwgt> 

</header> 

For each event, the weights are then written as follows: 
<rwgt> 

<wgt id="l">83214.7</wgt> 

<wgt id="2">61460</wgt> 
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<wgt id="3">47241.9</wgt> 
<wgt id="4">101374</wgt> 

<wgt id="64">34893.5</wgt> 
<wgt id="65">41277</wgt> 

</rwgt> 
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Conclusions 


The ATLAS/CMS Dark Matter Forum concluded its work in June 
2015. Its mandate was focused on identifying a priorifized, com- 
pacf sef of simplified model benchmarks fo be used for fhe design 
of fhe early Run-2 LHC searches for Bj +X final sfafes. Ifs parfici- 
panfs included many of fhe experimenfers from bofh collaborafions 
fhaf are involved in fhese searches, as well as many of fhe fheorisfs 
working acfively on fhese models. This reporf has documenfed fhis 
basis sef of models, as well as sfudies of fhe kinemafically-disfincf 
regions of fhe paramefer space of fhe models, fo aid fhe design of 
fhe searches. Table 6.1 summarizes fhe sfafe of fhe arf of fhe cal¬ 
culations, even! generafors, and fools fhaf are available fo fhe two 
LHC collaborafions fo simulafe fhese models af fhe sfarf of Run-2. 

If also describes some fhaf are known fo be under developmenf as 
fhe reporf was finalized. 

This documenf primarily presenfs sfudies relafed fo simplified 
models. The presenfafion of resulfs for EFT benchmark models is 
also discussed. The sfudies confained in fhis reporf are meanf fo 
highlighf fhe use of EFTs as a benchmark fhaf is complemenfary 
fo simplified models, and fo demonsfrafe how fhaf collider resulfs 
could be presenfed a function of fhe fraction of evenfs fhaf are valid 
wifhin fhe confacf inferacfion approximafion. 

A number of poinfs remain fo be developed beyond fhe scope 
of fhis Forum, in order fo fully benefif from LHC searches in fhe 
global quesf for Dark Matter. Firsf and foremosf, fo accomodafe 
fhe urgenf need of a basis sef of simplified models, fhis work has 
made many grounding assumpfions, as sfafed in fhe infroducfion. 
Deparfures from fhese assumpfions have nof been fully explored. 
As a consequence, fhe lisf of models and implemenfafions em¬ 
ployed by fhe ATLAS and CMS collaborafions for early LHC Run-2 
searches is nof meanf fo exhausf fhe range of possibilifies for medi- 
afing processes, lef alone cover all plausible mdoels of collider dark 
matter production. Rafher, if is hoped fhaf ofhers will continue fhe 
sysfemafic exploration of fhe mosf generic possibilifes for collider 
dark matter producfion, building upon fhe framework used in fhis 
reporf jusf as fhis reporf has relied heavily on fhe work of many 
ofhers. This also applies fo models fhaf exisf in liferafure buf do 
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Benchmark models for ATLAS and CMS Run-2 DM searches 

vector/axial vector mediator, s-channel (Sec. 2.1) 

Signature 

State of the art calculation and tools 

Implementation 

References 

jet + Bt 

NLO+PS (powheg, SVN r3059) 

NLO+PS (DMsimp UFO + MADGRArH5_AMC@NLO V2.3.0) 

NLO (mcfm V7.0) 

[Fori; Foro] 

[New] 
Upon request 

[HKR13; HR15; Ali+io; Naso4; 
FNO07] 

[Alw+14; All+14; Deg+12] 

[FW13; Har+15] 

W/Z /7 + Bt 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 

NLO+PS (DMsimp UFO + MADGRArH5_AMC@NLO V2.3.0) 

[Fora] 

[New] 

[Alw+14; All+14; Deg+12] 
[Alw+14; All+14; Deg+12] 

scalar/pseudoscalar mediator, s-channel (Sec. 2.2) 

Signature 

State of the art calculation and tools 

Implementation 

References 

jet + Bt 

LO+PS, top loop (powheg, r3059) 

LO+PS, top loop (DMsimp UFO + MadGraph5_aMC@NLO v.2.3.0) 
LO, top loop (mcfm V7.0) 

[Forn; Form] 

[New] 

Upon request 

[HKR13; HR15; Ali+io; Naso4; 
FNO07] 

[Alw+14; Hir+11; All+14; 
Deg+12] 

[FW13; Har+15] 

W/Z /7 + Bt 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 


[Alw+14; All+14; Deg+12] 

tt, bb+ Bt 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 

NLO+PS (DMsimp UFO + MADGRArH5_AMC@NLO V2.3.0) 

[Ford] 

[New] 

[Alw+14; All+14; Deg+12] 
[Alw+14; All+14; Deg+12] 


scalar mediator, f-channel (Sec. 2.3) 


Signature 


State of the art calculation and tools 


Implementation References 


jet(s) + £t (2-quark gens.) 
)et(s) + Bt (3-quark gens.) 

W/Z/j + Bt 
b + Bt 


LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) 

Specific simplified models with EW bosons (Sec. 3.1) 


[Forj] 

[Fori] 

TBC 

[Forg] 


[PVZ14; Alw+14; All+14; 
Deg+12] 

[Bel+12; Alw+14; All+14; 
Deg+12] 

[Bel+12; Alw+14; All+14; 
Deg+12] 

[LKW13; Agr+i4b; Alw+14; 
All+14; Deg+12] 


Signature and model 
Higgs + Bt, vector med. 
Higgs + Bt, scalar med. 
Higgs + Bt, 2HDM 

Signature and model 

W/Zfj + Bt, dim-7 

Higgs + Bt, dim-4/dim-5 
Higgs + Bt, dim-8 


State of the art calculation and tools Implementation 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) [Forh] 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) [Forh] 

LO+PS (UFO + MadGraph5_aMC@NLO V2.2.3) [Forb] 

Contact interaction operators with EW bosons (Sec. 3.1) 

State of the art calculation and tools Implementation 


LO+PS (UFO + MadGraph5_aMC@NLO 

LO+PS (UFO + MadGraph5_aMC@NLO 
LO+PS (UFO + MadGraph5_aMC@NLO 


V2.2.3) 

[Fore] 

V2.2.3) 

[Fore] 

V2.2.3) 

[Forh] 


References 

[Car+14; BLWi4b; Alw+14; 
All+14; Deg+12] 

[Car+14; BLWi4b; Alw+14; 
All+14; Deg+12] 

[BLWi4b; Alw+14; All+14; 
Deg+12] 


References 

[Cot+13; Car+13; CHH15; 
BLWi4b; Alw+14; All+14; 
Deg+12] 

[Car+14; PS14; BLWi4b; 
Alw+14; All+14; Deg+12] 
[Car+14; PS14; BLWi4b; 
Alw+14; All+14; Deg+12] 


Table 6.1: Summary table for available benchmark models considered within the works of this Forum. 
The results in this document have been obtained with the implementations in bold. 
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not have an implementation yet: we hope that this work will fur¬ 
ther encourage the theory and generator community to improve the 
implementation of new models as well as fhe precision of fhe cal¬ 
culations of exisfing ones. The role of consfrainfs on fhe mediafor 
particles from direcf pasf and presenf collider searches should also 
be developed furfher. 

Furfhermore, we see fhe need for broader discussion on fhe com¬ 
parison of experimenfal resulfs amongsf collider and non-collider 
searches for particle dark matter. This poinf will have fo be ad¬ 
dressed before fhe presenfafion of Run-2 resulfs: The uncerfainfies 
in fhe comparisons befween experimenfs should be discussed and 
conveyed, so fhaf fhe differenf resulfs can be placed in fheir correcf 
confexf, and so we can collecfively build a fair and comprehensive 
picfure of our undersfanding of particle Dark Matter. 
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A 

Appendix: Additional models for Dark Matter searches 


A.i Models with a single top—quark + tj 

Many different theories predict final states with a single top and 
associated missing transverse momentum (monotop), some of them 
including dark matter candidates. A simplified model encompass¬ 
ing the processes leading to this phenomenology is described in 
Refs. [AFMii; Agr+i 4 a; Bou+ 15 ], and is adopted as one of the 
benchmarks for Run 2 LHC searches. 

The simplified model is constructed by imposing that the model 
Lagrangian respects the electroweak SU{2)i x Lt(l)y gauge sym¬ 
metry and by requiring minimality in terms of new states to sup¬ 
plement to the Standard Model fields. As a result, two monotop 
production mechanisms are possible. In the first case, the monotop 
system is constituted by an invisible (or long-lived with respect to 
detector distances) fermion x and a top quark. It is produced as 
shown in the diagram of A.i (a) where a colored resonance cp lying 
in the triplet representation of SU{3)c decays into a top quark and 
a X particle. In the second production mode, the monotop state is 
made of a top quark and a vector state V cormected to a hidden 
sector so that it could decay invisibly into, e.g., a pair of dark mat¬ 
ter particles as studied in [Bou+ 15 ]. The production proceeds via 
flavor-changing neutral interactions of the top quark with a quark 
of the first or second generation and the invisible V boson (see the 
diagrams of A.i (b) and (c)). 

Resonant production 

In this case, a colored 2/3-charged scalar (cp) is produced and 
decays into a top quark and a spin-1 / 2 invisible particle, x- The dy¬ 
namics of the new sector is described by the following Lagrangian: 


C = 


cpd‘' 


^SR ■ 


^SrT5 


cpu 


«SR + ^SR 75 


it + h.c. 


(A.I) 


where u (d) stands for any up-type (down-type) quark, the nota¬ 
tion SR refers to the monotop production mechanism via a scalar 
resonance and all flavor and color indices are understood for clarity. 

In the notation of [Agr+iqa], the couplings of the new colored 
fields to down-t 5 rpe quarks are embedded into the 3x3 antisym- 
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Figure A.i: Feynman diagrams of 
leading order processes leading to 
monotop events: production of a 
colored scalar resonance (p decaying 
into a top quark and a spin-1 / 2 
fermion x (a)/ s— (b) and t-channel 
(c) non resonant production of a top 
quark in association with a spin-i 
boson V decaying invisibly. 


metric matrices flgjj (scalar couplings) and (pseudoscalar cou¬ 
plings) while those to the new fermion x and one single up-type 
quark are given by the three-component vectors and in 
flavor space. 

Under fhe form of Eq. (A.i), fhe Lagrangian is fhe one infro- 
duced in fhe original monofop search proposal [AFMii]. If has 
been used by fhe CMS collaborafion for Run I analyses affer ne¬ 
glecting all pseudoscalar componenfs of fhe couplings and adding 
fhe vecfor resonance case for which minimalify requiremenfs 
are difficulf fo accommodafe [CMSi5d]. In confrasf, fhe sfudy 
of Ref. [B0U-F15] has imposed elecfroweak gauge invariance and 
required minimalify. This enforces all new couplings fo be righf- 
handed so fhaf 

^SR ~ ^SR ~ 'jVs arid , (A.2) 

where fhe objecfs j/s and Ag are a fridimensional vecfor and a 3 x 3 
mafrix in flavor space respectively This class of scenarios is fhe 
one fhaf has been adopfed by fhe ATLAS collaborafion for ifs Run I 
monofop searches [ATLi5b] and will be considered by bofh collabo¬ 
rations for Run II analyses. 

The resulting model can be likened fo fhe MSSM wifh an R- 
parify violafing of a fop squark fo fhe Sfandard Model down-f5rpe 
quarks and an R-parify conserving inferacfion of a fop quark and a 
fop-squark fo a neufralino. 


Non-Resonant production 

For non-resonanf monofop producfion, fhe monofop sfafe is 
produced via flavor-changing neufral inferacfions of fhe fop quark, 
a lighfer up-f5rpe quark and a new invisible vecfor parficle V. This 
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is the only case considered, as having a new scalar would involve 
in particular a mixing with the SM Higgs boson and therefore a 
larger number of free paramefers. The Lagrangian describing fhe 
d5mamics of fhis non-resonanf monofop producfion case is: 


£ = 




apc+blcld 


u + h.c. 


(A.3) 


where fhe flavor and color indices are again understood for clarify. 
The sfrengfh of fhe inferacfions among fhese fwo sfafes and a pair 
of up-fype quarks is modeled via fwo 3 x 3 mafrices in flavor space 
for fhe vector couplings and bp^. for fhe axial vecfor couplings, 
fhe FC subscripf referring fo fhe flavor-changing neufral monofop 
producfion mode and fhe (1) superscripf fo fhe vecforial nafure of 
fhe invisible parficle. 

As for fhe resonanf case, fhe Lagrangian of Eq. (A.3) is fhe one 
fhaf has been used by CMS after reinfroducing fhe scalar opfion 
for fhe invisible sfafe and neglecfing all pseudoscalar inferac¬ 
fions [CMSi5d]. As already menfioned, a simplified sefup mofi- 
vafed by gauge invariance and minimalify has been preferred so 
thaf, as shown in Ref. [Bou+15], we impose all inferacfions fo in¬ 
volve right-handed quarks only, 

apc = bpQ = (A.4) 

where ar denotes a 3 x 3 matrix in flavor space. This implies fhe 
vecfor field fo be an 511 ( 2 )^ singlef. 


Model parameters and assumptions 

The models considered as benchmarks for fhe firsf LHC searches 
confain furfher assumpfions in terms of fhe flavor sfrucfure of fhe 
model wifh respecf fo fhe Lagrangians of fhe previous subsecfion. 

In order fo have an observable monofop signafure af fhe LHC, fhe 
Lagrangians infroduced above musf include nof too small couplings 
of fhe new parficles fo firsf and second generafion quarks. For 
simplicify, we assumed fhaf only charmels enhanced by parfon 
densify effecfs will be considered, so fhaf we fix 

(«r)i3 = («r)31 = « / , . , 

(A.5) 

(As)i2 = -(As) 21 = A and [y,)^=y, 

all ofher elemenfs of fhe matrices and vectors above being set to 
zero. 


Implementation In order to allow one for fhe Monte Carlo sim- 
ulafion of evenfs relevant for fhe monofop producfion cases de¬ 
scribed above, we consider fhe Lagrangian 


£ = 


aVpuYPRi + ^cpd‘^PRS + yfxPR^ + h-c., 


(A.6) 


where Fr sfands for fhe righf-handed chiralify projector and fhe 
new physics couplings are defined by fhe fhree paramefers a, A and 
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y. We additionally include a coupling of the invisible vector boson 
y to a dark sector (represented by a fermion \p) whose strength can 
be controlled through a parameter yoM/ 

-c = ■ (A.7) 

This ensures the option to make the V-boson effectively invisible by 
tuning gDM respectively to a. We implement the entire model in the 
FeynRules package [All+14] so that the model can be exported to a 
UFO library [Deg+12] to be linked to MadGraph5_aMC@NLO [Alw+14] 
for event generation, following the approach outlined in [Chr+11]. 

A. 1.1 Parameter scan 

Under all the assumptions of the previous sections, the parameter 
space of the resonant model is defined by four quantities, namely 
the mass of the new scalar field q), the mass of the invisible fermion 
X and the strengths of the interactions of the scalar resonance with 
the monotop system y and with down-t5rpe quarks A. One of both 
coupling parameters could however be traded with the width of the 
resonance. 

The parameter space of the non-resonant model is defined by 
two parameters, namely the mass of the invisible state V and its 
flavor-changing neutral coupling to the up-t5q)e quarks ar. 

In the case of the non-resonant model, the invisible vector is 
cormected to a hidden sector that could be, in its simplest form, 
parameterized by a new fermion [Bou+15]. This has effects on the 
width of the invisible V state. 

A consensus between the ATLAS and CMS collaborations has 
been reached in the case of non-resonant monotop production. The 
results have been described above. In contrast, discussions in the 
context of resonant monotop production are still on-going. The 
related parameter space contains four parameters and must thus 
be further simplified for practical purposes. Several options are 
possible and a choice necessitates additional studies that will be 
achieved in a near future. 

It has been verified that the kinematics do not depend on the 
width of the invisible state in the case where this width is at most 
10% of the U-mass. This is illustrated in Fig. A.2, where we show 
the transverse-momentum spectra of the U-boson when it decays 
into a top-up final state and for different U-boson masses. The 
results are independent of the visible or invisible decay modes as 
we are only concerned with the kinematic properties of the invisible 
state. 

A.1.2 Single Top Model implementation 

Card files for MadCraph5_aMC@NLO are provided on the Fo¬ 
rum SVN repository [Fort] and correspond to the Lagrangian that 
has been implemented in FeynRules. Each coupling constant of the 
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[GeV] p;? [GeV] p^! [GeV] 


model can be set via the block COUPX of the parameter card. Its en¬ 
tries 1, 2 and 3 respectively correspond to the monotop-relevant pa¬ 
rameters a, A and y, while the width (and in particular the invisible 
partial width) of fhe V-boson can be funed via fhe gum paramefer 
fo given in fhe enfry lo of fhe COUPX block. 

The masses of fhe particles are set in the MASS block of fhe pa¬ 
ramefer card, fhe PDG codes of fhe new sfafes being 32 (fhe vecfor 
sfafe V), 1000006 (fhe <y colored resonance), 1000022 (fhe invisible 
fermion x) and 1000023 (fhe fermion ip connecfing fhe V sfafe fo fhe 
dark secfor). The widfh of fhe new vecfor has fo be compufed from 
all open free-level decays (affer fixing yoM to a large value and sef- 
ting the relevant entry to Auto in the DECAY block of fhe paramefer 
card), while fhe way fo calculafe fhe widfh of fhe resonance (p is 
under discussion by bofh the ATLAS and CMS collaborations. The 
chi and psi fermions are faken sfable so fhaf fheir widfh vanishes. 


Figure A.2: 

Distributions of the transverse 
momentum of the V boson in the case 
of the process pp ^ tV ^ t[tu + c.c.). 
We have imposed that the V-boson is 
produced on-shell and have chosen 
its mass to be my — 200, 600 and 
1000 GeV (left, central and right 
panels). We have considered three 
possible cases for the total width of 
the V-boson, which has been fixed to 
0.61%, 0.1% and 10% of the mass. 


A.2 Further W+tj models with possible cross-section enhance¬ 
ments 


As pointed out in Ref. [Bel+i5b], fhe mono-W signafure can probe 
fhe iso-spin violafing inferacfions of dark maffer with quarks. The 
relevant operator after the electroweak symmetry breaking is 

^XlpX {uLjf'uL + . (A.8) 

Here, we only keep the left-handed quarks because the right- 
handed quarks do not radiate a W-gauge boson from fhe weak 
inferacfion. As fhe LHC consfrains fhe cufoff fo higher values, if 
is also imporfanf fo know fhe corresponding operators before fhe 
elecfroweak symmefry. Af fhe dimension-six level, fhe following 
operafor 

^XIpXQlA'Ql (A.q) 

conserves iso-spin and provides us f = 1 [Bel+i5b]. Af fhe dimension- 
eighf level, new operators appear fo induce iso-spin violafion and 
can be 

^XJpXiHQLlinQiH^) + ^XlpXiHQihnQLH^) ■ (A.io) 
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After inputting the vacuum expectation value of the Higgs field, we 
have 

^6 + 




C6 + clvl^/lA^ 


(A.ii) 


‘'8 ‘^EW/ 

For a nonzero Ce and dew ^ A, the iso-spin violation effects are 
suppressed. On the other hand, the values of Ce, Cg and Cg depend 
on the UV-models. 

There is one possible UV-model to obtain a zero value for cg and 
non-zero values for Cg and Cg. One can have the dark matter and 
the SM Higgs field charged under a new 11 ( 1 )' symmetry. There 
is a small mass mixing between SM Z-boson and the new Z' with 
a mixing angle of 0 (d|^/M|,). After integrating out Z', one has 
different effective dark matter couplings to and di fields, which 
are proportional to their couplings to the Z boson. For this model, 
we have cg = 0 and 




- 2 -F 5 sin^ Fw 




- 2.7 


(A.12) 


and order of unity. 


A.3 Simplified model corresponding to dimension-^ EFT oper¬ 
ator 

As an example of a simplified model corresponding to the dimension- 
5 EFT operator described in Section 3.2, we consider a Higgs portal 
with a scalar mediator. Models of this kind are among the most 
concise versions of simplified models that produce couplings of 
Dark Matter to pairs of gauge-bosons. Scalar fields may couple di¬ 
rectly to pairs of electroweak gauge bosons, but must carry part of 
the electroweak vacuum expectation value. One may thus consider 
a simple model where Dark Matter couples to a a scalar singlet 
mediator, which mixes with the fields in the Fliggs sector. 

Lc ^MsS^ a'S\H\^-\- ySxx (^.13) 

Where H is a field in the Higgs sector that contains part of the 
electroweak vacuum expectation value, S is a heavy scalar singlet 
and 3;; is a Dark Matter field. There is then an s-charmel diagram 
where DM pairs couple to the singlet field S, which then mixes 
with a Higgs-sector field, and couples to W and Z bosons. This 
diagram contains 2 insertions of EW symmetry breaking fields, 
corresponding in form to the effective dimension-5 operator in 
Section 3.2.1. 


A.4 Inert two-Higgs Doublet Model (IDM) 

Eor most of the simplified models included in this report, the mass 
of the mediator and couplings/width are non-trivial parameters of 
the model. In these scenarios, we remain agnostic about the theory 
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behind the dark matter sector and try to parameterize it in simple 
terms. 

We have not addressed how to extend the simplified models to 
realistic and viable models which are consistent with the symme¬ 
tries of fhe Sfandard Model. Simplified models offen violafe gauge 
invariance which is a crucial principle for building a consisfenf 
BSM model which incorporafes SM fogefher wifh new physics. For 
example, wifh a new heavy gauge vector boson mediafing DM in- 
feracfions, one needs nof jusf fhe dark maffer and ifs mediator, buf 
also a mechanism which provides mass fo fhis mediator in a gauge 
invarianf way. 

Considering bofh fhe simplified model and ofher elemenfs nec¬ 
essary for a consisfenf fheory is a nexf logical step. The aufhors 
of [Bel+i5c] ferm fhese Minimal Consisfenf Dark Maffer (MCDM) 
models. MCDM models are af fhe same fime sfill foy models fhaf 
can be easily incorporafed info a bigger BSM model and explored 
via complemenfary consfrainfs from collider and direcf/indirecf 
DM search experimenfs as well as relic densify consfrainfs. 

The idea of an inerf Two-Higgs Doublef Model (IDM) was in- 
froduced more fhan 30 years ago in Ref [DM78]. The IDM was firsf 
proposed as a Dark Maffer model in Ref. [BHR06] and ifs phe¬ 
nomenology furfher sfudied in Refs. [LH+07; Ham+09; LHYii; 

Gus+07; DS09; ATLi4d; ADK09; ATS09; NTV09; GCI13; GHSi3b; 

Bel+i5c]. If is an extension of fhe SM wifh a second scalar doublef 
(p2 wifh no direcf coupling fo fermions. This doublef has a discrete 
Z2 symmefry, under which ^2 is odd and all fhe ofher fields are 
even. The Lagrangian of fhe odd sector is, 

C=liDf,cp2f-V{cp,,cl>2) (A.14) 

wifh fhe pofenfial V confaining mass terms and cpi — cj)2 inferac- 
fions: 


V = -m^((plcpi)-m^(cp^cp 2 )+Ai(^l^i)^ + A 2 (cp^cp 2 f 

+ A3{cj)2(p2)i4’l^l) + A4{(p2(pl)i^l^2) + ^ + (</’2^l) 


where (pi and (p2 are SM and inerf Fliggs doublefs respecfively 
carrying fhe same hypercharge. These doublefs can be parameter¬ 
ized as 




^2 


1 / V 2 h+ \ 
\/2 \ hi + ih2J 


(A. 16) 


In addifion fo fhe SM, fhe IDM infroduces four more degrees 
of freedom coming from fhe inerf doublef in fhe form of a Z2-odd 
charged scalar and fwo neufral Z2-odd scalars hi and h2. The 
lighfesf neufral scalar, hi is idenfified as fhe dark maffer candidate. 
Aspecfs of fhe IDM collider phenomenology have been sfudied in 
[BPVoi; AHT08; Arh+14; Bel+i5c; BHR06; LGE09; CMR07; Dol+10; 


(A.15) 
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MSTio; Gus+12; ABG12; SK13; GHSi3a; Bel+i5a]. Its LHC signa¬ 
tures include dileptons [Dol+10; Bel+i5a], trileptons [MSTio] and 
multileptons [Gus+12] along with missing transverse energy, mod¬ 
ifications of fhe Higgs branching rafios [ABG12; SK13; GHSi3b], as 
well as £t -h jef, Z, and Higgs and Bj + VBF signals (see Figs. A.3- 
A.8). 



9 <siSLSiSLi!rs~is~s~s~s~^ a 


.'hi 


'-hi 

,''hi 

'-hi 



Figure A.3: Feynman diagrams for 
SS ^ ^ihi+ g process contributing 
to mono-jet signature, adapted from 
[Bel+i5c]. 



Figure A.4: Feynman diagrams for 
qq hih2 + g {gq ^ hih2 + q) process 
contributing to mono-jet signature, 
adapted from [Bel+i5c]. 
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a 



Figure A.5: Feynman diagrams for 


qq —> h\hi + Z process contributing 

i ^2 ,'h-[ 

to mono-Z signature, adapted from 


[Bel+i5c]. 


Based on fhe various LHC search charmels, DM phenomenology 
issues and fheorefical considerafions, numerous works have pro¬ 
posed benchmark scenarios for fhe IDM, see e.g. [Gus+12; GHSi3b] 
while a FeynRules implemenfafion (including MadGraph, CalcHEP 
and micrOMEGAs model files) was provided in [Gus+12]. An up- 
dafed analysis of fhe paramefer space has recenfly been performed 
in Ref. [Bel+i5c]. 

The aufhors suggesfed fo sfudy mono-X signafures fhaf are rel- 
evanf fo model-independenf collider DM searches, and evaluafed 
fheir rafes presenfed below. They have implemenfed and cross¬ 
checked fhe IDM model info CalcHEP and micrOMEGAs, wifh 
an implemenfafion publicly available on fhe HEPMDB dafabase, 
including loop-induced HHG and 77H models. They propose an 
addifional sef of benchmark poinfs, mosfly inspired by mono-X 
and VBE searches (Table. A.i). Though fhe overall paramefer space 
of IDM is 5-dimensional, once all relavanf consfrainfs are applied 
fhe paramefer space relevanf fo a specific LHC signafure fypically 
reduces fo 1-2 dimensional. In fhe mono-jef case, one can use fwo 
separafe simplified models, a yy —> hqhi + g process (via Higgs 
mediator) and aqq ^ hih2 + g{gq hih2 + q) process (fhrough 
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Figure A.6: Feynman diagrams for 
gS ^ H process contributing to 

mono-Higgs signature, adapted from 
[Bel+i5c]. 
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Figure A.y. Feynman diagrams for 
qq ^ hih2 + H process contributing to 
mono-Higgs signature, adapted from 
[Bel+i5c]. 


a Z-boson mediator) to capture the physics relevant to the search. 
The cross sections for the various mono-X and VBF signatures pro¬ 
duced by this model are displayed in Fig. A.9. 



Figure A.8: Diagrams for qq qqhih^ 
DM production in vector boson fusion 
process, adapted from [Bel+i5c]. 
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Figure A.g: LHC cross section at 
13 TeV for various signatures, from 
[Bel+i5c]. 
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B 

Appendix: Presentation of experimental results for rein¬ 
terpretation 


When collider searches present results with the recommended 
benchmarks, we suggest the following: 

• Provide limits in collider language, on fundamental parame¬ 
ters of the interaction: the couplings and masses of particles in 
simplified model. 

• Translate limits to non-collider language, for a range of assump- 
fions, in order fo convey a rough idea of the range of possibili- 
fies. The details of this point are left for work beyond the scope 
of fhis Forum. 

• Provide all necessary maferial for fheorisfs to reinterpret simpli¬ 
fied model results as building blocks for more complefe models 
(e.g. signal outflows, accepfances, etc). This point is detailed 
further in this appendix. 

• Provide model-independent results in terms of limifs on cross- 
section times efficiency fimes accepfance of new phenomena for 
all cases, buf especially when EFTs are employed as benchmarks 
This recommendation has been issued before: see Ref. [Kra+12] 
for detailed suggestions. 

• Provide easily usable and clearly labeled results in a digitized 
format, e.g. [Hep] entries, ROOT histograms and macros or 
tables available on analysis public pages. 

This appendix describes further considerations for reinterprefa- 
fion and reimplementafion of fhe analyses, as well as for the use of 
simplified model results directly given by the collaborations. 

B.i Reinterpretation of analyses 

In the case of reinterprefation for models differenf than those pro¬ 
vided by the experimental collaborations, the information needed 
primarily includes expected and observed exclusion lines along 
with their ztlu uncertainty, expected and observed upper limits in 
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case of simplified models, efficiency maps and kinematic disfribu- 
tions as reporfed in the analysis. If fhe kinematics of fhe new model 
fo be fesfed in fhe reinferprefafion is similar fo fhat of fhe original 
model provided by fhe collaborafion, if will be sfraighf-forward fo 
rescale fhe resulfs provided fo mafch fhe new model cross-secfion 
using this information. 

B.2 Reimplementation of analyses 

One of fhe imporfanf developmenfs in recenf years is an acfive de¬ 
velopment of software codes [Dum+15; Con+14; Kim+i5b; CYii; 
Kim+i5a; Bar+14] necessary for recasfing analyses. The aim of 
fhese codes is fo provide a public library of LHC analyses fhaf have 
been reimplemenfed and validafed, offen by fhe collaborafions 
fhemselves. Such libraries can fhen be used fo analyze validify of a 
BSM scenario in a sysfemafic and effective marmer. The availabilify 
of public libraries furfher facilifafes a unified framework and can 
lead fo an organized and cenfral sfrucfure fo preserve LHC infor¬ 
mation long ferm. The reimplemenfafion of an analysis consisfs of 
several sfages. Typically, fhe analysis nofe is used as a basis for fhe 
implemenfafion of fhe preselection and even! selection cufs in fhe 
user analysis code wifhin fhe recasfing frameworks. Signal evenfs 
are generafed, and passed fhrough a parameferized defecfor sim¬ 
ulation using software such as Delphes or PCS [Fav+14; Pgs]. The 
reconstructed objects are then analyzed using the code written in 
the previous step, and the results in terms of number of events are 
passed through a statistical analysis framework to compare with 
the backgrounds provided by the collaborations. 

In order to be able to effectively use such codes, it is important to 
get a complete set of information from the collaborations. 

For what concerns the generation of the models, it is desirable to 
have the following items as used by the collaborations: 

• Monte Carlo generators: Monte Carlo generators along with the 
exact versions used to produce the event files should be listed. 

• Production cross sections: The order of production cross sections 
(e.g. LO,NLO,NLL) as well as the codes which were used to 
compute them should be provided. Tables of reference cross 
sections for several values of particle masses are useful as well. 

• Process Generation: Details of the generated process, detailing 
number of additional partons generated. 

• LHE files: selected LHE files (detailing at least a few events 
if not the entire file) corresponding to the benchmarks listed 
in the analysis could also be made available in order to cross 
check process generation. Experimental collaborations may gen¬ 
erate events on-the-fly without saving the intermediate LHE 
file; we advocate that the cross-check of process generation is 
straight-forward if this information is present, so we encourage 
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the generation of a few selecfed benchmark poinfs allowing for 
a LHE file fo be saved. Special attention should be paid to list 
the parameters which change the production cross section or 
kinematics of fhe process e.g. mixing angles. 

• Process cards: Process cards including PDF choices, defails of 
matching algorithms and scales and details of process genera¬ 
tion. If process cards are nof available, fhe above ifems should be 
clearly identified. 

• Model files: For models which are nof already implemenfed 
in MadGraph5_aMC@NLO, fhe availabilify of fhe corre¬ 
sponding model files in fhe UFO formaf [Deg+12] is highly 
desired. This formaf defails fhe exacf nofation used in fhe 
model and hence sefs up a complefe framework. In case Mad- 
Graph5_aMC@NLO is nof used, enough information should 
be provided in order fo clearly identify fhe underlying model 
used for inferprefafions and reproduce fhe generation. 

The ATLAS/GMS Dark Matter Forum provides most of fhe infor¬ 
mation needed within its SVN repository [Fork] and on a dedicated 
HEPData [Hep] page dedicated to the results in this report. 

Efficiency maps and relevanf kinemafic disfribufions as reporfed 
in fhe analysis should be provided, in a digifized formaf wifh 
clearly specified unifs. If selecfion criferia cannof be easily simu- 
lafed fhrough parameferized defecfor simulation, fhe collaborations 
should provide fhe efficiency of such cufs. Overall reconsfrucfion 
and identification efficiencies of physics objecfs are given as an in- 
puf fo fhe defecfor simulation software. If is fhus very useful fo gef 
paramefrized efficiencies for reconsfrucfed objecfs (as a funcfion of 
fhe rapidify t] and/or fransverse momenfum pj), along wifh fhe 
working poinfs af which fhey were evaluafed (e.g. loose, fighf se¬ 
lecfion). Objecf definifions should be clearly identifiable. Digifized 
kinemafic disfribufions are often necessary for fhe validation of fhe 
analysis so fhaf fhe resulfs from fhe collaborafion are obfained, and 
so are fables containing the events passing each of fhe cufs. 

The availabilify of digifized dafa and backgrounds is one of fhe 
primary requiremenfs for fasf and efficient recasting. Platforms 
such as HepDafa [Hep] can be used as a cenfralized reposifory; 
alfernafively, analysis public pages and fables can be used for dis¬ 
semination of resulfs. Bofh dafa and Sfandard Model backgrounds 
should be provided in fhe form of birmed hisfogram fhaf can be 
inferpolafed if needed. 

A defailed descripfion of fhe likelihood used in order fo derive 
fhe limifs from fhe comparison of dafa fo signal plus background 
should be given. This can be inferred from fhe analysis documenfa- 
fion ifself, however direcf availabilify of fhe limif seffing code as a 
workspace in RooSfafs or HisfFitfer [Baa+15] is highly desirable. 

Finally, fhe collaborations can also provide an analysis code 
directly implemented in one of fhe public recasfing codes defailed 
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above. Such codes can be published via INSPIRE [Ins] in order to 
track versioning and citations. 

B.3 Simplified model interpretations 

Dark Matter searches at the LHC will include simplified model 
interpretations in their search results. These interpretations are 
simple and can be used for a survey of viabilify of paramefer space. 
Codes such as [Kra+i4a; Kra+i4b; Pap+14] can make use of fhe 
simplified model resulfs given in fhe form of 95% Confidence Level 
(CLs) upper limif or efficiency maps in order fo fesf Beyond fhe 
Sfandard Model paramefer space. As menfioned above, if will fhus 
be exfremely useful if fhe resulfs are given in a digifized form fhaf 
is easily usable by fhe fheory communify. 

The paramefer space of fhese models should be clearly specified. 
For example, for a simplified model confaining dark matter mass 
m^, mediator mass M^ed ^rid couplings g^, if will be very useful 
fo have upper limifs on fhe producf of couplings ^/gxgq or cross 
secfion times branching rafio as a funcfion of m^, Mmed- Limifs 
on visible cross sections of fhe simplified models considered for 
inferprefafions should be made available. 

The usage of simplified model resulfs relies on inferpolafing be- 
fween upper limif values. In order fo facilifafe fhe inferpolafion, 
regions where large variation of upper limifs is observed should 
confain denser grid, if a uniform grid over fhe entire plane is nof 
possible. For simplified model involving more fhan fhree parame- 
fers (two masses and product of couplings), slices of upper limifs in 
fhe additional dimensions will be necessary for reinferprefafion. 

As already menfioned in fhe infroducfion fo fhis Chapter, accep- 
fance and efficiency maps for all fhe signal regions involved in fhe 
analysis should be made available. These resulfs are nof only useful 
for model fesfing using simplified models buf also fo validate im- 
plemenfafion of fhe analysis. Information abouf fhe mosf sensitive 
signal regions as a funcfion of new parficle masses is also useful in 
order fo defermine fhe validify of approximafe limif setting proce¬ 
dures commonly used by fheorisfs. 



c 

Appendix: Additional details and studies within the 
Forum 


Further information for baryonic 7 J Model 
Cross-section scaling 

The dependence of the cross section of the pp ^ Hxx + ^ pro¬ 
cess on ghz'Z' is shown in Figure C.i. The curves have been fit to 
second-order pol5momials, where y is the cross-section and x is the 
coupling ghz’Z'- 

For = 100 GeV, the fit function is 

y = - 0.12 - 3.4 X lO-^x -H 2.7 x 10 ” V 
. For = 1 TeV, the fit function is is 

y = 0.0012 - 2.4 x 10 "^x -h 1.5 x 10 "V 

/ 

y = - 0.12 - 3.4 X 10 "^x -F 2.7 x 10 " V. (C.i) 

For Mined = 1 TeV, the fit function is is: 


y = 0.0012 - 2.4 x lO’^x -F 1.5 x 10 ” V. 


(C.2) 


o [pb] o [pb] 
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9.zz.tGeV] 


Figure C.i: Cross section of the pp —> 
Hxx process as a function of ghz'Z' fot 
m2' = 100 GeV (left) and m^' = 1 TeV 
(right). The fit functions are shown in 
the text. 



g,,z,[GeV] 
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